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Chapter  1 


INTRODUCTION 

Air  Traffic  Control  (ATC)  is  a broad  concept  designed 
to  provide  information  and  assistance  to  insure  the  safe, 
orderly  flow  of  aircraft.  It  is  a concept  rooted  in  a com- 
plex system  of  flight  rules,  regulations,  and  procedures 
which  becomes  tangible  in  the  form  of  a multitude  of 
facilities  using  varying  degrees  of  communications,  radar 
and  computer  technology.  Positive  ATC,  however,  is  achieved 
through  the  interaction  of  people;  in  particular,  through 
the  interaction  of  pilots  and  air  traffic  controllers. 

To  accomplish  his  job,  the  air  traffic  controller  is 
placed  at  the  apex  of  an  information  system.  His  information 
sources  are  many;  radio  communications,  radar  presentations, 
flight  plans,  flight  progress  strips,  and  telephone  communi- 
cations. The  information  he  receives  describes  the  dynamic 
air  traffic  environment  at  any  given  instant  in  terms  of 
numbers,  locations,  directions,  velocities  and  altitudes 
of  aircraft.  The  controller's  responsibilities  are  to 
continuously  assess  the  information,  to  recognize  potential 
conflicts,  and  to  formulate  and  implement  corrective  action 
in  a timely  manner.  As  the  density  of  air  traffic  increases, 
a corresponding  increase  in  information  input  to  the  con- 
troller and  in  the  probability  of  traffic  conflict  occurs. 


Simultaneously,  the  time  available  for  the  controller  to 
evaluate  data,  to  recognize  potential  conflicts,  and  to 
formulate  corrective  action  decreases. 

Safety  dictates  that  he  concentrate  his  attention 
on  the  resolution  of  potential  conflict.  To  achieve  this 
end,  he  will  take  appropriate  corrective  action.  In  the 
context  of  conflict  resolution,  a controller's  corrective 
action  can  be  characterized  as  preventative  in  nature  and 
may  impose  an  economic  penalty  in  terms  of  the  costs 
incurred  by  the  various  delaying  techniques  employed  to 
preclude  conflict. 

Economic  considerations  are  less  important  than 
safety  considerations,  but  could  both  be  optimally  dealt 
with  simultaneously,  regardless  of  the  air  traffic  situation? 
In  answering  this  question  two  important  elements  must  be 
considered:  The  amount  of  time  the  controller  has  for  data 

analysis,  and  the  amount  of  information  available  to 
develop  a plan  for  optimum  traffic  flow. 

Statement  of  the  Problem 

The  problem  is  twofold.  One,  the  air  traffic  con- 
troller does  not  always  have  sufficient  time  to  develop  a 
plan  to  maximize  both  safety  and  economy;  and  two,  he  does 
not  have  adequate  information  to  formulate  such  a plan. 

Improvements  in  the  design  of  the  air  route  system 
could  provide  the  controller  with  additional  time  and 
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increase  his  effectiveness  by  simplifying  and  streamlining 
the  flow  of  traffic.  The  greatest  potential  for  design 
improvements  lies  in  the  terminal,  or  approach  control 
area  (ACA) . It  is  within  the  ACA  that  arrival  routes, 
after  having  separated  from  the  enroute  network  of  air 
routes,  converge  on  the  destination  airport  and  that 
departure  routes  diverge  from  the  airport  to  link  up  with 
the  enroute  structure.  In  the  ACA,  arrival  and  departure 
routes  are  integrated  into  a single  network  of  flight  paths 
normally  characterized  by  progressive  altitude  restrictions. 
Aircraft  flight  separation  requirements,  both  vertical  and 
horizontal,  are  driving  forces  in  designing  the  ACA  and  in 
developing  procedures  to  handle  air  traffic  in  the  ACA. 

Significant  changes  to  the  ACA  cannot  be  implemented 
on  an  experimental  basis,  however,  without  running  the  risk 
of  disrupting  operations.  Furthermore,  changes  can  be  quite 
costly  and  time  consuming  (17:1),  depending  upon  the  extent 
of  physical  re-design  and  upon  the  amount  of  coordination 
among  air  traffic  agencies  required  by  the  Federal  Aviation 
Administration  (FAA) . 

An  alternative  to  actually  experimenting  with  an  ACA 
in  seeking  design  improvements  is  to  use  computer  assisted 
simulation.  Tuan  defines  a computer-simulation  model  as 

...  a procedural-logical-mathematical  represen- 
tation of  a real-world  system  programmed  for  digital 
computers  within  which  experiments  can  be  conducted 
over  specific  periods  of  time  [17:1]. 
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Computer-simulation  is  an  ideal  instrument  to  experiment 
with  a variety  of  changes  to  the  ACA,  and  to  evaluate  their 
effect,  without  creating  actual  traffic  disruption  and  with- 
out incurring  excessive  costs.  Potential  changes  to  the  ACA 
that  can  be  explored  through  simulation  include:  (1)  Changes 

in  the  routing  of  flight  paths;  (2)  Changes  in  the  allocation 
of  controlled  airspace;  (3)  Changes  in  the  location  of  flight 
facilities;  (4)  Changes  in  procedures  and  regulations  to 
handle  traffic;  and  (5)  Changes  in  the  decision  rules  used 
by  controllers  to  sequence  air  traffic  (17:2).  Those  changes 
that  prove  to  be  of  benefit  through  simulation  experiments 
could  then  be  considered  for  actual  implementation  in  the  ACA. 

Besides  assisting  in  the  evaluation  of  changes  in  the 
ACA  which  provide  the  controller  more  time,  simulation  could 
be  used  to  provide  the  controller  additional  information. 

The  effect  of  various  controller  actions  on  simulated  traffic, 
representative  of  a controller's  real  traffic  at  any  given 
time,  could  be  evaluated  on  a near  real-time  basis  (17:5). 

In  this  respect  the  controller  would  have  assistance  in 
planning  since  he  could  evaluate  a course  of  action  prior 
to  its  implementation.  The  planning  assistance,  hopefully, 
would  result  in  an  optimal  flow  of  air  traffic  that  not 
only  considers  flight  safety  but  ecomony  as  well. 


4 


Literature  Review 


The  application  of  computer  technology  to  the  ter- 
minal ATC  environment  is  not  a new  concept.  A report  by 
the  Mitre  Corporation  presented  four  articles  (15:9-14) 
which  summarized  the  present  state  of  automation  of  terminal 
and  enroute  ATC  facilities.  Technical  features  of  the 
Automated  Radar  Terminal  System  (ARTS)  were  described  in 
some  detail.  The  ARTS  provides  alphanumeric  data  blocks 
(containing  aircraft  identity,  altitude,  and  ground  speed) 
which  are  superimposed  on  the  controllers  radar  display 
(15:1).  The  development  of  the  ARTS  has  provided  the  con- 
troller with  an  improved  picture  of  the  existing  traffic 
situation,  and  allows  him  to  monitor  aircraft  flight  progress 
to  insure  compliance  with  control  instructions. 

The  Mitre  report  also  discussed  the  Radar  Data  Pro- 
cessing (RDP)  System,  Conflict  Alert,  and  Intermittent 
Positive  Control,  all  of  which  have  application  to  the  ter- 
minal area.  RDP  is  directed  towards  improvement  of  the 
accuracy  and  reliability  of  radar  surveillance  of  aircraft. 
Conflict  Alert  and  Intermittent  Positive  Control  are  both 
being  developed  to  prevent  midair  collisions  in  congested 
areas  (15:19-26).  The  procedure  involves  the  continuous 
surveillance  of  a block  of  airspace  projected  in  front  of 
each  aircraft.  If  the  projected  airspace  of  two  or  more 
aircraft  intersect,  the  computer  provides  an  alert  signal 
to  the  controller. 
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The  articles  mentioned  above  were  concerned  with 


functionally  oriented  technological  capabilities.  This 
technology  provides  the  basic  building  blocks  onto  which 
can  be  added  control  aiding  functions.  Holland  and  Garceau 
presented  a summary  of  selected  efforts  in  the  area  of  con- 
trol aiding  functions  (7:3-1  to  3-123).  These  efforts 
represent  an  evolutionary  development  which  has  culminated 
in  recent  efforts  to  completely  automate  the  sequencing 
and  spacing  of  aircraft.  The  more  important  of  these 
projects  are:  Final  Approach  Spacing  for  ARTS  (FASA) , 

and  Computer  Aided  Approach  System  (CAAS) . 

The  FASA  project  involved  the  preliminary  sequencing 
of  aircraft  based  on  predicted  arrival  times  at  the  runway. 

The  sequence  is  modified,  if  necessary,  prior  to  the  Final 

★ 

Approach  Fix  (FAF)  and  from  that  point  the  sequence  is  firm 
(7:3-94  to  3-109).  Although  the  method  employed  in  this 
project  appears  to  have  potential,  a field  test  was  incon- 
clusive because  of  difficulties  in  implementation. 

The  CAAS  project  attempted  to  achieve  final  approach 
spacing  by  using  a computer  to  specify  departure  times  from 
navigational  fixes  outside  the  ACA.  Departure  times  were 
then  "made  qood"  by  enroute  or  transition  controllers.  A 
field  test  indicated  the  CAAS  system  provided  more  con- 
sistent and  accurate  landing  intervals  than  the  manual 
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control  system.  Other  potential  benefits  were  not  com- 
clusively  identified  due  to  constraints  imposed  by  the  test 
environment.  In  addition,  the  system  was  not  favorably 
received  by  the  controllers  since  it  caused  an  increase  in 
workload  (7:3-109  to  3-123). 

The  use  of  simulation  as  a means  to  evaluate  alter- 
native control  methods  was  presented  by  Gabrielli.  In  his 
study,  Gabrielli  experimented  with  approach  geometries 
(geometric  arrangements  of  flight  paths  leading  to  a runway) . 
Two  approach  geometries  were  compared  using  the  statistical 
output  from  a computer  simulation.  The  comparison  was  used 
to  support  his  recommendation  of  one  geometry  over  the  other 
(6:1-21)  . 

A similar  method  was  used  by  Mohleji  and  Horowitz  in 
their  analysis  of  Denver's  Stapleton  Airport.  Their  study 
compared  various  approach  geometry  configurations  via  com- 
puter simulation  and  recommended  an  optimal  geometry  for 
use  in  conjunction  with  the  ARTS  system.  The  recommendation 
was  not  adopted,  however,  because  of  the  heavy  reliance  on 
visual  separation  of  aircraft  in  the  Stapleton  ACA.  Visual 
separation  allows  reduce  intervals  between  aircraft  and  last 
minute  pilot  initiated  corrections  near  final  approach.  Con- 
sequently, the  automated  control  system  could  not  exceed  the 
effectiveness  of  the  manual  control  system  by  a significant 
margin  (10:1-1  to  5-4;  9:5-1  to  5-8). 
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The  FASA,  CAAS , and  approach  geometry  studies  pro- 
vide evidence  that  simulation  can  be  a useful  tool  in  the 
analysis  of  traffic  flow.  A synthesis  of  the  approaches 
taken  in  these  studies  has  provided  the  basis  for  this 
research.  In  addition  the  Holland  and  Garceau  report  pro- 
vides a listing  of  control  tools  that  can  be  incorporated 
into  the  development  of  an  ATC  simulation  model  (7:2-1  to 
2-18)  . 


Justification 


In  the  summer  of  1968,  the  Department  of  Transpor- 
tation Air  Traffic  Control  Advisory  Committee  was  formed 
".  . . for  the  purpose  of  recommending  an  air  traffic 
control  (ATC)  system  for  the  1980's  and  beyond  [19:3]." 

The  committee's  technical  staff  was  composed  of  approxi- 
mately 150  people  drawn  from  all  segments  of  the  aviation 
industry  and 


. . . concentrated  on  control  of  aircraft  through 
the  airspace,  from  takeoff  to  landirg.  Emphasis  was 
placed  on  the  denser  portions  of  the  airspace  where 
the  danger  of  midair  collisions  and  the  need  for 
efficient  use  of  scarce  resources  (principally  runways 
and  terminal  airspace)  make  sophisticated  ATC  mandatory 
if  safety  is  to  be  assured  without  sacrifice  of  capacity 
and  without  unacceptable  delays  or  interference  with 
freedom  of  flight  [19:3], 

The  report  of  the  committee  pointed  out  air  traffic  was 
experiencing  severe  problems  and  noted  the  crisis  at  a few 
high  volume  airports  was  due  to  the  failure  of  airport 
capacity  and  of  air  traffic  control  capacity  to  keep  pace 


8 


with  the  growth  of  the  aviation  industry.  The  committee 
further  noted  the  demand  for  ATC  service  was  estimated  to 
almost  triple  by  1980  and  to  triple  again  by  1995  (19:5). 

Simulation  provides  a logical  and  economic  means  to 
experimentally  test  new  concepts  designed  to  cope  with 
increasing  demand  for  ATC  services.  In  addition  it  is 
possible  to  test  ATC  designs  for  efficiency  at  present 
levels  of  traffic  and  at  the  levels  predicted  to  occur  in 
1980  or  1995. 

There  is  additional  justification  for  this  research 
in  the  potential  economy  resulting  from  improved  traffic 
handling.  Simulation-based  analysis  of  the  air  traffic 
controller's  environment  provides  a means  to  reduce  air 
transportation  costs.  The  inherent  assumption  in  this 
approach  is  that  in  any  given  air  traffic  situation  there 
exists  an  optimal  method  to  handle  the  traffic  which  will 
keep  in-flight  time,  flight  distance,  and  fuel  consumption 
at  minimum  levels.  In  additior,  the  efficient  handling 
of  aircraft  will  effectively  reduce  the  workload  imposed 
on  the  controller.  These  savings,  in  particular  the  reduc- 
tion of  fuel  consumption,  are  of  direct  interest  to  the 
Department  of  Defense  (1:39-41;  20:7-10). 

Objective 

The  objective  of  this  research  was  to  apply  com- 
puter simulation  technology  to  a specific  military  ACA  as 
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a means  to  explore  and  to  evaluate  potential  improvements  in 
ACA  design  or  ATC  procedures.  This  objective  was  to  be 
pursued  in  three  sequential  phases. 

The  first  phase  called  for  the  development  of  a simu- 
lation model.  The  model,  once  validated,  would  then  become 
the  basic  tool  to  be  used  in  each  of  the  remaining  three 
phases.  In  order  to  be  a useful  tool,  the  model  had  to  dupli- 
cate an  actual  ACA  as  close  as  possible  in  both  the  physical 
design  features  of  the  ACA  and  in  the  flow  of  air  traffic 
within  the  ACA.  Once  the  model  was  constructed  and  validated, 
phase  two  could  commence. 

Phase  two  required  that  the  validated  model  be  used 
in  an  experimentation  process  wherein  potential  changes  in 
design  or  procedure  could  be  tested.  Experimentation  would 
involve  altering  the  validated  model,  to  reflect  potential 
improvements,  and  comparing  the  output  of  the  altered  model 
against  that  of  the  validated  model. 

Following  experimentation,  phase  three  dictated 
the  integration  of  the  validated  model  with  a computer 
algorithm  to  identify  potential  air  traffic  congestion  or 
conflict.  The  idea  underlying  phase  three  was  to  use  simu- 
lation to  evaluate  various  controller  actions  that  could  be 
implemented  to  preclude  impending  congestion  or  conflict. 

Such  an  evaluation  procedure  could  be  used  to  test  the 
feasibility  of  providing  the  controller  with  near  real-time 
information  to  assist  him  in  conflict  resolution. 
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Research  Questions 


This  research  sought  the  answer  to  the  following 
four  questions:  (1)  Can  the  variables  comprising  the 

military  ATC  environment  be  adequately  represented  by  a 
computer  model?  (2)  Can  changes  in  the  terminal  environ- 
ment be  effectively  evaluated  by  the  simulation  model? 

(3)  Can  the  model  be  used  to  predict  future  congestion  or 
conflict?  (4)  Can  the  model  be  used  to  compare  available 
solutions  to  a predicted  conflict  situation  and  select  a 
course  of  action  to  resolve  the  conflict? 
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Chapter  2 
BACKGROUND 

Thie  research  specifically  focuses  on  a single  ATC 
facility,  the  Dayton  Approach  Control,  at  Wr ight-Patterson 
AFB , Ohio.  A computer  simulation  model  was  developed  to 
represent  the  air  traffic  environment  within  the  Dayton  ACA. 
Simulated  air  traffic,  having  the  same  frequency,  volume  and 
mixture  characteristics  as  actual  air  traffic  in  the  Dayton 
ACA,  is  generated  by  the  model  and  flows  through  the  simula- 
tion. The  simulation  model  considers  the  interaction  of  all 

\ 

air  traffic  but  is  limited  to  the  sequencing  of  aircraft  to 
and  from  the  runway  at  Wright-Patterson  AFB;  the  model  has 
not  been  expanded  to  include  other  airfields  located  within 
the  Dayton  ACA. 

The  Dayton  ACA 

The  geographical  (horizontal)  boundaries  of  the  Dayton 
ACA  are  depicted  in  Figure  1.  Within  these  boundaries,  Dayton 
Approach  Control  has  jurisdiction  over  the  airspace  from  the 
surface  up  to  9,000  feet  except  in  the  vicinity  of  the  Rose- 
wood and  the  Richmond  navigation  facilities  (i.e.  VORTACs) . 
Over  the  Rosewood  and  the  Richmond  VORTACs,  Dayton  Approach 
Control's  airspace  extends  from  the  surface  upward  to  5,000 

) 

feet  and  to  4,000  feet,  respectively. 
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Figure  1.  Dayton  Approach  Control 


There  are  12  airports  located  in  the  Dayton  ACA.  The 
two  major  airports  are  Wright-Patterson  AFB  and  James  M.  Cox — 
Dayton  Municipal  Airport.  Eleven  low  altitude  airways  pene- 
trate the  Dayton  ACA  and  six  electronic  navigational  aids, 
used  to  define  airways  and  instrument  approaches,  are  con- 
tained within  the  boundaries  of  the  Dayton  ACA.  Aircraft 
navigate  within  the  ACA  by  use  of  the  navigational  aids,  by 
radar  vectors  assigned  by  Dayton  Approach  Control  or  by 
visual  references.  Aircraft  operate  in  accordance  with 
Instrument  Flight  Rules  (IFR)  or  Visual  Flight  Rulqs  (VFR) . 

The  air  traffic  in  the  Dayton  ACA  is  a mixture  of  military, 
commercial,  and  civilian  aircraft.  The  Dayton  ACA  recorded 
approximately  237,500  air  operations*  in  Fiscal  Year  1975, 
and  approximately  227,300  air  operations  in  Fiscal  Year 
1976  (3).  The  percentage  of  air  operations  in  the  Dayton 
ACA  generated  by  Wright-Patterson  AFB  has  been  calculated 
by  the  statistical  analysis  of  air  traffic  data  (see  Data 
Collection) . 

Wright-Patterson  AFB 

Wright-Patterson  AFB  has  a single  runway  oriented 
to  a magnetic  bearing  of  050/230  Degrees.  The  runway  is 
12,600  feet  long  by  300  feet  wide  and  is  adequate  to  handle 


The  phrase  "Air  Operation"  includes  such  activities  as: 
takeoff,  landing,  through-flight,  execution  of  an  instrument 
approach  procedure,  and  execution  of  a missed  approach  pro- 
cedure. For  a more  technical  definition,  see  FAA  Regulation 
7210.3,  pages  179-181. 
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the  largest  aircraft.  Ten  approach  procedures  are  available 
for  navigation  to  the  runway.  These  include  both  low  and 
high  altitude  non-precision  instrument  approach  procedures, 
Approach  Surveillance  Radar  (ASR) . Precision  Approach 
Radar  (PAR) , and  Instrument  Landing  System  (ILS)  approaches. 
Figure  2 is  a diagram  of  the  air  field  at  Wright-Patterson  AFB. 

Computer  Language 

The  simulation  programming  language  used  in  the 
development  of  the  simulation  model  of  the  Dayton  ACA  is 
the  General  Purpose  Simulator  System  (GPSS) . GPSS  is 
both  a computer  language  and  a computer  program.  As  a lan- 
guage it  provides  the  capability  to  describe  a system  with 
a well  defined  vocabulary  and  grammar.  With  GPSS  a discrete- 
event  simulation  computer  model  of  a system  can  be  built  to 
reproduce  the  dynamic  behavior  of  the  system  over  time  (14: 
vii) . As  a computer  program,  GPSS  interprets  the  simulation 
model  of  the  system  described  in  the  GPSS  language,  thereby 
making  it  possible  to  conduct  experiments  with  the  simulation 
model  on  a computer.  The  GPSS  program  schedules  each  event 
for  each  traffic  element  in  the  system,  to  occur  in  future 
simulated  time  by  reference  to  a self-contained  simulation 
clock  (8:2-1).  It  also  causes  the  events  to  occur  in  the 
proper  time-ordered  sequence  and  provides  a means  of  assign- 
ing relative  priorities  to  be  used  in  resolving  time  con- 
flicts. 
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Definitions 


The  following  list  of  definitions  is  provided  for 
clarification  of  terms  used  in  this  paper. 

Airspeed : The  velocity  of  an  aircraft  through  the  air  mass. 

In  the  context  of  this  paper,  airspeed  is  considered  to  be 
the  same  as  groundspeed;  that  is,  the  aircraft's  speed  over 
the  ground.  (1)  Cruise  Airspeed  - the  airspeed  maintained 
in  level  flight  at  an  enroute  altitude.  (2)  Maneuvering 
Airspeed  - the  airspeed  used  within  the  Approach  Concrol 
Area  to  maneuver  the  aircraft  from  an  Arrival  Fix  to  tho 
Final  Approach  Fix.  (3)  Final  Approach  Airspeed  - the 
airspeed  used  during  final  descent  for  landing. 

Alternate  Vectoring  Route  (AVR) : Any  route  other  than  a 

Standard  Vectoring  Route. 

Approach  Control  Area  (ACA) : The  airspace  under  the  control 

of  personnel  at  the  Radar  Approach  Control  Facility  (includ- 
ing tower  controllers) . The  horizontal  (geographic)  and 
vertical  (altitude)  boundaries  for  the  Dayton  ACA  are  pre- 
sented at  the  beginning  of  this  chapter. 

Approach  Geometry:  A designated  route  for  vectoring  air- 

craft from  an  Arrival  Fix  to  a Final  Approach  Fix. 

Arrival  Fix  (AF) : A point  in  close  proximity  to  the  perim- 
eter of  the  Approach  Control  Area.  Aircraft  inbound  to  the 
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Dayton  ACA  are  transitioned  from  the  enroute  structure  to  the 
AFs.  During  the  transitioning  process,  control  of  an  air- 
craft is  transferred  from  the  enroute  controller  to  the 
Approach  controller. 

Conflict : A situation  in  which  two  or  more  aircraft  will 

arrive  at  the  same  point  at  the  same  time  without  required 
vertical  or  horizontal  separation. 

Congestion:  A situation  wherein  the  rate  of  arriving  air- 

craft is  such  that  one  or  more  aircraft  must  be  delayed  in 
order  to  prevent  conflict. 

Final  Approach  Fix  (FA F) : The  point  from  which  an  aircraft 

initiates  its  final  descent  for  landing.  An  aircraft  is 
maneuvered  either  by  a Standard  Vectoring  Route  or  an  Instru- 
ment Approach  Procedure  so  as  to  arrive  at  the  FAF  at  the 
required  altitude,  airspeed,  and  heading. 

Flight  Information  Publications  (FLIP) : A variety  of  publi- 

cations for  use  by  pilots  and  navigators  in  planning  and 
executing  a flight.  FLIP  is  provided  by  the  Department  of 
Defense  (DOD)  primarily  for  use  by  military  aircrews.  An 
analogous  set  of  publications  is  provided  by  the  Federal 
Aviation  Administration  (FAA)  for  use  by  civilian  pilots 
and  navigators. 
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Flow  Control:  A process  of  controlling  aircraft  in  order  to 


meet  a scheduled  sequence.  Flow  control  may  involve  delaying 
or  accelerating  particular  aircraft  in  a given  traffic 
situation . 

Initial  Approach  Fix  ( IAF ) ; A geographic  point  at  which  an 
Instrument  Approach  is  initiated.  The  IAF  is  normally  defined 
electronically  by  navigational  aids. 

Instrument  Approach:  For  the  purposes  of  this  paper  an 

instrument  approach  is  a non-radar  controlled  approach 
procedure  for  transitioning  an  aircraft  from  an  Initial 
Approach  Fix  to  a Final  Approach  Fix,  thence  to  a runway. 

In  flying  an  instrument  approach , a pilot  follows  a published 
procedure  and  proceeds  via  his  own  navigation.  Instrument 
approach  procedures  approved  by  the  DOD  for  military  air- 
craft are  published  in  FLIP. 

Instrument  Flight  Rules  (IFR):  A set  of  rules  provided  by 

the  FAA  for  the  conduct  of  flight.  These  rules  direct 
specific  behavior  by  both  pilot  and  air  traffic  controller. 

A pilot  operating  in  accordance  with  IFR  is  provides  separa- 
tion from  other  IFR  aircraft  by  an  air  traffic  control 
agency. 

Jet  Route:  One  of  a system  of  flight  paths  specified  for 

enroute  navigation  by  the  FAA.  Jet  routes  are  designed  for 
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operation  of  aircraft  between  the  altitudes  of  18,000  feet 
(mean  sea  level)  and  flight  level  450  (45,000  feet  mean  sea 
level) . Jet  routes  are  defined  by  specific  compass  bearings 
from  ground  based  transmitters  and  are  published  in  FLIP. 

Low  Altitude  Airway:  One  of  a system  of  flight  paths 

specified  for  enroute  navigation  by  the  FAA.  Low  altitude 
airways  are  designed  for  aircraft  operating  between  the 
altitudes  of  1200  feet  (above  ground  level)  up  to,  but  not 
including,  18,000  feet  (mean  sea  level).  The  airways  are 
defined  by  specific  compass  bearings  from  navigational  aids 
and  are  published  in  FLIP. 

Mean  Time  to  Fly  (MTF) ; The  average  time  required  for  an 
aircraft  of  a given  category  to  fly  between  two  points  with- 
out being  delayed  by  conflicting  traffic. 

Missed  Approach  (MA) : Failure  to  land  after  making  an 

approach.  A missed  approach  may  be  pre-planned,  as  in  the 
case  of  a training  maneuver,  or  spontaneous,  as  in  the  case 
of  improper  aircraft  alignment  with  the  runway.  Occupation 
of  the  runway  by  another  aircraft  is  another  condition 
requiring  a landing  aircraft  to  execute  a missed  approach. 

Runway  Occupancy  Time  (ROT) : The  amount  of  time  that  an  air- 

craft spends  on  the  runway  during  takeoff  or  landing  plus 
any  additional  time  on  the  runway  prior  to  takeoff  or  sub- 
sequent to  landing. 
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Scheduling : The  establishment  of  an  arrival  sequence  for 

aircraft  based  on  pre-flight  planning  data. 

Sequencing ; The  ordering  of  arriving  aircraft,  with  respect 
to  their  landing  position,  from  first  to  last. 

Standard  Vectoring  Routes  (SVR) : The  routes  normally  used 

by  controllers  to  maneuver  aircraft  from  Arrival  Fixes  to 
the  Final  Approach  Fix.  These  routes  may  be  prescribed  by 
published  directives  or  they  may  be  established  by  con- 
vention . 

Tower : A control  facility,  located  in  close  proximity  to  a 

runway,  which  provides  visual  monitoring  and  control  of 
arriving  and  departing  aircraft,  and  which  exercises  final 
clearance  authority  for  use  of  the  runway. 

Visual  Flight  Rules  (VFR) : A set  of  rules  provided  by  the 

FAA  for  the  conduct  of  flight.  Under  VFR  the  pilot  has 
more  flexibility  as  to  his  actions  than  under  IFR,  and  he 
retains  total  responsibility  for  the  avoidance  of  other 
aircraft. 

Assumptions  and  Limitations 

The  following  assumptions  and  limitations  apply  to 
the  development  of  the  Dayton  ACA  simulation  model: 

(1)  Consideration  of  the  final  approach  and  landing 
phases  of  flight  is  limited  to  Wr ight-Patterson 
AFB . 
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(2)  The  runway  in  use  during  any  simulation  run  is 
Runway  23,  the  primary  instrument  runway  at 
Wright-Patterson  AFB . 

(3)  The  effect  of  wind  on  the  behavior  of  an  air- 
craft is  not  considered. 

(4)  Altitude  is  considered  to  occur  in  discrete 
increments;  that  is  vertical  airspace  can  be 
thought  of  as  a series  of  surfaces  piled  one 
on  top  of  another.  An  aircraft  can  maneuver 

of  a "surface"  without  conflicting  with  an  air- 
craft on  another  "surface." 

(5)  Vectoring  routes  are  limited  to  a discrete 
number  of  alternatives. 

(6)  Airspeed  changes  are  considered  to  occur  in 
discrete  units. 

(7)  Separation  distances  between  aircraft  are  in 
accordance  with  FAA  and  DOD  criteria.  They 
vary  depending  on  aircraft  category. 

(8)  Human  error  on  the  part  of  pilots  or  air 
traffic  controllers  is  not  considered  in  the 
simulation  model. 

(9)  Unusual  situations  or  emergency  conditions  are 
not  considered  in  the  simulation  model. 
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Chapter  3 
METHODOLOGY 


The  methodology  in  this  research  involves  the  con- 
struction of  a simulation  model  of  the  Dayton  ACA  and  its 
subsequent  use  as  a tool  for  experimentation.  The  first 
step  in  model  construction  was  to  collect  data  to  satisfy 
information  requirements.  Information  was  required  in  two 
primary  areas:  information  on  the  physical  composition  of 

the  Dayton  ACA  and  information  on  the  nature  of  the  air 
activity  within  the  Dayton  ACA.  The  second  step  was  to  ana- 
lyze the  data  and  to  validate  the  results  in  preparation  for 
model  construction.  The  third  step  was  to  develop  the  simu- 
lation model,  itself,  in  the  form  of  a GPSS  computer  program 
(8).  The  development  process  was  based  directly  upon  the 
information  obtained  from  the  data  collection  and  analysis. 
The  fourth  step  was  to  validate  the  model  to  insure  that  the 
model  was  a true  representation  of  the  Dayton  ACA.  Valida- 
tion of  the  model  marked  the  completion  of  phase  one  as 
described  in  the  objective  section.  The  last  step  in  the 
methodology  was  to  apply  the  simulation  model  to  the  remain- 
ing phases  (two  and  three)  established  in  the  objective 
section. 

Phase  two  called  for  experimentation  with  the  simu- 
lation model.  This  was  accomplished  by  altering  the  original 
validated  model  in  a manner  that  would  test  the  significance 
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of  changes  in  ACA  design  and  in  ATC  procedures.  The  computer 
output  of  the  altered  models,  i.e.  the  experimental  models, 
was  compared  to  that  of  the  validated  model.  The  results  of 
these  experiments  are  recorded  in  Chapter  4.  Phase  three 
was  not  accomplished.  This  phase  is  discussed  in  Chapter  5 
under  Recommendations  for  Future  Research. 

This  methodology  chapter,  therefore,  is  organized  in 
accordance  with  the  five  steps  described  above  and  is  a 
systematic  elaboration  of  each  step,  beginning  with  data 
collection . 

Data  Collection 

As  noted  above,  information  was  required  in  two  pri- 
mary areas:  the  physical  composition  of  the  Dayton  ACA  and 

the  nature  of  the  air  activity  in  the  Dayton  ACA.  In  both 
areas  the  primary  emphasis  was  to  collect  data  relevant  to 
the  role  of  Wr ight-Patterson  AFB  within  the  Dayton  ACA  since 
the  simulation  model  was  limited  to  handling  air  traffic  to 
and  from  Wright-Patterson . Data  collection  on  the  physical 
composition  of  the  Dayton  ACA  is  presented  first. 

Composition  of  the  Dayton  ACA.  For  purposes  of  model  con- 
struction, the  data  for  the  Dayton  ACA  dealt  primarily  with 
the  physical  location  of  the  entities  that  make  up  the  Dayton 
ACA:  runways,  navigational  aids,  air  routes,  arrival  and 

departure  routes,  and  geographical  boundaries.  Given  the 
exact  location  of  these  entities,  precise  measurements  of 
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distance  relative  to  the  runway  at  Wright-Patterson  AFB  were 
made.  These  measurements,  in  nautical  miles,  were  used  in 
constructing  the  simulation  model. 

Data  pertaining  to  the  physical  characteristics  of 
the  Dayton  ACA  were  extracted  from  FLIP  and  from  other  docu- 
ments, maps,  and  charts  obtained  from  Dayton  Approach  Control 
(3).  FLIP  Enroute  High  Altitude  Chart  H-3  was  used  as  guid- 
ance for  the  jet  route  structure,  Figure  3,  that  overlies 
the  Dayton  ACA.  FLIP  Enroute  Low  Altitude  Charts  L-ll, 

L-21,  and  L-23  were  used  in  plotting  the  low  altitude  air- 
way structure  which  penetrates  the  controlled  airspace  of 
Dayton  Approach  Control.  Figure  4 is  a presentation  of  the 
low  altitude  airway  structure  in  relation  to  the  Dayton  ACA. 

Aircraft  inbound  to  the  Dayton  ACA,  on  an  IFR  flight 
plan,  proceed  via  the  jet  route  structure  or  via  the  low 
altitude  airway  network  and  are  transitioned  to  one  of  six 
arrival  fixes  (AFs)  that  service  Wright-Patterson  AFB.  In 
the  transitioning  process,  the  responsibility  of  air  traffic 
control  of  individual  aircraft  is  transferred  from  the 
enroute  air  traffic  control  facility  to  Dayton  Approach 
Control  who  then  controls  each  aircraft's  flight  path  to 
the  runway  at  Wright-Patterson.  The  location  of  the  six 
AFs  and  the  associated  air  routes  to  Wright-Patterson  were 
obtained  in  discussions  with  the  Deputy  Chief  of  the  Dayton 
Approach  Control  (3) . The  air  routes  from  the  AFs  to  the 
runway  are  the  Standard  Vectoring  Routes  (SVRs)  in  the 


Figure  3.  Jet  Routes. 


Figure  4.  Low  Altitude 


simulation  model  and  are  depicted,  along  with  the  six  AFs , 
in  Figure  5.  All  SVRs  converge  on  Runway  23  (see  Assumptions 
and  Limitations) . 

Air  traffic  characteristics.  In  order  for  the  simulation 
model  to  be  an  accurate  representation  of  the  Dayton  ACA, 
simulated  air  traffic  must  flow  through  the  model  as  real 
traffic  flows  through  the  Dayton  ACA.  Information  describing 
the  characteristics  of  the  air  traffic  within  Dayton  ACA  was 
needed.  The  source  of  this  information  was  the  flight 
progress  strip. 

A flight  progress  strip  (a  one  inch  by  eight  inch 
strip  of  paper)  contains  flight  information  on  an  individual 
aircraft.  The  amount  of  information  on  each  strip  varies 
and  is  a function  of  two  variables:  (1)  the  flight  rules 

under  which  the  aircraft  is  operating,  and  (2)  the  nature  of 
the  flight  activity.  Aircraft  are  operated  under  IFR  or 
under  VFR  and  air  activities  are  departure,  arrival  or 
through-flight  operations.  Normally,  flight  progress 
strips  contain  such  encoded  information  as  aircraft  call 
sign,  type  of  aircraft,  route  of  flight,  destination,  and 
proposed  arrival  or  departure  time.  The  regulations  govern- 
ing the  use  of  flight  progress  strips,  the  format  instruc- 
tions, and  information  codes  are  contained  in  FAA  Handbook 
7110.65.  Some  examples  of  flight  progress  strips  appear  in 
Figure  6. 
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Figure  6.  Sample  Flight  Progress  Strips 


For  flights  operating  in  accordance  with  IFR , flight 
progress  strips  are  computer-generated  and  supplied  to  an 
air  traffic  controller  before  the  aircraft  identified  on  the 
strip  enters  his  area  of  jurisdiction.  The  controller,  him- 
self, generates  flight  progress  strips  for  aircraft  operating 
in  accordance  with  VFR  when  he  initially  establishes  radio 
communications  with  the  pilot.  Air  traffic  controllers  use 
the  flight  progress  strips  for  planning  purposes  and  for 
recording  the  time  and  nature  of  particular  events  associated 
with  an  aircraft.  Examples  of  such  recordings  are  the  actual 
time  of  departure,  the  time  an  aircraft  penetrates  the  con- 
troller's airspace,  the  time  of  initial  radio  contact,  and 
the  time  of  transition  tb  another  control  facility. 

Since  flight  progress  strips  are  generated  for  every 
aircraft  operating  under  the  control  of  an  air  traffic  con- 
trol facility,  they  are  a source  of  raw  data  to  accurately 
indicate  the  characteristics  of  air  activity  in  a given 
area  (13).  Analysis  of  flight  progress  strips  can  provide 
information  on:  (1)  the  volume  of  air  traffic  for  any  given 

time  period,  (2)  the  frequency  and  the  types  of  air  opera- 
tions, (3)  the  types  of  aircraft  that  make  up  the  air 
traffic,  and  (4)  the  volume  of  traffic  using  navigational 
facilities  and  selected  airports  located  within  a given  area. 

Sampling  procedure.  All  the  flight  progress  strips  generated 
during  a single  day  at  a given  air  traffic  control  facility 
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are  packaged,  filed  and  retained  for  90  days.  In  order  to 
collect  data  to  identify  and  analyze  the  characteristics  of 
the  air  activity  in  the  Dayton  ACA,  packages  of  flight 
progress  strips  were  selected  from  the  files  of  Dayton 
Approach  Control.  Seven  packages  were  randomly  selected 
from  the  files  for  the  month  of  September,  1976  (one  package 
was  selected  from  among  all  of  the  "Sundays,"  one  from  among 
all  the  "Mondays,"  . . . , one  from  among  all  the  Saturdays"). 
In  this  respect  the  selected  flight  progress  strips  repre- 
sented an  entire  week  of  air  operations.  According  to 
Puhala,  September  was  typical  month  for  air  activity  in  the 
Dayton  ACA.  Nothing  unusual  occurred  in  September  that  would 
have  created  a greater  or  lesser  amount  of  air  traffic  than 
is  normally  experienced  (13) . 

The  seven  packages  contained  4526  flight  progress 
strips  (a  figure  obtained  by  summing  recorded  daily  totals) . 
From  this  population  a sample  of  311  flight  progress  strips 
was  selected.  The  selection  plan  for  the  sample  was  in 
accordance  with  a systematic  sampling  technique  as  outlined 
by  Emory  (5:153-154).  Every  fifteenth  flight  progress  strip 
was  selected  after  the  first  was  randonly  drawn. 

The  formula  used  to  calculate  the  sample  size  was 
taken  from  Emory  (5:150-152): 

±.06  = desired  maximum  error; 

1.09op  = .95  confidence  level  for  estimating  the  interval 
within  which  to  expect  the  population  proportion; 
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ap  = .03061  = standard  error  of  the  proportion  (.06/1.96) 
tt(I-tt)  = measure  of  population  dispersion,  here  estimated 
as  maximum  = .25; 

N = 4526  = population  size; 


op  = / TT  ( l-Ttf 

n-1 


N-n 

N-1 


n-,nr-,  / *25  / 4526-n 

.03061  = / r-  • / . c — r 

n-1  4526-1 


(. 03061) 2 = - 4526_n 


TT3T 


4525  (n-1)  (.03061)  = .25(4526-n) 


4 . 2398 (n) -4 . 2398  = 11 31 . 5- . 25 ( n) 


4.4898  (n)  = 1135.7398 


n = 252.96  = sample  size. 

Determination  of  the  minimum  sample  size  was  based  on 
an  estimated  dispersion  of  .25:  [.5(1— .5)],  since  there  was 

no  a priori  knowledge  of  tt,  and  on  a desired  confidence  inter*- 
val  of  12%.  Therefore,  the  minimum  sample  size  to  obtain  the 
desired  precision  is  approximately  253.  Since  there  existed 
considerable  uncertainty  as  to  the  frequency  of  occurrence  of 
the  various  attributes  of  population  aircraft,  it  was  decided 
that  the  sample  should  be  increased  above  the  minimum.  Thus, 
the  sample  size  was  set  at  300  to  insure  sufficient  accuracy 
would  be  obtained.  To  determine  the  selection  interval  for 
the  systematic  sample: 

4526  = 15.08 
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Selection  of  every  15th  flight  progress  strip  resulted  in  an 
actual  sample  of  311. 

Data  file  construction.  The  data  extracted  from  each  flight 
progress  strip  was  coded  in  the  following  manner  and  entered 
into  a data  file. 

Type  Aircraft:  "C-172,"  "F-100,"  "B-727,”  etc. 

Arrival,  Departure,  Through-flight:  "A,"  "D,"  "T,"  respec- 

tively . 

Time:  In  Greenwich  Mean  Time  (to  obtain  local  time  substract 

4 hours  since  September  occurred  during  Daylight  Savings 
Time) . 

Destination,  Departure  Airport:  "D"  = Dayton  Municipal 

Airport,  "W"  = Wr ight-Patterson  AFB,  "0"  = other. 

Last  Fix  (the  last  navigational  fix  an  arrival  aircraft 
crossed  prior  to  landing):  Examples  of  coding  are: 

"DAY"  = the  Dayton  VORTAC;  "FFO"  = the  Wr ight- 
Patterson  VORTAC;  etc. 

Previous  Fix  (the  next  to  last  navigational  fix  an  arrival 
aircraft  crossed  prior  to  landing) : Encoded  as  above. 

Departure  Fix  (the  first  navigational  fix  crossed  after 
takeoff  by  a departure  aircraft):  Encoded  as  above. 
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Next  Departure  Fix  (the  second  navigational  fix  crossed 


after  takeoff  by  a departure  aircraft) : Encoded  as  above. 

Flight  Rules:  "I"  = Instrument  Flight  Rules,  "V"  = Visual 

Flight  Rules. 

Day : "1"  = Sunday,  "2"  = Monday,  . . . , "7"  = Saturday. 

Category  (based  on  classification  by  airspeed:  Cruise, 

maneuvering  and  final  approach  airspeeds  in  knots 
groundspeed) : 

"V"  = 550,250,180 
"W"  = 500,220,170 
"X"  = 300,170,150 
"Y"  = 250,130,100 
"Z"  = 140,110,800 

The  data  file  contained  311  line  entries,  one  line 
for  each  flight  progress  strip.  Each  data  line  contained 
the  coded  information  described  above.  Figure  7 is  a portion 
of  the  resultant  data  file.  The  data  extracted  from  the 
flight  progress  strips  is  nominal  level  data  and  the  data 
file  was  formatted  to  permit  the  use  of  the  computer  pro- 
gram Statistical  Package  for  the  Social  Sciences  (SPSS)  in 
data  analysis. 

SPSS  is  an  integrated  system  of  computer  programs 
designed  to  provide  a unified  and  comprehensive  package 
that  enables  the  user  to  perform  different  types  of  data 
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analysis  in  a simple  and  convenient  manner  (12:1).  SPSS 
Version  6 was  used  to  evaluate  the  data  file  on  the  Honeywell 
635  computer  located  at  Headquarters,  Air  Force  Logistics 
Command,  Wr ight-Pat ter son  AFB , Ohio. 

Analysis  of  Flight  Progress  Strips 

Analysis  of  the  sample  data  (flight  progress  strips) 
was  accomplished  in  three  phases.  The  first  phase  was  to 
observe  the  frequency  of  occurrence  of  various  aircraft 
attributes.  SPSS  subprogram  FREQUENCIES  (12:194)  was  used 
initially  to  obtain  descriptive  presentations  of  the  data. 

The  data  was  then  analyzed  using  the  Chi-Square  test  for 
homogeneity  of  occurrence  ratios  (4:369-371,  380-382). 

The  second  phase  of  analysis  was  to  determine  the 
degree  to  which  one  characteristic  of  system  traffic  is 
dependent  upon  some  other  characteristic.  SPSS  subprogram 
CROSSTABS  (12:230)  was  used  to  produce  2-way  crosstabulations 
of  variables.  The  decree  of  dependency  of  the  variables 
based  on  the  distribution  of  frequency  counts  in  the  cross- 
tabulation tables  was  ascertained  by  the  use  of  the  Chi- 
Square  test  of  dependency  (4:376-378).  This  analysis  was 
particularly  useful  in  developing  the  simulation  model  since 
it  permitted  the  assignment  of  various  characteristics,  on 
a percentage  basis,  to  traffic  being  generated  in  the  simu- 
lation . 
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The  third  phase  of  analysis  was  to  identify  the 
distribution  of  air  activity.  Information  derived  through 
this  analysis  was  used  to  govern  the  generation  of  the 
various  types  of  air  activity  within  the  simulation  model. 

Data  analysis  follows  the  three  phases  described 
above.  Phase  one,  occurrence  ratios,  is  presented  first. 

Occurrence  ratios  of  attributes.  Each  aircraft  in  the  sample 
was  identified  by  seven  characteristics:  Operation,  time, 

airport,  route,  flight  rules,  day  and  aircraft  type.  The 
frequency  of  occurrence  of  the  various  characteristics  are 
shown  below  in  Table  1.  In  addition  to  the  characteristics 
shown  in  Table  1,  the  route  of  each  aircraft  was  identified, 
in  a general  sense,  by  the  last  navigation  fix  prior  to 
landing  or  the  first  navigation  fix  after  departure.  Based 
on  the  observed  characteristics,  several  hypotheses  were 
formulated  to  test  the  homogeneity  of  occurrence  ratios. 

The  first  hypothesis  concerns  types  of  operations. 

The  hypothesis  is  that  the  typed  of  operations  have  homo- 
geneous occurrence  ratios. 

H q : P (Operation^)  = 1/3 

: Not 

Operation:  A D 

Observed:  130  134 

Expected:  103.?"  103. 6 

<0i-Ei>2 

r /•  n n n 
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47 

103.  (T 
30.97 
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Table  1 


Characteristics  of  Traffic 


Through- 


Operations  : 
Arrival 
Departure 
Thro ugh- Flight 


Airport : 


Dayton 

153 

58.0 

60.8 

55.2 

- 

Wright- Pat ter son 

73 

27.7 

21.5 

33.6 

- 

Other 

38 

14.4 

17.7 

11.2 

- 

Flight  Rules: 

IFR 

226 

72.7 

76.2 

76.9 

51.1 

VFR 

85 

27.3 

23.8 

23.1 

48.9 

Day: 

Sunday 

50 

16.1 

14.6 

16.4 

19.1 

Monday 

38 

12.2 

13.1 

9.0 

19.1 

Tuesday 

54 

17.4 

16.2 

19.4 

14 . 9 

Wednesday 

51 

16.4 

12.3 

20.1 

17.0 

Thursday 

35 

11.3 

13.8 

8.2 

12.8 

Friday 

45 

14.5 

16.2 

14.2 

10.6 

Saturday 

38 

12.2 

13.8 

12.7 

6.4 

Aircraft  Type: 

V 

14 

4.5 

5.4 

5.2 

0 

W 

56 

18.0 

23.1 

17.9 

4.3 

X 

29 

9.3 

9.2 

12.7 

0 

Y 

34 

10.9 

11.5 

10.4 

10.6 

Z 

178 

57.2 

50.8 

53.7 

85.1 

Notes  on  Table  1: 

There 

were  47 

through- 

■flights  that 

did  not 

use  any  airport. 

Under 

Airport , 

"Other" 

means  any  of 

the 

10  small  civilian 

airfields  within  the  boundaries  of 

the 

Dayton  ACA. 


Total  Arrivals  Departures  Flights 
# % % % % 


130 

134 

47 


41.8 

43.1 

15.1 
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I 


= 46 . 54  >X2  = 9.21  a=.01  d.f.=c-l=2 

Therefore  is  rejected:  The  types  of  operations  do 

not  have  homogeneous  occurrence  ratios.  Lacking  any  other 
information,  the  best  estimate  of  the  actual  occurrence 
ratios  is  the  sample  observations. 

The  airports  at  which  arrivals  and  departures  occur 
is  an  important  consideration  in  developing  the  simulation 
model.  The  hypothesis  was  formulated  that  use  of  Dayton, 
Wright-Patterson,  and  "other"  airports  occur  in  homogeneous 
ratios : 


Hq : P(Airport^)  = 1/3 

H : Not 


Airport : 

D 

W 

0 

Observed 

38 

73 

153 

Expected 

88 

88 

88 

( 0 . -E . ) 2 

i l 

E. 

l 

28.41 

2.56 

48.01 

v2  = 78.98>X2 

= 9.21 

a= . 01 

d . f . =c-l=2 

''S  c 

Therefore 

Ho  is 

rejected 

: The  airports  are 

not  used 

in  homogeneous  ratios. 

The  best 

estimate  available 

of  the 

actual  occurrence 

ratios 

is  the 

sample  data. 

The  flight  rules  under  which  an  aircraft  is  operating 
are  important  because  they  may  effect  the  manner  in  which  an 
aircraft  is  handled.  The  hypothesis  is  that  IFR  and  VFR 
occur  with  equal  probability. 
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Hq  : P(IFR)  = 

P(VFR)  = 

1/2 

H1:  Not 

Rules : 

IFR 

VFR 

Observed : 

226 

85 

Expected : 

155.5 

155 

( | 0i-Ei|-l/2) 

2 

Ol  Cl 

: 31.51  31.51 


J 


6.63  ot= . - 01  d.f.=l 


Therefore  HQ  is  rejected:  The  probability  of  occur- 

rence is  not  equal. 

Finally,  the  occurrence  ratios  of  aircraft  types 
were  hypothesized  to  be  equal. 


Hq : P (Aircraft  Type 

H1:  Not 

i)  = 1/5 

Type : 1 

2 

3 

4 

5 

Observed:  14 

56 

29 

34 

178 

Expected:  62.2 

62.2 

62.2 

62.2 

62.2 

(0  -E  )2 

— i : 37.34 

Ei 

.618 

17.72 

12.78 

215.58 

X2  = 284 . 05>X2  = 13. 
As  c 

27  o=.01 

d.f .=4 

Therefore  is 

rejected : 

Aircraft 

types 

do  not  have 

homogeneous  occurrence  ratios. 


♦This  formula  was  required  because  of  the  number  of  cells. 
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The  results  of  these  tests  indicate  that  in  all  cases 


the  frequency  of  occurrence  was  not  homogeneous.  The  fre- 
quencies observed  in  the  sample  are  assumed  to  provide  the 
best  estimate  of  the  actual  occurrence  ratios. 

Dependency  between  attributes.  The  second  phase  of  data 
analysis  involved  tests  for  dependency  conducted  on  the  fol- 
lowing air  traffic  characteristics: 

1.  Operation  - Aircraft  Type 

2.  Airport  - Aircraft  Type 

3.  Flight  Rules  - Aircraft  Type 

4.  Flight  Rules  - Airport 

5.  Airport  - Arrival  Route 

6 . Airport  - Departure  Route 

7.  Aircraft  Type  - Day  of  the  Week 

8.  Aircraft  Type  - Time  of  Day 

9.  Airport  - Operation 

The  first  hypothesis  tested  concerned  operation  and 
aircraft  type.  The  hypothesis  is  that  operation  and  aircraft 
type  occur  independently. 

Hq:  P (Operational  Aircraft  Type^)  = P (Operation^) 

Hx:  Not 

X2  = 22.99>X^  = 20.09  o».01  d.f.=8 

As  c 

Therefore  is  rejected:  Dependency  exists  between 

2 

type  of  aircraft  and  type  of  operation.  The  x contingency 
table  used  to  revolve  this  hypothesis  is  shown  in  Table  2. 
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Table  2 


Crosstabulation:  Aircraft  Category 

and  Operation 

ARRDER 

COUNT  I 

ROW  PCT  IARRIVAL  DEPARTR  THRU  FLT  ROW 


COL  PCT 
TOT  PCT 

I 

I 

1 

.1 

2 

r 

• X 
_ r _ 

3 

.1 
_ T 

TOTAL 

LAI  “ ~ “ 

0. 

I 

7 

I 

7 

i 

0 

I 

14 

FAST  MOVER 

I 

50.0 

I 

50.0 

i 

0. 

I 

4.5 

I 

5.4 

I 

5.2 

i 

0. 

I 

I 

2.3 

I 

2.3 

i 

0. 

I 

-I- 

-I- 

-I- 

-I 

0. 

I 

30 

I 

24 

i 

2 

I 

55 

AIR  CARRIER 

I 

53.6 

I 

42.9 

i 

3.6 

I 

18.0 

I 

23.1 

I 

17.9 

i 

4.3 

I 

I 

9.6 

I 

7.7 

i 

0.6 

I 

-I- 

-I- 

-i- 

-I 

0. 

I 

12 

I 

17 

i 

0 

I 

29 

MEDIUM 

I 

41.4 

I 

58.6 

i 

0. 

I 

9.3 

I 

9.2 

I 

12.7 

i 

0. 

I 

I 

3.9 

I 

5.5 

i 

0. 

I 

-I- 

-I- 

-i- 

-I 

0. 

I 

15 

I 

14 

i 

5 

I 

34 

SMALL 

I 

44.1 

I 

41.2 

i 

14.7 

I 

10.9 

I 

11.5 

I 

10.4 

i 

10.6 

I 

I 

4.8 

I 

4.5 

i 

1.6 

I 

-I- 

-I- 

-i- 

-I 

0. 

I 

66 

I 

72 

i 

40 

I 

178 

SLOW 

I 

37.1 

I 

40.4 

i 

22.5 

57.2 

I 

50.8 

I 

53.7 

i 

85.1 

I 

I 

21.2 

I 

23.2 

i 

12.9 

I 

-I- 

-I- 

-i- 

-I 

COLUMN 

130 

134 

47 

311 

TOTAL 

41.8 

43.1 

15.1 

100.0 

CHI  SQUARE  = 22.98888  WITH  8 DEGREES  OF  FREEDOM 
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Given  an  aircraft  type,  the  bivariate  probabilities  shown  in 
the  contingency  table,  therefore,  are  the  best  information 
available  for  designating  an  aircraft  generated  in  the  model 
as  an  arrival,  a departure  or  a through-flight. 

The  second  hypothesis  is  that  airport  usage  occurs 
independently  of  aircraft  type. 

Hq : P (Airport^ | Aircraft  Type^)  = P (Airport^) 

Hx : Not 

Xg  = 139. 11>X2  = 20.10  a= . 01  d.f.=8 

Therefore  Hq  is  rejected:  Dependency  exists  between 

2 

airport  usage  and  aircrart  type.  Table  3 is  the  x con- 
tingency table  used  to  revolve  the  hypothesis.  The  infor- 
mation in  the  contingency  table  will  be  used  to  assign 
aircraft,  generated  in  the  model,  to  an  airport. 

The  third  hypothesis  is  that  flight  rules  are  inde- 
pendent of  aircraft  type. 

Hq:  P ( IFR  j Aircraft  Type^)  = P ( IFR  j Airf  raf  t Type_.) 

= P ( IFR) 

H,:  Not 

X2  = 63.4  4 >X2  = 13.277  a=.01  d.f.=4 

A s c 

Therefore  Hq  is  rejected:  Dependency  exists  between 

flight  rules  and  aircraft  type.  Given  an  aircraft  type,  the 
best  estimate  available  on  the  percentage  of  the  type  of 
aircraft  that  will  be  IFR  is  the  percentage  obtained  from 
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Table  3 


Crosstabulation : Aircraft  Category 

and  Airport 


CAT 

COUNT 
ROW  PCT 
COL  PCT 
TOT  PCT 

I 

I OTHER 

I 

I 

. T 

WPAFB 

I 

I 

__  T 

DAYTON 

I 

_ T 

ROW 

TOTAL 

0. 

I 

0 

I 

11 

I 

3 

I 

14 

FAST  MVR 

I 

0. 

I 

78.6 

21.4 

I 

5.3 

I 

0. 

I 

15.1 

I 

2.0 

I 

I 

0. 

I 

4.2 

I 

1.1 

I 

- 

•I  — 

— 

-I  — 

-I 

-I 

0. 

I 

0 

I 

4 

I 

50 

I 

54 

CARRIER 

I 

0. 

I 

7.4 

I 

92.6 

I 

20.5 

I 

0. 

I 

5.5 

I 

32.7 

I 

I 

0. 

I 

1.5 

I 

18.9 

I 

- 

■I— 

-I  — 

— 

-I 

-I 

0. 

I 

1 

I 

28 

I 

0 

I 

29 

MEDIUM 

I 

3.4 

I 

96.6 

I 

0. 

I 

11.0 

I 

2.6 

I 

38.4 

I 

0. 

I 

I 

0.4 

I 

10.6 

I 

0. 

I 

- 

I — 

-I— 

-I 

-I 

0. 

I 

9 

I 

11 

I 

9 

I 

29 

SMALL 

I 

31.0 

I 

37.9 

I 

31.0 

I 

11.0 

. 

I 

23.7 

I 

15.1 

I 

5.9 

I 

I 

3.4 

I 

4.2 

I 

3.4 

I 

- 

I — 

— 

-I  — 

— 

-I 

-I 

0. 

I 

28 

I 

19 

I 

91 

I 

138 

SLOW 

I 

20.3 

I 

13.8 

I 

65.9 

I 

52.3 

I 

73.7 

I 

26.0 

I 

59.5 

I 

I 

10.6 

I 

7.2 

I 

34.5 

I 

- 

•I  — 

-I  — 

— 

-I' 

-I 

COLUMN 

38 

73 

153 

264 

TOTAL 

14.4 

27.7 

58.0 

100.0 

CHI  SQUARE 

= 139.11394 

WITH 

8 

DEGREES 

OF  FREEDOM 
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the  sample.  The  contingency  table  used  to  resolve  this 
hypothesis  is  shown  in  Table  4. 

The  third  hypothesis  dealing  with  flight  rules  is 
that  flight  rules  are  independent  of  airport. 

Hq  : P (IFR|  Airport^  = P ( IFR  | Airport  ^ ) = P(IFR) 

: Not 

X2  = 9.25>X2  = 9.21  a-. 01  d.f.=2 

s c 

Therefore  is  rejected:  Dependency  exists  between 

airport  used  and  flight  rules.  Since  the  conditional  prob- 
abilities are  not  equal,  the  best  estimate  of  the  actual 
probabilities  are  those  obtained  from  the  sample  data.  This 
information  will  be  used  to  assign  flight  rules  to  an  air- 
craft based  on  its  airport  of  origin  or  destination.  The 
2 

X contingency  table  used  to  resolve  this  hypothesis  is 
shown  in  Table  5. 

Airport  was  next  crosstabulated  with  route  for  both 
arrivals  and  departures.  It  should  be  noted  that  certain 
combinations  of  airport  and  route  are  precluded  by  airspace 
divisions  within  the  Dayton  ACA.  For  example,  an  arrival  via 
70U  cannot  go  to  Dayton.  Similarly,  a departure  from  Wright- 
Patterson  cannot  go  to  9WJ.  These  limitations  result  in 
contingency  tables  that  do  not  meet  the  accepted  minimum 
cell  size  criteria  and  are  technically  invalid.  Since  most 
of  the  zero  value  cells  are  known  to  have  zero  value  in  the 
population  of  actual  air  traffic,  this  technical  limitation 
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Table  4 


Crosstabulation:  Aircraft  Category 

and  Flight  Rules 


COUNT 

I 

ROW  PCT 

IIFR 

VFR 

COL  PCT 

I 

ROW 

TOT  PCT 

I 

20 

. I 

21 

.1 

TOTAL 

CAT  

•I- 

-I  — 

— 

-I 

0. 

I 

14 

I 

0 

I 

14 

FAST  MOVER 

I 

100.0 

I 

0. 

I 

4.5 

I 

6.2 

I 

0. 

I 

I 

4.5 

I 

0. 

I 

- 

■I- 

-I  — 

-I 

0. 

I 

56 

I 

0 

I 

56 

AIR  CARRIER 

I 

100.0 

I 

0. 

I 

18.0 

I 

24 . 8 

I 

0. 

I 

I 

18.0 

I 

0. 

I 

- 

•I- 

-I  — 

-I 

0. 

I 

28 

I 

I 

29 

MEDIUM 

I 

96.6 

I 

3.4 

I 

9.3 

I 

12.4 

I 

1.2 

I 

I 

9.0 

I 

0.3 

I 

- 

I- 

-I  — 

— 

-I 

0. 

I 

29 

I 

5 

I 

34 

SMALL 

I 

85.3 

I 

14.7 

I 

10.9 

I 

12.8 

I 

5.9 

I 

I 

9.3 

I 

1.6 

I 

- 

•I- 

— 

-I  — 

— 

-I 

0. 

I 

99 

I 

79 

I 

178 

SLOW 

I 

55.6 

I 

44.4 

I 

57.2 

I 

43.8 

I 

92.9 

I 

I 

31.8 

I 

25.4 

I 

- 

I- 

-I— 

— 

-I 

COLUMN 

226 

85 

311 

TOTAL 

72.7 

27.3 

100.0 

CHI  SQUARE  = 63.44043 

WITH 

4 

DEGREES 

OF  FREEDOM 
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Table  5 


Crosstabulation:  Airport 

and  Flight  Rules 


COUNT  I 


ROW 

PCT 

I IFR 

VFR 

ROW 

COL 

PCT 

I 

TOTAL 

TOT 

PCT 

I 

20 

.1 

21 

.1 

DESDEP 

— 

■I- 

-I  — 

-I 

0. 

I 

34 

I 

4 

I 

38 

OTHER 

I 

89.5 

I 

10.5 

I 

14.4 

I 

16.8 

I 

6.5 

I 

I 

12.9 

I 

1.5 

I 

- 

•I- 

-I-- 

-I 

0. 

I 

61 

I 

12 

I 

78 

WRIGHT 

PAT 

I 

83.6 

I 

16.4 

I 

27.7 

I 

30.2 

I 

19.4 

I 

I 

23.1 

I 

4.5 

I 

- 

•I- 

-I-- 

-I 

0. 

I 

107 

I 

46 

I 

153 

DAYTON 

I 

69.9 

I 

30.1 

I 

58.0 

I 

53.0 

I 

74.2 

I 

I 

40.5 

I 

17.4 

I 

- 

I- 

-I-- 

-I 

COLUMN 

202 

62 

264 

TOTAL 

76.5 

23.5 

100.0 

CHI  SQUARE  = 9.25  WITH  2 DEGREES  OF  FREEDOM 


is  considered  irrelevant.  The  contingency  tables  shown  in 
Tables  6 and  7 are  therefore  considered  to  be  simple  prob- 
ability tables  and  the  resulting  percentages  are  treated  as 
representative  of  population  values. 

The  type  of  aircraft  is  an  important  characteristic 
of  a traffic  element  as  it  will  determine  speed,  following 
distances,  and  times  to  fly  between  geographic  point.  It 
was  valuable  to  determine  if  aircraft  type  is  dependent  on 
other  characteristics.  Thus  far,  the  dependency  of  aircraft 
type  has  been  evaluated  with  respect  to  type  of  operation, 
airport,  and  flight  rules.  Two  additional  hypotheses  were 
tested  regarding  aircraft  type. 

The  first  hypothesis  is  that  aircraft  types  occur 
independent  of  the  day  of  the  week. 

Hq:  P (Aircraft  Type^jDay^)  = P (Aircraft  Type^Day^) 

= P (Aircraft  Type^) 

H^:  Not 

X2  = 18 . 88<X2  = 34.81  a=.01  d.f.=18 

As  c 

Therefore  cannot  be  rejected:  There  is  insuffi- 

cient evidence  to  assess  dependency  between  aircraft  type 
and  day  of  the  week.  Consequently,  aircraft  in  the  simu- 
lation will  be  assigned  an  aircraft  type  without  regard  to 

2 

the  day  of  the  week.  The  x contingency  table  used  to 
resolve  this  hypothesis  is  shown  in  Table  8. 


49 


H 

M 

i 

H 

M 

M 

H 

M 

i 

M 

H 

H 

M 

H 

H 

M 

H 

IH 

H 

i 

O 

1 

1 

1 

1 

1 

1 

1 

1 

H 

1 

IN 

*H 

cn 

O 

1 

cn 

rH 

r- 

rH 

1 

o 

1 

in 

CN 

1 

• 

• 

• 

1 

rH 

• 

• 

• 

1 

• 

• 

• 

1 

H 

• 

1 

rH 

m 

CN 

1 

CO 

cn 

1 

O 

o 

o 

1 

m 

2 

1 

H 

1 

H* 

00 

H 

1 

1 

H 

O 

H 

1 

1 

H 

i 

H 

H 

M 

M 

1 

1 

M 

M 

H 

M 

H 

1 

1 

H 

1 

H 

M 

M 

H 

1 

1 

H 

<J\ 

1 

1 

1 

rH 

KD 

r- 

O 

1 

1 

1 

o 

1 

1 

1 

CN 

CN 

cn 

rH 

1 

1 

1 

m 

rH 

cn 

1 

• 

• 

• 

1 

• 

• 

• 

1 

rH 

• 

• 

• 

1 

rH 

• 

oj 

1 

m 

r- 

*H 

1 

o 

o 

o 

! 

CN 

CN 

CN 

1 

m 

■u 

23 

1 

» 

1 

CN 

CT> 

rH 

1 

rH 

a 

D 

1 

1 

1 

1 

0 

r~ 

1 

1 

1 

1 

os 

H 

H 

H 

H 

M 

M 

H 

H 

H 

M 

M 

M 

M 

M 

H 

H 

M 

rH 

in 

1 

i 

1 

1 

1 

1 

1 

1 

cfl 

i 

CN 

rH 

o 

1 

CO 

00 

O 

1 

O 

00 

rH 

1 

m 

rH 

> 

i 

• 

• 

• 

1 

• 

• 

• 

1 

• 

• 

• 

1 

rH 

• 

•H 

i 

i H 

in 

CN 

1 

o 

*1* 

1 

cn 

on 

1 

m 

E-| 

i 

rH 

iH 

1 

iH 

m 

1 

rH 

in 

1 

rH 

n 

< 

i 

1 

1 

1 

< 

M 

H 

i 

M 

M 

H 

H 

M 

1 

H 

M 

M 

M 

H 

1 

M 

M 

w 

M 

H 

1 

H 

T> 

• 

1 

1 

1 

1 

C 

1 

1 

1 

1 

(0 

1 

O 

1 

•H 

r*- 

o 

o 

1 

O 

1 

rH 

o 

1 

• 

• 

• 

1 

• 

• 

• 

1 

• 

• 

• 

1 

• 

• 

4J 

1 

O 

o 

O 

1 

cn 

o 

rH 

1 

o 

o 

O 

1 

H 

ID 

23 

1 

1 

o 

1 

1 

0 

03 

1 

1 

rH 

1 

1 

a; 

a. 

m 

1 

1 

1 

1 

rH 

n 

M 

M 

M 

H 

H 

M 

W 

M 

H 

H 

IH 

H 

IH 

H 

H 

H 

H 

J3 

•H 

• 

1 

1 

1 

1 

(0 

< 

cn 

1 

1 

1 

1 

e* 

1 

O 

1 

in 

in 

o 

rH 

1 

O 

1 

in 

rH 

1 

• 

• 

• 

1 

• 

• 

• 

1 

• 

• 

• 

1 

• 

•• 

1 

O 

O 

O 

1 

00 

o 

in 

1 

o 

O 

o 

1 

m 

c 

33 

1 

1 

rH 

o 

1 

1 

0 

O 

1 

1 

rH 

1 

1 

•H 

1 

1 

1 

1 

4J 

H 

M 

M 

M 

M 

M 

H 

H 

H 

H 

H 

M 

H 

H 

H 

H 

H 

03 

• 

1 

1 

1 

1 

rH 

CM 

1 

1 

1 

1 

3 

1 

cn 

in 

O 

1 

rH 

r- 

CN 

o 

1 

O 

o 

cn 

CN 

1 

CN 

XI 

1 

• 

• 

• 

1 

• 

• 

• 

1 

CN 

• 

# 

• 

1 

CN 

• 

<3 

1 

CN 

1 

H 

1 

r- 

m 

O 

1 

'T 

-u 

a 

1 

*H 

*H 

1 

1 

cn 

00 

CN 

1 

CN 

in 

M 

1 

1 

1 

1 

en 

OS 

1 

1 

1 

1 

0 

H 

M 

M 

H 

M 

H 

M 

M 

HH 

M 

H 

H 

H 

M 

M 

H 

« Q 

e-  o 

os  os 


E-i  E-*  E-* 
E-i  U U U 
z a*  a < eu 
D 

O 2 J E- 
UOOO 
« u e- 


o r'.  4— i 

in  m o 
in  m h 


Mrtr«o 
• • • 

rH  O OO 


in  co  vo  cn 


r-*  m m 
(N  in  h 


50 


CHI  SQUARE  = 67.94834  WITH  12  DEGREES  OF  FREEDOM 


Table 


M 

►H 

*H 

b-i 

t— t t— • 

b-i 

b-i 

M 

M M 
■ 

*-« 

H 

M 

M b-i 

CO 

1 

I 

1 

rH 

O 

1 

O 1 

r—i 

a> 

O 

1 

1 

O 1 

CO 

O 

1 

1 

X 1 

X 

1 

• 

• 

• 

• 

o 1 

• 

• 

r—i 

• 

1 

X 

o 

H | 

CN 

o 

1— 1 1 1 

>n 

o 

r*  i 

<7* 

O 

1 

*H 

1 

l—i 

1 

H 

CO 

l 

D 

1 

1 

1 

i 

r- 

hi 

1 

b-i 

b-i 

*— « 

1 

b-i  HI 

b-i 

b-i 

b-i 

1 

b-i  b-i 

b-i 

b-i 

i 

b-i  b-i 

r* 

1 

1 

1 

<H 

r* 

o 

1 

1 

O 1 

O 

1 

1 

1 

—i 

O' 

o 

1 

1 

O 1 

CN 

O' 

1 

• 

• 

• 1 

• 

• 

• 

• 

• 

cn 

1 

X 

o 

i-H  | 

o 

o 

O 1 

c 

— . 1 

D 

»"3 

1 

LD 

1 

1 

n 

1 

4J 

2 

1 

1 

1 

1 

3 

1 

1 

1 

1 

0 

M 

w 

b-i 

M 

M 

M M 

b-i 

M 

w 

W M 

b-i 

M 

M 

*-<  b-i 

a: 

• 

1 

1 

1 

1 

X 

1 

1 

I 

1 

<D 

1 

CN 

ro 

ro 

<T>  1 

X 

rH 

O 

CO  1 

CN 

r* 

n i 

"T 

ro 

M 

1 

• 

• 

0 

• 

• 

• 

CN 

• 

3 

1 

ro 

SO 

r—i  | 

r- 

n 

in  i 

o 

X 

n i 

ro 

4J 

X 

1 

<H 

1 

rH 

CN 

i 

ro 

X 

rH  1 

rj 

z 

1 

1 

i 

1 

(Q 

r- 

1 

1 

i 

\ 

a 

M 

M 

►H 

IH 

►H 

t-H  b-i 

I 

b-i 

M 

M 

W M 

b-i 

Hi 

M 

M b-i 

S 

U1 

1 

1 

1 

ro 

o 

ro 

1 

1 

<y*  i 

X 

rH 

X 

1 

1 

CO  I 

CN 

lO 

1 

1 

r*»  i 

T3 

\ 

• 

• 

• i 

0 

• 

• 1 

f ) 

• 

• 

• | 

CN 

• 

C 

1 

o 

*T 

CN  | 

r" 

CO 

in  i 

CN 

r- 

•H  | 

O 

fTJ 

1 

CN 

*H 

1 

rH 

CN 

i 

CN 

n 

r—i  j 

CN 

< 

1 

1 

i 

1 

4J 

W 

1 

1 

i 

1 

^4 

<rs 

M 

b-i 

1 

M 

M 

M 

M b-i 

b-i 

M 

b-i 

W M 
■ 

M 

b-i 

M 

M M 

8. 

rr 

1 

1 

1 

1 

1 

1 

1 

1 

1 

O 

1 

<T> 

O 

O 1 

O 

1 

r—i 

O 

•H 

1 

• 

• 

• 

• 

• 

• 

• 

< 

1 

o 

o 

C 1 

cn 

o 

rH  | 

o 

c 

O 1 

H 

IS 

1 

1 

c 

1 

1 

N 

1 

1 

r—i 

1 

1 

•• 

ro 

1 

1 

1 

1 

c 

M 

M 

M 

►H 

M 

W M 

M 

b-i 

b-i 

b-i  M 

M 

b-i 

IH 

W-l  b-i 

0 

• 

1 

1 

1 

1 

•H 

ro 

1 

1 

1 

1 

4J 

1 

O 

1 

o 

0>  1 

O 

1 

3 

1 

• 

• 

• 1 

» 

• 

• 1 

• 

• 

• 

r—i 

1 

o 

o 

O 1 

f—i 

o 

ro  | 

o 

o 

O 1 

ro 

3 

D 

1 

1 

r—i 

a 

1 

1 

A 

O 

1 

1 

r—i 

1 

1 

<q 

r* 

1 

1 

1 

1 

•*j 

M 

M 

M 

M 

M 

b-i  M 

b-i 

b-i 

M 

W b-i 

M 

w 

b-i 

b-i  W-« 

cn 

• 

1 

1 

1 

1 

<n 

CN 

1 

1 

1 

1 

0 

1 

in 

ro 

CO 

i 

o 

X 

sQ 

i 

& 

ro 

X 

CC  1 

r—i 

'-r 

H 

1 

• 

• 

— 4 

• 

• 

• | 

• 

• 

CN 

• 

U 

1 

ro 

ro 

tT  | 

CO 

r* 

l 

i — 1 

X 

n j 

C 

a 

H 

1 

1 

ro 

CN 

i 

i 

CN 

1 

1 

H 

CN 

I 

i 

CN 

a 

M 

1 

HH 

HH 

M 

M 

i 

b-i  b-i 

M 

M 

M 

1 

b-i  b-i 

b-i 

b-i 

b-i 

i 

*H  b-i 

<H 

1 

1 

1 

ro 

O 

O 

1 

1 

O'  1 

r* 

o 

o 

i 

1 

CO  1 

o 

o 

1 

r*  i 

O 

TT 

1 

• 

• 

• 1 

• 

• 

rH 

• 

• 

• i 

CN 

• 

1 

O 

(N  I 

o 

in 

lO  1 

00 

o 

ON  1 

RTE 

ROD 

1 

1 

1 

CN 

«H 

1 

1 

1 

CN 

ro 

1 

1 

1 

n 

1 

1 

1 

rH 

M »H 

M HH 

b-i 

M 

M 

w 

M b-i 

b-i 

b-i 

M 

M b-i 

M 

M 

b-i 

b-i  b-i 

Z 0-0-3- 

p 5 j e 

O o o c 

Qi  O r 


a: 

CQ 

z 

o 

M 

Ci- 

p 

X 

te 

> 

6 

o 

Si 

< 

Q 
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CHI  SQUARE  = 18.88173  WITH  18  DEGREES  OF  FREEDOM 


The  second  hypothesis  is 

independently  of  the  time  of  day 

H„:  P (Aircraft  Type.! Time.) 

0 11 

= P (Aircraft  Type J 


that  aircraft  type  occurs 
(one  hour  periods) . 

= P (Aircraft  Type^  , Time . ) 


H.:  Not 

X?  = 6.25<xl  = 20.09  a=.01  d.f.=8 

s 

Therefore  HQ  cannot  be  rejected:  There  is  insuffi- 

cient evidence  to  assess  dependency  between  aircraft  type 
and  time  of  day.  The  contingency  table  used  to  resolve  this 
hypothesis  is  shown  in  Table  9.  The  table  required  that 
"time"  be  collapsed  into  three  periods.  The  time  period 
labeled  "Daytime"  covers  the  nine  busiest  hours  of  the  day. 

The  final  hypothesis  tested  deals  with  the  relation- 
ship between  operation  and  airport.  The  hypothesis  is  that 
operation  is  independent  of  airport. 

Hq:  P (Operation^ | Airport^ ) = P (Operation^ , Airport^ ) 

= P (Operation^) 

H1 : Not 

X2  = 1 . 73  <X2  = 9.21  a*. 01  d.f.=2 

s c 

Therefore  HQ  cannot  be  rejected:  There  is  insuffi- 

cient evidence  to  assess  dependency  between  operation  and  air- 
port. The  contingency  table  used  to  resolve  this  hypothesis 
is  shown  in  Table  10. 

Dependency  was  found  in  all  but  the  last  three  tests. 
The  fact  that  aircraft  type  occurs  independent  of  either  day 
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CH_  SQUARE  = G. 24562  WITH  8 DEGREES  OF  FREEDOM 


Table  10. 


Crosstabulation:  Operation 

and  Airport 


DESDEP 

COUNT 

I 

ROW 

PCT 

IOTHER 

WPAFB 

DAYTON 

ROW 

COL 

PCT 

I 

TOTAL 

TOT 

PCT 

I 

I 

I 

I 

ARRDEP 

— 

— 

■I 

-I 

-I 

-I 

1. 

I 

19 

I 

26 

I 

54 

I 

99 

ARRIVAL 

I 

19.2 

I 

26.3 

I 

54.5 

I 

49.0 

I 

55.9 

I 

42.6 

I 

50.5 

I 

I 

9.4 

I 

12.9 

I 

26.7 

I 

- 

I' 

-I 

-I 

-I 

2. 

I 

15 

I 

35 

I 

53 

I 

103 

DEPARTR 

I 

14.6 

I 

34.0 

I 

51.5 

I 

51.0 

I 

44 . 1 

I 

57.4 

I 

49.5 

I 

I 

7.4 

I 

17.3 

I 

26.2 

I 

- 

-I 

-I 

-I 

COLUMN 

34 

61 

107 

202 

TOTAL 

16.8 

30.2 

53.0 

10 0.  0 

CHI  SQUARE  = 1.72927  WITH  2 DEGREES  OF  FREEDOM 


of  the  week  or  time  of  day  indicates  the  sample  drawn  from  an 
entire  week  of  air  activity  reflects  the  occurrence  of  air- 
craft types  during  any  time  period  within  that  week. 

The  dependency  relationships  established  by  the  above 
tests  were  used  to  assign  characteristics  to  air  traffic 
generated  in  the  simulated  system.  The  probability  that  an 
aircraft  will  possess  a specific  combination  of  character- 
istics can  be  viewed  as  the  extremity  of  a branching  network; 
Figure  8 shows  how  the  occurrence  ratios  and  the  dependency 
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Route 


Figure  8.  Branching  Network 


relationships  were  used  to  assign  characteristics.  The  num- 
ber of  aircraft  possessing  the  combination  of  characteristics 
contained  in  each  branch  is  the  product  of  the  probabilities 
contained  in  that  branch  and  the  number  that  entered  the 
network. 

Distribution  of  air  activity.  An  important  element  in  the 
simulation  model  is  the  rate  at  which  aircraft  arrive  at 
and  depart  from  Wright-Patterson  AFB.  As  a means  of  identify- 
ing the  distribution  of  air  traffic,  two  procedures  were 
employed.  First,  the  distribution  of  IFR  departures  from 
Wright-Patterson  was  determined  using  curve-fit  procedures 
(4:371-376).  Second,  the  total  number  of  air  operations 
occurring  during  certain  peak  hours  was  used  to  determine 
the  number  of  IFR  departures  from  Wright-Patterson  based  on 
the  probability  of  an  air  operation  having  that  combination 
of  characteristics.  The  results  of  these  two  procedures  were 
then  compared. 

Analysis  of  the  sample  of  311  aircraft  previously 

described  indicated  that  the  period  from  1000  to  1900  had 

the  greatest  amount  of  aircraft  activity.  In  addition,  the 

2 

traffic  during  that  period  was  relatively  constant.  A x 
test  for  homogeneity  of  occurrence  ratios  was  conducted  as 
follows : 


V 

P (Hour ^ ) 

= 1/9 

H1 : 

at  least 

one  not  equal 
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, ! (O.-E.  ) 

2 _ r ( l 1 I 


Hour 

Obs 

Exp 

(°i-Ei) 

E . 
1 

1 

23 

22.11 

. 0358 

2 

18 

22.11 

.7640 

3 

26 

22.11 

. 6844 

4 

22 

22.11 

.0005 

5 

18 

22.11 

. 7640 

6 

18 

22.11 

. 7640 

7 

23 

22.11 

. 0358 

8 

25 

22.11 

.3777 

9 

26 

22.11 

.6844 

4.1106 

2 

*c 

20.090 

a=.01  d. 

f .=c-l=8 

2 

*s 

= 20.09 

2 

Therefore  xs  falls  within  the  acceptance  region  and 
Hq  cannot  be  rejected.  This  result  indicates  that  during  the 
nine-hour  period  sampled,  the  hours  are  relatively  homo- 
geneous with  respect  to  the  rate  of  aircraft  activity. 

Based  on  the  result  of  the  foregoing  analysis,  the 
flight  progress  strips  were  obtained  for  all  departures  from 
Wright-Patterson  during  the  priod  from  1000  to  1900  on  the 
busiest  day  in  September,  1976.  There  were  51  such  departures 
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recorded.  The  hypothesis  was  then  formulated  that  these 
observed  departures  were  obtained  from  a population  following 
a Poisson  distribution. 

Hq : x~Poisson  (X) 

: Not 

The  random  variable  x is  defined  as  the  number  of 
departures  in  a ten  minute  period.  There  were  54  such  10 
minute  periods  observed.  The  parameter  X of  the  Poisson 
distribution  was  estimated  from  the  sample  mean  arrivals 
per  period  (x) . The  hypothesis  test  was  resolved  as  shown 
in  Table  11. 


Table  11. 

Chi-Square  Goodness-of-Fit  Test 


X 

f (X) 

Exp 

Obs 

(0i~El)  /E 

0 

. 3889 

21.00 

24 

. 4286 

1 

. 3673 

19.83 

16 

.7397 

2 

.1734 

9.36 

10 

. 0438 

3 

.0704 

3.80 

4 

.0104 

1.2225 

or 

more 

x = .9444 

= X 

= 9.210  a=. 01  d.f.=4-l-l=2 

X*  = 1.222<Xg  = 9.21 
2 

Therefore  xs  falls  within  the  acceptance  region  and 
Hq  cannot  be  rejected.  This  result  indicates  that  the  nine- 
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hour  sample  of  departures  was  taken  from  a population  of 
departures  following  a Poisson  distribution  with  a parameter 
X of  . 9444  . 

2 

The  f (x)  value  in  the  x table  was  obtained  using 
>x  -X 

the  formula  f(x)  = e (4:173).  The  computed  f(x) 

represents  the  probability  that  x will  take  on  a specific 
value.  The  probability  value  for  x = 3 or  more  is  the 
cumulative  probability  of  all  10  minute  periods  having  3 
or  more  departures. 

To  confirm  the  conclusion  stated  above,  a second 
hypothesis  was  formulated:  The  mean  times  between  depar- 

tures observed  in  the  nine-hour  period  follow  an  exponential 
distribution. 

H^:  x~Exp  ( X ) 

: Not 

The  random  variable  x is  defined  as  the  mean  time  between 

departures.  The  parameter  X of  the  exponential  distribution 

was  estimated  from  the  sample  mean  time  between  arrivals  (x) . 

The  hypothesis  test  was  resolved  using  the  Kolmogorov-Smirnov 

test  on  maximum  differences  (4:382-386).  The  critical  value 

for  the  maximum  difference  (MAX  D ) was  obtained  from 

c 

Lilliefor's  table  of  computed  values.  The  sample  maximum 
difference  value  (MAX  Ds)  was  computed  as  shown  in  Table 
12,  where: 

x = 10.657 
1/x  = .09383  = X 


60 


Table  12 


Kolmogorov-Smirnov  Test 


Group 

Boundary 

F (x) 

#/ n 

S (x) 

D 

Less  than  . 

5 

. 0488 

0 

0 

. 0488 

.5 

1-2 

2.5 

.2134 

14/15 

.2800 

.0666 

3-4 

4.5 

. 3430 

21/50 

.4200 

*.0770* 

5-6 

6.5 

.4567 

25/50 

. 5000 

. 0433 

7-8 

8.5 

.5507 

29/50 

. 5800 

.0293 

9-10 

10.5 

.6285 

32/50 

.6400 

.0115 

11-12 

12.5 

. 6896 

34/50 

. 6800 

.0096 

13-14 

14.5 

.7433 

37/50 

.7400 

. 0033 

15-16 

16.5 

.7878 

39/50 

. 7800 

. 0078 

17-18 

18.5 

. 8245 

42/50 

. 8400 

.0155 

19-20 

20.5 

.8534 

43/50 

. 8600 

. 0066 

21-22 

22.5 

.8788 

44/50 

. 8800 

.0012 

23-24 

24.5 

. 8997 

44/50 

. 8800 

.0197 

25-26 

26.5 

.9171 

45/50 

. 9000 

. 0171 

27-28 

28.5 

. 9308 

46/50 

. 9200 

. 0108 

More  than 

28.5 

1.0000 

50/50 

1.0000 

0 

MAX  D 

X 

= .0770 

MAX  D 

c 

/ 1.25 
/ N ' 

= .1768 

The 

results 

obtained  from 

Table  12 

indicate  that 

MAX  D falls  within 
s 

the  acceptance 

region 

and  Hq 

cannot  be 

rejected. 

Failure 

to  reject  the  H 

g indicates  the 

sample 
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of  times  between  departures  was  taken  from  an  exponentially 
distributed  population  with  a mean  time  between  departures 
of  10.657  minutes. 


The  F(x)  value  in  the  Kolmogorov-Smirnov  table  was 
obtained  using  the  formula:  F(x)  = l-e-'x  (4:214),  where  x 

is  the  group  boundary.  F(x)  is  the  cumulative  probability 
an  observed  time  between  departures  will  be  less  than  x. 

For  example,  F(14.5)  is  the  probability  an  observed  mean 
time  between  departures  will  be  less  than  14.5 
F ( 14 . 5 ) = l-e-<-°9383)(14.5) 


= l-e-1*360  = 1-.2566  = .7433 


It  was  established  in  the  previous  section  that  air 
operations  are  equally  likely  to  occur  during  any  day  of  the 
week.  It  was  also  established  earlier  in  this  section  that 
air  operations  are  equally  likely  during  the  nine  busiest 
daylight  hours.  Based  on  these  contentions,  the  entire  63 
hours  (9  hours  for  7 days)  were  treated  as  a homogeneous  mass 
and  the  mean  number  of  operations  per  hour  was  determined 
from  the  sample.  The  computation  was  as  follows: 


Day 

AO 

AO 

1 

34 

3.77 

2 

23 

2.55 

3 

31 

3.44 

4 

29 

3.27 

5 

20 

2.22 

6 

25 

2.77 

7 

27 

3.00 

21.00  t 7 = 3 AO/hr 
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Where  AO  = Air  Operation.  AO/hr  was  increased  by  a scaling 
factor  of  14.553  to  give  an  estimate  of  the  mean  number  of 
operations  per  hour  in  the  population. 

14.553  x 3 AO/hr  = 43.66  AO/hr 


The  scaling  factor  is  based  on  the  sample  size  as  a percent- 
age of  the  population. 

4526 


311 


= 14.553 


Figure  9 shows  the  distribution  of  daily  means  after  scaling. 
During  the  63  hours  being  considered,  an  air  operation  occurs 
every  1.374  minutes  or  every  82.4  seconds.  This  number  was 
truncated  and  used  as  the  parameter  lambda  of  an  exponential 
distribution.  In  the  simulation  model,  the  computer  generates 
5000  exponentially-distributed  air  operations  over  a period 
of  409,918  seconds.  These  air  operations  were  assigned  char- 
acteristics in  accordance  with  the  probabilities  previously 
developed.  Consequently,  497  IFR  departures  from  Wright- 
Patterson  were  generated.  This  means  there  were  824  seconds 
or  13.7  minutes  between  departures. 

The  curve-fit  procedures  presented  earlier  in  this 
section  indicated  a mean  time  between  departures  of  10.66 
minutes.  However,  this  value  was  obtained  by  sampling  day 
1 which,  as  shown  above,  had  a sample  mean  number  of 
operations  of  3 .77.  Therefore,  14.553  x 3.77  AO/hr  = 

54.98  AO/hr  or  one  operation  every  1.092  minutes  or  every 
65.5  seconds.  Using  65  for  the  same  5000  generations  gives 


654  seconds  or  10.90  minutes  between  departures.  10.66=10.90. 
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Figure  9.  Aircraft  per  Hour 
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On  the  basis  of  the  foregoing  analysis  it  was  con- 
cluded that  times  between  air  operations  are  exponentially 
distributed  with  a mean  time  between  operation  of  82  seconds. 
Air  activity  was  generated  in  the  simulation  model  in  accor- 
dance with  this  distribution. 

Development 

Morris  described  the  development  of  a model  as  a 
procedure  involving  three  concepts.  The  first  concept  can 
be  usefully  viewed  as  ".  . .a  process  of  enrichment  or 
elaboration  [ll:B-709]."  That  is,  the  model  is  initially 
very  simple  but  is  elaborated  upon  until  it  evolves  into  a 
reflection,  as  close  as  possible,  of  the  actual  situation 
being  modeled. 

The  second  concept  deals  with  the  starting  point  for 
the  elaboration  process.  The  starting  point  should  be  deter- 
mined by  " . . . analogy  or  association  with  previously  well 

developed  logical  structures  . . . [11 : B-709] . " 

The  final  concept  is  an  extension  of  the  first  and 
suggests  that  the  "...  elaboration  or  enrichment  involves 
at  least  two  looping  or  alternation  procedures  [ll:B-709]." 
These  looping  procedures  are  explained  as  follows.  The  first 
loop  occurs  between  modification  of  the  model  and  confronta- 
tion by  the  data.  The  model,  confronted  by  the  data,  is 
tested  and  a modification  results  which  is  again  confronted 
by  the  data. 
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The  second  loop  occurs  between  the  deductive  tract- 
ability  of  the  model  and  the  assumptions  upon  which  it  is 
based.  If  the  model  is  tractable,  one  returns  to  the  assump- 
tions and  elabo. ates  upon  them.  If,  on  the  other  hand,  the 
model  is  not  tractable,  one  returns  to  the  assumptions  and 
simplifies  them.  The  two  looping  procedures  progressively 
refine  both  the  model  and  its  underlying  assumptions  (11: 

B-7  09 ) . 

The  development  of  a simulation  model  for  the  Dayton 
ACA  adhered  to  the  formula  described  by  Morris.  The  start- 

ing point  for  the  model  was  based  on  established  concepts. 
Research  described  earlier  under  Literature  Review  indicates 
simulation  has  been  applied  successfully  in  evaluating  ACA 
configurations.  In  particular,  simulation  based  queuing 
theory  has  been  used  by  Gabrielli  and  Mohleji  (6:1-4;  10: 

1-1  to  1-4) . 

The  initial  model  was  very  simple  and  ignored  many 
aspects  of  the  terminal  environment  that  were  eventually 
included.  It  consisted  of  two  arrival  routes  and  a single 
destination.  Aircraft  flowing  through  this  model  were 
characterized  as  homogeneous  in  nature  and  were  not  distin- 
guished, one  from  another,  by  differences  in  altitude,  air- 
speed, or  ROT.  Furthermore,  variation  in  air  traffic  delay 
times  and  aircraft  departures  were  not  considered  for  the 
initial  model. 
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From  this  starting  point,  the  model  was  elaborated 
through  a stepwise  inclusion  of  the  complicating  factors 
(addressed  below)  until  a basic  scenario  was  developed.  At 
each  step  in  the  elaboration  process,  the  model  was  evaluated 
through  the  application  of  the  two  looping  procedures  recom- 
mended by  Morris  (11:3-709).  The  model  was  modified  and 
tested  repeatedly  at  each  step  to  insure  its  behavior  was 
as  expected.  In  addition,  the  results  obtained  at  each  step 
were  evaluated  to  determine  whether  the  assumptions  support- 
ing the  model  required  elaboration  or  simplification. 

The  first  step  in  the  elaboration  procedure  was  the 
inclusion  of  multiple  arrival  routes.  There  are  six  AFs 
associated  with  Wright-Patter.  son  which  are  "fed"  by  several 
enroute  low  altitude  airways  and  jet  routes  as  illustrated 
in  Figures  3 and  4.  Each  AF  has  associated  with  it  a 
Single  Standard  Vectoring  Route  (SVR) , except  for  ROD  which 
has  two.  The  SVRs  are  illustrated  in  Figure  5. 

The  two  routes  associated  with  ROD  are  required 
because  arrivals  from  that  AF  must  transition  a departure 
corridor.  If  an  aircraft  is  in  the  departure  corridor, 
arrivals  from  ROD  must  maintain  9000  feet  MSL  until  clearing 
the  corridor.  This  requires  the  arrival  aircraft  to  be 
vectored  further  East  (ROD  A)  to  allow  room  for  descent. 

If  the  corridor  is  clear,  or  if  coordination  can  be  arranged 
with  the  departure  controller,  the  arriving  aircraft  can  be 
descended  through  the  corridor  and  therefore  can  be  vectored 
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closer  to  the  FAF  (ROD  B) . Personnel  at  Dayton  Approach 
Control  indicated  the  probability  of  being  cleared  through 
the  departure  corridor  is  .50  (13). 

Aircraft  arriving  at  the  AFs  are  routed  by  an  SVR 
to  the  FAF  and  then  proceed  to  the  runway.  Arrivals  at 
the  AFs  are  generated  by  the  simulation  model  in  accordance 
with  the  distribution  identified  through  analysis  of  flight 
progress  strips. 

The  next  factor  to  be  incorporated  into  the  model  was 
differences  in  airspeeds.  The  performance  characteristics  of 
the  various  aircraft  types  to  be  included  in  the  model  were 
extracted  from  Jane's  All  the  World's  Aircraft  (16:235-462). 
These  characteristics  were  analyzed  to  determine  natural 
groupings  of  aircraft  into  categories,  with  each  category 
possessing  three  airspeeds:  cruise  airspeed,  maneuvering 

airspeed,  and  final  approach  airspeed. 

Five  categories,  "V"  through  "Z,"  were  used. 

Category  V included  high  performance  jet  fighter  type  air- 
craft (examples:  F-4,  T-38,  F-15,  F-100,  etc.)  and  the 

jumbo  jets  (examples:  C-5A,  L-1011,  B-747) . Category  W 

included  the  air-carrier  type  jet  aircraft  that  are  not 
large  enough  to  be  designated  as  jumbo  (examples:  B--707, 

DC-9,  B-737,  DC-8,  C-141,  KC-135) . Category  X included 
the  medium  size  transport  type  aircraft  (examples:  C-130, 

Logair  L-382) . Category  Y included  the  small,  high  per- 
formance type  aircraft  such  as  the  business  jet  (examples: 
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Lear  Jet,  T-39,  King  Air,  T-37,  A-37,  T-33,  OV-IO).  The 
last  category,  Category  Z,  included  the  small,  low  speed 
group  of  aircraft  (examples:  Cessna  172,  180,  Cherokee  140, 

Bonanza) . The  specific  airspeeds  corresponding  to  each 
category  are  listed  under  Data  File  Construction. 

These  airspeeds  are  important  for  the  computations 
of  the  MTF  between  two  points.  The  MTF  varies  for  each 
category,  hence  a separate  computation  for  each  category 
was  required.  The  computed  time  represents  the  average  or 
expected  time  to  fly  between  two  points.  In  the  simulation, 
however,  actual  times  are  stochastic  events  with  as  assumed 
normal  distribution  around  the  MTF.  The  computed  MTF  also 
represents  the  minimum  time  to  fly  between  two  points;  that 
is,  with  no  other  conflicting  traffic  in  the  area,  the  air- 
craft will  be  expected  to  transit  two  points  at  the  MTF. 

Time  greater  than  the  MTF  is  interpreted  as  delay  time, 
and  is  accumulated  by  the  GPSS  computer  program  for  use  in 
evaluating  the  effectiveness  of  the  simulation  model  scenario 
being  tested. 

The  inclusion  of  Discrete  Variable  Delay  Times  ( DVDT) , 
was  the  next  step  in  the  elaboration  process.  Delays  may 
take  three  forms:  holding  patterns,  path  stretching  (vec- 

toring), or  speed  changes.  In  the  case  of  a holding  pattern, 
an  aircraft  is  usually  committed  to  hold  for  five  minutes 
once  it  enters  the  pattern  (the  time  required  to  execute  a 
standard  holding  pattern) . At  the  minimum,  an  aircraft  is 
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committed  for  three  minutes  (the  time  required  to  execute 
a 360  degree  turn)  (18:5-3  to  5-8).  Eecause  of  these 
limitations,  other  delaying  techniques  are  preferred. 

"Path  stretching"  is  a technique  wherein  an  aircraft 
is  directed  along  some  route  other  than  the  SVR.  These 
Alternate  Vectoring  Routes  (AVR)  can  increase  the  time  to 
fly  (delay)  or  decrease  the  time  to  fly  (advance)  between 
two  points.  The  seven  SVRs  used  to  route  aircraft  to  Wright- 
Patterson  offer  very  little  opportunity  to  advance  aircraft 
since  they  cover  the  shortest  distance  required  to  position 
an  aircraft  at  the  FAR.  Because  of  this,  the  basic  scenario 
is  limited  to  delay  vectors.  The  delay  vector  issued  by  the 
controller  is  given  in  anticipation  of  a future  conflict.  In 
this  situation,  the  controller  can  only  estimate  the  amount 
of  delay  required  to  bring  about  the  desired  spacing.  As  a 
result,  the  AVRs  associated  with  each  AF  are  limited  to  two 
or  less  depending  on  the  AF  and  on  airspace  limitations.  The 
AVRs  are  shown  in  Figure  10.  If  the  delay  provided  by  an  AVR 
is  inadequate  to  provide  required  spacing,  the  controller 
will  also  issue  airspeed  instructions. 

Another  delaying  technique  is  airspeed  control.  By 
directing  an  increase  or  a decrease  in  airspeed,  the  con- 
troller can  change  an  aircraft's  arrival  time  at  a given 
point  within  a limited  range.  Airspeed  can  be  incremented 
or  decremented  over  a continuous  scale.  In  automated 
terminal  areas,  airspeed  is  rigidly  controlled  (9:3-1). 
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Figure  10.  Alternate  Vectoring  Routes 


However,  under  normal  circumstances,  airspeed  control  is 
relative;  that  is,  the  controller  will  direct  the  pilot  to 
"reduce  to  final  approach  airspeed"  or  "maintain  present 
airspeed."  The  point  at  which  this  directive  will  be  issued 
is  based  on  the  controller's  estimate  of  what  change  is 
required.  As  a result,  in  the  basic  scenario  airspeed 
reductions  reflect  the  difference  between  maneuvering  air- 
speed (MAS)  and  final  approach  speed  (FAS)  over  a specified 
distance  prior  to  intercepting  the  final  approach  course. 
Similarly,  airspeed  increases  occur  when  the  pilot  main- 
tains MAS  longer  than  normal,  i.e.  until  just  prior  to 
intercepting  final  approach  course. 

The  inclusion  of  DVDTs  in  the  model  is  based  on  the 
characteristics  of  available  delay  techniques  discussed 
above.  Delay  times  resulting  from  these  techniques  are 
accumulated  and  used  as  an  indication  of  scenario  effec- 
tiveness . 


I 
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Variations  in  ROT  was  the  next  factor  to  be  included 
in  the  model.  The  time  an  aircraft  spends  on  the  runway  is 


an  important  variable  in  airport  capacity.  Holland,  in  a 
study  conducted  for  the  Mitre  Corporation,  defines  airport 
capacity  as  ".  . .a  measure  of  the  number  of  operations 
(arrivals,  arrivals  and  departures,  departures)  per  unit  of 
time  [7:2]."  ROT  is  a limiting  factor  in  capacity  because, 
in  general,  simultaneous  use  of  the  runway  by  two  or  more 
aircraft  is  not  permitted.  Once  an  aircraft  is  authorized 
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the  use  of  the  runway,  either  for  takeoff  or  landing,  the 
use  of  the  runway  by  other  aircraft  is  precluded  until  the 
runway  is  once  again  unoccupied.  ROT,  then,  determines  the 
maximum  number  of  operations  the  runway  can  handle  per  unit 
of  time,  and,  consequently,  places  a constraint  upon  the  ATC 
system  capacity. 

ROT  is  not  a significant  variable  in  light  traffic 
situations,  but  its  importance  increases  at  airports  where 
there  exists  a single  runway,  such  as  at  Wright-Patterson  AFB 
In  traffic  situations  wherein  density  approaches  saturation, 
ROT  becomes  a critical  variable  due  to  the  relationship 
between  runway  capacity  and  the  ATC  system  capacity. 

Holland  included  ROT  as  a subset  of  "Aircraft/Pilot 
Performance"  in  his  schema  of  "Factors  Affecting  Airport 
Capacity  [7:5]."  The  ROT  for  a landing  aircraft  begins 
when  the  aircraft  crosses  the  runway  threshold,  just  prior 
to  touching  down,  and  terminates  when  the  aircraft  has 
taxied  physically  clear  of  the  runway.  The  duration  of  the 
ROT  for  a landing  aircraft,  then,  is  a function  of  numerous 
variables  that  include  more  than  pilot  or  aircraft  perfor- 
mance: aircraft  speed  at  touch  down,  aircraft  weight; 

aircraft  braking  capability;  wind  conditions  on  the  runway 
(direction  and  velocity) ; runway  conditions  (icy,  wet,  dry); 
and  the  number  of  available  exits  positioned  along  the  runway 

The  ROT  for  departure  aircraft  begins  when  the  con- 
trol tower  issues  clearance  for  takeoff  and  ends  when  the 
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departing  aircraft  is  airborne.  Factors  that  affect  the  ROT 
of  departure  aircraft  include:  operating  requirements  of 

the  agency  controlling  the  use  of  the  aircraft;  aircraft 
performance  characteristics;  runway  slope;  meteorological 
conditions  (wind  direction  and  velocity,  density  altitude, 
etc . ) . 

In  order  to  incorporate  ROT  into  the  simulation 
model,  data  reflecting  the  duration  of  ROT  by  aircraft 
category  was  obtained.  Observations  were  made  of  departing 
and  landing  aircraft  at  Wright-Patterson  AFB  and  their 
respective  ROTs  timed.  ROTs  for  individual  aircraft  were 
averaged  to  obtain  ROTs  by  category  of  aircraft.  The  cate- 
gory ROTs  derived  through  observation  were  used  to  provide 
general  guidance  only  and  the  computed  values  were  evaluated 
by  tower  personnel  and  modified  as  necessary  to  be  compat- 
ible with  their  experience  (3). 

It  was  determined  that  each  aircraft  category  has  a 
single  average  ROT  that  applies  whether  an  aircraft  in 
the  category  is  landing  or  departing.  The  ROTs  for  the 
various  aircraft  categories  are: 


1. 

Category 

V - 90 

seconds 

2. 

Category 

W - 12 

0 seconds 

3. 

Category 

X - 60 

seconds 

4. 

Category 

Y - 50 

seconds 

5. 

Category 

Z - 20 

seconds 
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These  times  have  been  assigned  to  the  various  aircraft  cate- 
gories in  the  model. 

The  next  step  in  the  elaboration  process  was  the 
inclusion  of  departing  aircraft.  Departures  represent  an 
additional  demand  on  the  runway  and  must  be  sequenced  with 
landing  traffic.  Departing  aircraft  are  generated  by  the 
simulation  in  accordance  with  the  distribution  derived  in 
the  data  analysis.  Departures  normally  have  a lower  priority 
for  use  of  the  runway  than  arrivals  and  are  not  permitted  to 
interfere  with  landing  traffic.  However,  if  several  air- 
craft are  waiting  for  departure,  departure  priority  may  be 
increased.  For  example,  if  five  or  more  aircraft  are  waiting 
for  departure,  the  priority  can  be  increased  to  the  same 
level  as  arrival  aircraft;  if  ten  or  more  are  waiting, 
priority  can  be  promoted  to  a level  higher  than  that  of 
arrivals . 

The  Dayton  ACA  has  four  departure  corridors.  After 
takeoff,  a departing  aircraft  is  no  longer  considered  to  be 
a factor  in  the  simulation,  except  arriving  aircraft  will 
not  be  permitted  to  transit  departure  corridors. 

Multiple  approaches  constituted  the  next  step  in  the 
elaboration  procedure.  Multiple  approaches  include  aircraft 
executing  a Missed  Approach  (MA)  as  well  as  training  flights 
making  a series  of  planned  approaches  and  landings.  An  MA 
occurs  when  the  pilot  of  an  aircraft  on  final  approach 
elects  not  to  land  or  is  denied  permission  to  land  by  the 
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control  tower.  In  this  case  the  arrival  aircraft  is  directed 


back  to  the  FAF  by  way  of  an  appropriate  SVR.  Similarly,  an 
aircraft  making  multiple  approaches  is  treated  essentially 
the  same  as  other  traffic,  the  only  distinction  being  that 
its  SVR  originates  at  the  runway  rather  than  at  an  AF. 

Another  type  of  air  activity  that  falls  into  the  multiple 
approach  category  is  the  touch-an-go.  A touch-and-go 
occurs  when  the  aircraft  lands  and  then  takes  off  again 
without  stopping.  Since  practice  multiple  approaches  and 
touch-and-go  landings  do  not  appear  on  a flight  progress 
strip,  they  were  incorporated  into  the  model  on  the  basis 
of  estimates  by  Dayton  Approach  Control  Personnel  (13) . 

Missed  approaches  caused  by  traffic  conflict  occur  when- 
ever the  simulation  dictates. 

The  seventh  step  in  the  process  of  developing  the 
simulation  model  was  the  allowance  for  instrument  approaches 
other  than  by  radar  vector.  In  executing  an  instrument 
approach,  an  aircraft  flies  to  the  FAF  via  a published  maneu- 
vering procedure  rather  than  an  SVR.  Each  category  of  air- 
craft has  associated  with  it  a time  that  represents  the  MTF 
from  an  IAF  to  the  FAF  via  a published  approach.  Arrival 
time  at  the  FAF  is  considerably  more  variable  for  aircraft 
flying  a published  approach  than  for  aircraft  on  SVRs.  This 
increased  variability  could  cause  delays  to  other  aircraft 
during  high  density  traffic  periods.  The  two  published 
approaches  used  in  the  model  are  shown  in  Figures  11  and  12. 
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HIGH  ALTITUDE 
INSTRUMENT  APPROACH  PROCEDURE 
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Once  cleared  for  an  instrument  approach,  the  pilot  proceeds 
by  his  own  navigation.  Therefore  the  controller  normally 
does  not  control  the  aircraft's  arrival  time  through  vector- 
ing or  airspeed  changes.  Normally,  then,  the  aircraft  will 
be  placed  in  a holding  pattern  until  the  traffic  in  the 
system  provides  an  opening.  These  aspects  of  an  instrument 
approach  have  been  included  in  the  model. 

Another  step  was  the  inclusion  of  traffic  operating 
under  VFR.  Although  VFR  aircraft  are  not  provided  with  the 
same  control  vectoring  or  separation  as  IFR  aircraft,  it  is 
necessary  to  allow  spacing  so  that  VFR  aircraft  can  land. 

In  the  Dayton  ACA  simulation  model,  VFR  aircraft  are  inserted 
into  the  flow  of  IFR  traffic  on  the  basis  of  visual  inspec- 
tion of  the  final  approach  course  by  the  tower  controller. 
Once  inserted,  the  VFR  aircraft  has  a priority  equal  to  IFR 
traffic  and  may  cause  or  incur  delay. 

The  last  factor  to  be  included  in  the  model  was 
differences  in  aircraft  altitudes.  Altitude  is  a third 
dimension,  and  its  inclusion  significantly  complicates  the 
model.  Altitude  separation  allows  faster  aircraft  to  pass 
slower  aircraft  and  permits  simultaneous  occupancy  of  a 
geographic  point  by  more  than  one  aircraft.  The  GPSS 
language  provides  simulation  of  these  characteristics  by 
default  unless  specifically  precluded.  Therefore,  it  was 
only  necessary  to  provide  specific  altitude  separation  in 
the  area  between  a point  of  common  use  (PCU)  and  the  FAF 
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(PCU  is  a point  at  which  two  or  more  SVRs  converge) . The 
area  between  a PCU  and  the  FAF  is  designated  a "facility"  in 
the  simulation  model,  and  therefore  can  only  be  occupied  by 
one  aircraft.  Since  the  traffic  controller  has  the  capabil- 
ity to  maintain  altitude  separation  between  aircraft  in  the 
common  use  areas,  two  "facilities"  were  placed  one  on  top 
of  the  other  along  these  routes.  The  result  is  that  altitude 
separation  is  provided  throughout  an  aircraft's  flight  pro- 
file, from  the  enroute  structure  to  the  FAF.  Figure  13 
presents  a profile  of  the  altitudes  involved. 

The  various  complicating  factors  described  above 
have  been  incorporated  into  the  model  using  GPSS  program 
entities  and  instructions. 

GPSS  Simulation  Program 

Appendix  A contains  a complete  listing  of  the  GPSS 
program  used  to  simulate  the  Dayton  ACA  basic  scenario.  In 
addition,  the  following  information  is  provided  to  assist 
the  reader  in  understanding  the  program. 

Transactions  (aircraft)  which  enter  the  simulation 
are  assigned  parameters  which  reflect  system  status  or 
characteristics  possessed  by  a specific  transaction.  The 
parameters  which  may  be  assigned  to  a transaction  are: 
Floating  Points  (PL^) : 
i Contents 


1. 


Flight  Rules 


Optimum 


Figure  13.  Altitude  Profile 


2.  Category 

3.  Cruise  Airspeed 

4 . Maneuvering  Airspeed 

5.  Final  Approach  Airspeed 

6.  Operation 

7.  Airport 

8.  Route  in  miles 

9.  Alternate  Vectoring  Route  in  miles 

10.  Second  Alternate  Vectoring  Route  in  miles 

14.  Time  to  Fly  one  mile  at  manuevering  airspeed 

15.  Number  of  the  savevalue  which  contains  the 
computed  total  delay  time 

16.  Number  of  facilities  used 

17.  Initial  value  of  split  sequence 

18.  Route  number 

19.  Storage  used 

22.  Queue  number 

23.  Distance  to  high  initial  approach  fix 

24.  Distance  to  low  initial  approach  fix 


25. 

Total 

time 

to 

fly  from  Point  of  Common  Use 

(PCU) 

to  runway 

26. 

Delay 

time 

for 

AVR  number  1 

27. 

Delay 

time 

for 

AVR  number  2 

28. 

Delay 

time 

for 

airspeed  reduction 

29. 

Delay 

time 

for 

holding  pattern 
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30.  Total  time  to  fly  high  altitude  instrument 


| 

* 


I 


approach. 

31.  Total  time  to  fly  low  altitude  instrument 
approach. 

Halfword  (PH^) : 
i_  Contents 

2 . Runway  Occupancy  Time 

3.  Extra  delay  imposed  by  delay  technique. 

4.  High  facility  reference 
6 . Loop  index 

15.  Time  to  fly  one  mile  on  final. 

17.  Split  sequence  number 
Fullword  (PF^) : 
i.  Contents 
1.  Clock  time  at  AF 
10.  Clock  time  at  exit  from  runway. 

15.  Clock  time  at  PCU. 

The  area  from  the  PCU  to  the  runway  consists  of  a 
series  of  19  one-mile  facilities  which  are  labeled  BDFl- 
BDF13  and  HIGH1-HIGH6 . The  HIGH,  facilities  allow  altitude 

l 

separation.  Figure  14  provides  a conceptual  view  of  the 
arrangement  of  facilities.  Other  facilities  include  the 
runway  (RWY) , the  high  altitude  instrument  approach  course 
(HAA) , and  the  low  altitude  instrument  approach  course  (LAA) . 
The  purpose  of  this  arrangement  of  facilities  is  to  provide 
aircraft  separation.  A transaction  seizes,  in  turn,  the 
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Figure  14.  Arrangement  of  GPSS  Entities 


number  of  facilities  matching  the  required  following  distance 
foT~i t s~ ca tego ry . — S±rrcre~  a facility  cran  only~be~  occupied— by- — 
one  aircraft,  the  following  aircraft  is  prevented  from  vio- 
lating the  required  following  distance.  If  it  is  delayed, 
the  delay  time  is  aggregated  and  the  delay  technique  used  is 
identified. 

The  PCU  area  also  includes  two  storages  (storage  10 
and  storage  11) . A storage  may  be  used  simultaneously  by 
more  than  one  transaction.  Their  purpose  is  to  allow  exam- 
ination of  the  included  facilities  to  determine  if  they  are 


in  use.  Departure  entries  to  the  runway  examine  storage  11. 
If  it  is  empty,  they  seize  the  runway.  VFR  entries  examine 
storage  10.  If  it  is  empty,  they  seize  BDF12 . 

The  split  block  is  used  to  create  duplicates  of 
arriving  aircraft.  The  duplicates  enter  a specified  queue 
and  remain  there  for  a time  equal  to  the  computed  delay  for 
the  "parent"  transaction.  A separate  queue  is  assigned  to 
each  type  of  delay  the  "parent"  transaction  may  incur.  Usage 
of  the  various  queues  is  computed  by  GPSS  and  presented  in 
the  standard  output.  To  supplement  the  standard  output,  a 
series  of  matrix  savevalues  were  established  to  record 
the  category  and  route  of  each  transaction  experiencing 
delay  and  the  total  system  time  of  all  transactions.  This 
information  was  used  to  determine  which  areas  of  the  simu- 
lated ACA  offered  potential  for  improvement,  and  to  allow 
comparison  of  total  system  times  during  experimentation. 


85 


— 


Functions  were  used  in  the  program  for  three  purposes. 
-The-  first  was  to  create  the  distribution  of  aircraft  being 
generated  in  the  model.  This  function  (FN1)  uses  a random 
number  as  its  argument  and  produces  exponentially  distributed 
deviates  with  a mean  of  1.0.  This  value  is  multiplied  by  the 
mean  time  between  aircraft  to  provide  exponentially  distrib- 
uted times  between  aircraft. 

The  second  purpose  of  a function  was  to  create  the 
distribution  of  times  to  fly  between  any  two  points.  This 
was  accomplished  through  the  use  of  two  functions.  The 
first  (FN2)  uses  a random  number  as  its  argument  and  produces 
deviates  that  are  normally  distributed  between  5 and  -5  with 
a mean  of  0.  The  resulting  value  is  then  used  as  the  argu- 
ment of  a second  function  (FN14)  which  produces  deviates 
between  .90  and  1.10  with  a mean  of  1.00.  This  value  is 
multiplied  by  the  mean  time  to  fly  between  two  points  to 
provide  normally  distributed  times.  FN15,  when  used  in 
combination  with  FN2,  gives  similar  results  but  with  greater 
variation  in  times. 

The  third  use  of  a function  was  as  the  modifier  in 
a TRANSFER  statement.  A function  used  for  this  purpose 
(FN3-FN13)  uses  a random  number  as  its  argument  and  provides 
a symbolic  program  address  as  its  output.  Each  program 
address  in  the  function  occurs  with  some  specified  prob- 
ability. The  transaction  at  the  TRANSFER  block  is  then 
routed  to  the  program  address  provided  by  the  function. 
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The  probabilities  specified  in  the  functions  are  those 
described  in  the  Analysis  of  Flight  Progress  Strips  section. 

Other  blocks  and  entities  were  used  in  the  simulation 
as  required  to  provide  a logical  flow  of  transactions . The 
programming  is,  in  general,  in  accordance  with  basic  GPSS 
programming  techniques  (8,  14) . One  section,  however, 
requires  additional  explanation.  MACRO  statement  number 
eight  (EIGHT:  MACRO;  a,  b,  c)  provides  the  capability  for 
a transaction  to  reduce  its  MTF  (advance)  to  avoid  conflict. 
This  reduction  simulates  the  air  traffic  controller's 
ability  to  recognize  in  advance  a potential  conflict  and 
to  avoid  that  conflict  by  directing  one  of  the  conflicting 
aircraft  to  maintain  MAS  when  it  would  normally  be  reducing 
to  FAS.  Such  a directive  would  have  the  effect  of  acceler- 
ating one  aircraft  to  pass  another  so  both  could  continue 
without  delay. 

The  series  of  statements  contained  in  MACRO  number 
eight  are  referenced  each  time  an  aircraft  arrives  at  one 
of  the  AFs . Estimated  times  of  arrival  (ETAs)  are  computed 
for  the  aircraft's  arrival  at  the  PCU  using  both  normal  and 
increased  speed.  The  computed  ETAs  are  compared  with  ETAs 
of  other  aircraft  converging  on  the  same  PCU.  If  the  com- 
parison indicates  the  current  aircraft  will  obtain  adequate 
clearance  by  flying  at  the  increased  speed,  then  it  is 
assigned  the  revised  speed.  The  ETA  reflecting  the  change 
is  then  recorded  for  future  comparisons.  If  the  comparison 
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indicates  the  current  aircraft  will  not  obtain  adequate 
clearance  by  flying  at  the  increased  speed,  or  if  it  already 
has  adequate  clearance  at  the  normal  speed,  then  it  retains 
its  normal  speed.  In  this  case  the  ETA  reflecting  normal 
airspeed  is  recorded  for  future  comparisons.  These  compari- 
sons are  accomplished  on  all  other  aircraft  converging  on 
the  same  PCU  up  to  a maximum  of  five.  It  is  felt  that  this 
algorithm  provides  a good  approximation  of  the  effect  on  air 
traffic  of  the  controller's  decision-making  process. 

The  preceding  discussion,  combined  with  a basic 
understanding  of  GPSS  programming  techniques,  and  reference 
to  Appendix  A should  provide  sufficient  understanding  of  the 
techniques  employed  in  the  Dayton  ACA  simulation. 

Validity 

During  the  development  process  the  validity  of  the 
Dayton  ACA  simulation  model  was  evaluated  at  each  step.  The 
evaluation  was  directed  towards  examining  internal  validity 
and  external  validity. 

Internal  validity  of  a research  design  is  its  "... 
ability  to  measure  what  it  aims  to  measure  [5:120],"  or  in 
the  context  of  a simulation  model,  its  ability  to  represent 
what  it  aims  to  represent.  One  measure  of  internal  validity 
is  the  extent  to  which  important  aspects  of  the  situation 
under  study  are  covered  by  the  model.  As  a practical  matter, 
certain  aspects  of  the  Dayton  ACA  were  excluded  from 
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consideration  (see  Assumptions  and  Limitations)  based  on  the 
assumption  that  their  exclusion  would  not  detract  signifi- 
cantly from  the  validity  of  the  model  and  that  their  inclu- 
sion would  represent  unnecessary  complication. 

The  ability  of  the  model  to  represent  what  it  aims 
to  represent  depends  heavily  on  the  data  upon  which  it  is 
based.  The  data  collection  and  analysis  procedures  described 
earlier  were  carefully  developed  to  meet  this  test  of  valid- 
ity. The  aircraft  activity  in  the  fully  developed  model  was 
then  monitored  continuously,  during  a variety  of  runs,  and 
compared  with  the  desired  aircraft  activity  identified  from 
the  sample.  The  fully  developed  model  was  found  to  be  an 
accurate  reflection  of  the  actual  situation  with  respect  to 
its  inputs. 

The  internal  validity  of  the  model  also  depends  on 
its  consistency  and  logic.  Consistency  can  be  tested  by 
using  the  same  scenario  on  subsequent  runs  of  the  simulation 
with  a different  random  number  sequence.  The  stochastic 
events  in  the  simulation  will  occur  in  a different  way,  but 
the  overall  outcome  of  each  run  should  be  equivalent.  The 
Dayton  ACA  was  tested  in  this  manner  and  found  to  be  con- 
sistent. 

Changes  in  characteristics  of  the  model  should  pro- 
vide logically  anticipated  results.  For  example,  the  intro- 

I 

duction  of  variable  ROTs  could  logically  be  expected  to 
produce  increased  delays  for  arriving  aircraft.  Failure 


of  the  model  to  produce  this  expected  result  would  raise 
questions  as  to  its  validity.  During  the  development  of 
the  Dayton  ACA  model,  the  logic  of  the  results  was  evaluated 
at  each  step  in  the  elaboration  procedure.  Illogical  or 
unexpected  results  were  analyzed  and  corrected  before  pro- 
gressing to  the  next  step. 

The  sensitivity  of  the  model  to  changes  in  input 
parameters  (e.g.  the  characteristics  of  arriving  aircraft) 
is  another  measure  of  internal  validity.  For  example,  a 
decrease  in  the  mean  time  between  arrivals  should  result 
in  a proportionate  increase  in  the  delays  experienced  by 
traffic  in  the  system.  Similarly,  a decrease  in  the  per- 
centage of  arriving  aircraft  of  one  category  (decreased 
homogeneity  of  traffic)  should  result  in  a proportionate 
increase  in  delays  experienced  by  traffic  in  the  system. 
Various  times  between  arrivals  were  used  to  test  the  models 
sensitivity.  The  effect  on  system  traffic  was  proportionate 
to  the  change  in  time  between  arrivals,  and  was  consistent 
with  anticipated  results.  Division  of  arriving  aircraft 
into  five  categories  resulted  in  increased  average  delay 
to  system  traffic,  particularly  to  high  speed  aircraft. 

This  result  is  as  would  be  logically  expected. 

Times  to  fly  between  two  points  have  been  assumed 
to  be  normally  distributed  around  the  MTF  between  those 
points.  The  variance  of  the  normal  distribution  was  set  at 
several  values  to  test  the  sensitivity  of  the  system  to  this 
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parameter.  The  effect  was  found  to  be  negligible  when  the 
nominal  time  between  air  activity  of  82  seconds  was  used. 

This  result  is  attributed  to  the  fact  that  the  final  approach 
course  at  Wright-Patterson  is  operating  far  below  capacity. 

As  a result,  the  occurrence  of  a conflict,  and  therefore  of 
a delay,  is  a function  of  the  randomness  of  air  traffic. 

This  randomness  is  not  effected  by  the  stochastic  nature  of 
times  to  fly  between  various  points  in  the  ACA. 

The  internal  validity  of  the  Dayton  ACA  simulation 
model  has  been  adjudged  to  be  satisfactory  based  on  the 
criteria  discussed  above.  Further  tests  were  conducted  on 
the  models  external  validity. 

External  validity  is  concerned  with  the  degree  to 
which  findings  of  the  model  can  be  generalized  to  the  real 
world  situation  that  it  represents  (5:121).  External 
validity  can  be  determined  by  comparison  of  model  output 
with  observed  occurrences  in  the  real  world,  and  by 
evaluation  of  model  output  by  experts.  The  output  from 
the  Dayton  ACA  model  was  evaluated  by  these  two  methods. 
Actual  traffic  flow  was  observed  and  compared  to  simulated 
traffic  flow  in  the  model.  In  addition,  the  chief  of  the 
Dayton  Radar  Approach  Control  and  his  assistant  were  asked 
to  evaluate  the  model  at  various  stages  in  its  development 
and  to  compare  the  model's  behavior  with  actual  traffic 
flow  (3,  13).  Thus,  the  external  validity  of  the  model 
was  confirmed.  However,  the  ultimate  validity  of  the 
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model  and  any  solutions  generated  by  it  can  only  be  con- 
firmed by  field  test  in  the  real  world  situation.  Such  a 

v 

test  is  beyond  the  scope  of  this  research. 


Application  of  the  Model 

Validation  of  the  computer-simulation  model  of  the 
Dayton  ACA  marked  the  completion  of  phase  one  as  established 
in  the  objective  section.  The  validated  model  was  the  basic 
tool  needed  to  conduct  further  research  through  experimenta- 
tion with  changes  in  ACA  design  and  ATC  procedures.  The 
computer  output  of  the  validated  model,  presented  in  Chapter 
4,  is  the  reference  or  standard  against  which  all  experi- 
ments were  compared  for  evaluation. 

Analysis  of  the  output  of  the  validated  model  sug- 
gested three  general  areas  in  which  to  experiment.  The 
first  two  areas  involve  changes  in  ATC  procedures,  and  the 
third,  ACA  design:  (1)  service  priority  bared  on  AFs, 

(2)  service  priority  based  on  aircraft  speed  category,  and 

(3)  arrival  route  re-design. 

Service  priority  based  on  AF.  The  computer  output  of  the 
validated  model  indicated  wide  disparity  in  the  utilization 
of  AFs  by  arrival  aircraft  destined  for  Wright-Patterson . 

For  example,  approximately  45%  of  Wright-Patterson  bound 
aircraft  enter  the  ACA  at  OPW  (see  Figure  5) . The  queue 
discipline  used  in  the  validated  model  is  "first  come, 
first  served,"  and  no  service  priority  is  given  to  arrival 
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aircraft  except  as  indicated  by  MACRO  statement  number  eight 
(see  GPSS  Simulation  Model) . Assignment  of  service  priority, 
as  a function  of  AF,  was  expected  to  decrease  average  system 
time  (the  average  time  from  the  arrival  of  an  aircraft  into 
the  system  at  an  AF  to  its  departure  from  the  system  at  a 
runway  exit) . 

This  speculation  represented  a departure  from  the 
normal  ATC  procedure  of  handling  air  traffic  on  a first 
come,  first  served  basis.  Experimentation  with  queue 
discipline,  however,  was  justified  on  the  basis  of  its 
potential  to  illuminate  the  dynamic  interaction  of  air- 
craft in  the  system.  Following  this  line  of  reasoning, 
several  experiments  were  conducted  in  which  the  validated 
model's  queue  discipline  was  changed  from  first  come,  first 
served  to  prioritized  discipline  based  on  AFs.  The  experi- 
ments, and  the  results,  are  presented  in  Chapter  4. 

Service  priority  based  on  speed  category.  Aircraft  speed 
categories  is  another  area  wherein  use  of  priority  offered 
potentially  rewarding  experimentation.  Assignment  of  service 
priority  to  aircraft  on  the  basis  of  speed  category  was 
hypothesized  to  be  a means  of  reducing  average  system  time 
because  of  the  number  of  high  speed  aircraft  arriving  at 
Wright-Patterson.  Experiments  were  conducted  to  test  the 
hypothesis  by  altering  the  validated  model  to  assign  service 
priority  by  speed  category  in  accordance  with  several  unique 
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arrangements.  These  experiments  explored  the  effect  of 
changing  the  standard  ATC  "first  come,  first  served"  pro- 
cedure as  did  the  experiments  dealing  with  service  priority 
by  AF.  The  experiments  and  results  are  also  reported  in 
Chapter  4. 

Arrival  route  re-design.  The  last  experiment  accomplished 
with  the  simulation  model  was  in  the  area  of  arrival  route 
re-design.  The  design  considered  for  experimentation  was 
a single  straight-in  arrival  route  to  the  runway  in  lieu  of 
the  existing  route  structure.  The  validated  model  was 
restructured  to  permit  all  aircraft  entering  the  ACA  at 
the  six  AFs  to  converge  on  a single  navigational  fix  16 
nautical  miles  from  the  approach  end  of  Wright-Patterson ' s 
runway  23.  All  aircraft  proceed  from  this  fix  to  the  run- 
way via  a straight-in  approach  course.  This  experiment 
did  not  represent  a radical  departure  from  the  existing 
route  structure  (all  seven  arrival  routes  in  the  validated 
model  eventually  combine  to  form  a single  straight-in  course 
as  depicted  in  Figure  5) , but  it  did  probe  the  potential  of 
ACA  re-design  to  reduce  average  system  time  This  experi- 
ment is  reported  in  Chapter  4 including  a detailed  descrip- 
tion of  the  design  change  and  experiment  results. 
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Chapter  4 


RESULTS 

This  chapter  describes  the  results  of  several  experi- 
ments conducted  on  the  basic  simulation  model.  The  purpose 
of  these  experiments  was  to  improve  understanding  of  air- 
craft interaction  in  the  Dayton  ACA.  Three  groups  of  experi- 
ments were  conducted:  (1)  service  priority  based  on  AF, 

(2)  service  priority  based  on  aircraft  speed  category,  and 

(3)  arrival  route  re-design.  The  computer  output  from  the 
basic  model  will  be  described  first,  followed  by  a discussion 
of  each  group  of  experiments.  The  chapter  will  conclude  with 
a summary  of  experiment  results. 

Computer  Output  of  the  Basic  Model 

The  computer  output  from  a GPSS  simulation  is  too 
extensive  to  be  herein  presented.  The  significant  informa- 
tion has  been  extracted,  summarized  and  presented  in  Tables 
13  and  14.  Table  13  summarizes  the  results  obtained  with 
RNl  (random  number  sequence  number  1) , and  Table  14  sum- 
marizes the  results  obtained  with  RN5.  The  notation  used 
in  the  tables  is  explained  as  follows:  D = average  delay, 

<Tr  = average  system  time,  D = overall  average  delay,  and 
IT  = overall  average  system  time.  MA(UP)  is  the  number  of 
unplanned  missed  approaches,  and  MA(P)  is  the  number  of 
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Table  13 


Model  Data:  Basic  (RN1) 


TOTAL  ARRIVALS:  493 

TOTAL  DEPARTURES:  450 

MA(UP)  : 13 

MA(P) : 70 


ROUTE 

# 

D 

ST 

CATEGORY 

# 

ST 

DEP 

450 

51.9 

— 

V 

32 

723.4 

VFR 

26 

28.0 

133.2 

W 

142 

815.8 

OPW 

170 

28.3 

880.0 

X 

43 

934.8 

RODA 

15 

65.0 

1596.0 

Y 

68 

1203.1 

RODB 

19 

38.3 

1005.9 

Z 

182 

1479.2 

RID 

7 

40.9 

2175.3 

70U 

62 

32.5 

1080.9 

D = 

97.23 

3ZH 

13 

31.7 

1468.2 

EAT 

64 

34.5 

1764.4 

ST  = 

1082. 

50 

MA 

83 

545.0 

729.7 

HAA 

16 

321.7 

2032.9 

LAA 

18 

321.6 

1678.4 

Table  14. 

Model  Data:  Basic  (RN5) 


TOTAL  ARRIVALS:  520 

TOTAL  DEPARTURES:  415 

MA  (UP)  : 11 

MA  (P)  : 71 


ROUTE 

« 

D 

ST 

CATEGORY 

# 

ST 

DEP 

415 

74.0 

_ 

V 

29 

773.8 

VFR 

45 

71.8 

177.9 

W 

130 

861.5 

OPW 

163 

26.6 

892.8 

X 

66 

927.6 

RODA 

18 

59.2 

1517.5 

Y 

68 

1174.0 

RODB 

15 

34.7 

958.4 

Z 

183 

1472.3 

RID 

6 

33.3 

2038.8 

70U 

74 

36.8 

1112.8 

D = 

114. 

6 

3ZH 

17 

29.6 

1407.4 

EAT 

52 

35.3 

1783.2 

ST  = 

1063 

.3 

MA 

82 

555.3 

730.8 

HAA 

16 

313.6 

1971.1 

LAA 

33 

347.8 

1666.5 
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planned  missed  approaches.  Comparison  of  Tables  13  and  14 
reveals  variation  in  the  number  of  aircraft  in  each  category 
and  in  the  number  of  aircraft  that  used  each  route.  Because 
of  this  variation,  experimental  models  will  be  compared  with 
the  basic  model  which  used  the  same  random  number  sequence. 
The  basic  model  is  the  "standard,"  with  which  all  experiments 
were  compared. 

Experiments  Involving  Priority  by  AF 

The  first  experiment  revised  the  priority  of  aircraft 
arriving  from  OPW,  70U,  and  EAT.  The  geographic  locations  of 
these  three  AFs  are  shown  in  Figure  5,  page  29.  Aircraft 
arriving  from  OPW,  70U,  and  EAT  were  assigned  PR6  (priority 
level  six),  PR5,  and  PR4,  respectively,  instead  of  their 
normal  PR3 . These  three  AFs  account  for  79  percent  of  the 
aircraft  arriving  at  Wright-Patterson  AFB . It  was  antici- 
pated that  increasing  their  priority  would  reduce  delays  in 
the  system. 

The  computer  output  from  this  experiment  is  summa- 
rized in  Tables  15  and  16.  The  data  in  Table  15  resulted 
from  using  RN1,  while  the  data  in  Table  16  resulted  from 
using  RN5.  The  output  from  this  model  was  compared  with  the 
output  from  the  basic  model  in  an  attempt  to  identify  pat- 
terns, directions,  and  magnitudes  of  change.  Comparison  of 
D and  ST  reveals  no  consistent  pattern  of  change.  The 
inconsistent  behavior  of  D and  ST  was  assumed  to  be  partially 
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Table  15 


Model  Data:  0PW=PR6/70U=PR5/EAT=PR4 (RNl) 


TOTAL 

ARRIVALS : 

497 

TOTAL 

DEPARTURES 

: 4 36 

MA(UP) 

: 18 

MA ( P ) : 

78 

ROUTE 

# 

D 

ST 

CATEGORY 

# 

ST 

DEP 

436 

57.5 

_ 

V 

42 

726.3 

VFR 

42 

79.2 

183.2 

W 

144 

793.1 

OPW 

147 

26.8 

887.8 

X 

32 

978.5 

RODA  ' 

17 

54.5 

1394.0 

Y 

63 

1272.1 

RODB 

24 

42.0 

974.0 

Z 

174 

1410.2 

RID 

10 

45.1 

2064.2 

70U 

68 

35.8 

1099.5 

D = 

107.6 

3ZH 

14 

37.0 

1539.3 

EAT 

55 

37.3 

1718.5 

ff  = 

1024  . 

6 

MA 

96 

556.0 

720.8 

HAA 

7 

300.0 

1847.9 

LAA 

17 

307.0 

1784.2 

Table  16. 

Model  Data:  0PW=PR6/70U=PR5/EAT=PR4 (RN5) 


TOTAL 

ARRIVALS : 

497 

TOTAL 

DEPARTURES 

: 432 

MA  ( UP ) 

: 19 

MA  (P)  : 

66 

ROUTE 

# 

D 

ST 

CATEGORY 

# 

ST 

DEP 

432 

58.4 

— 

V 

38 

787.4 

VFR 

28 

29.8 

135.6 

W 

155 

833.1 

OPW 

175 

27.4 

896.3 

X 

35 

1022.3 

RODA 

19 

63.8 

1530.1 

Y 

56 

1284.6 

RODB 

19 

60.1 

1106.9 

Z 

185 

1400.3 

RID 

7 

48.4 

2067.4 

70U 

57 

35.0 

1101.2 

D 

= 106. 

4 

3ZH 

14 

34.5 

1482.8 

EAT 

54 

34.4 

1705.6 

ST 

= 1065 

.6 

MA 

85 

560.9 

734.5 

HAA 

16 

314.3 

1835.1 

LAA 

23 

357.3 

1607.7 
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caused  by  changes  in  the  volume  of  traffic  on  each  route,  and 
the  volume  in  each  category.  These  changes  were  caused  by 
the  PRIORITY  block  inserted  into  the  program  to  achieve  the 
desired  priority  alterations.  The  new  PRIORITY  blocks 
re-initiate  the  GPSS  scan  of  the  current  events  chain  (8:3- 
21/22) . The  result  is  a change  in  the  order  of  arrival  of 
transactions  at  various  blocks  in  the  program  and  a change 
in  the  characteristics  assigned  to  each  transaction.  In  an 
effort  to  reduce  the  effect  of  these  changes  in  the  volume  of 
traffic,  ST  by  route,  D by  route,  and  ST  by  category  were 
compared  with  the  counterpart  values  in  the  basic  model. 

These  comparisons  are  portrayed  graphically  in  Figures  15 
through  20.  In  each  figure,  the  solid  line  is  the  basic 
model  and  the  dashed  line  is  the  experimental  model.  Again, 
there  is  no  strong  pattern  of  change,  and  the  changes  that 
occur  are  inconsistent. 

Analysis  was  conducted  to  determine  the  relative 
efficiency  of  the  model.  An  efficiency  factor  was  developed 
which  takes  into  account  the  change  in  D by  route,  by  volume 
of  traffic  on  each  route,  and  by  mix  of  aircraft  categories 
of  traffic  on  each  route.  The  efficiency  factor  was  com- 
puted for  each  experimental  model  as  follows: 


12  G 

E - l 1,1 

i-1 


d.  f . 
i i 


a . 

L ij 


where : 


Route  number 
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BASIC  MODEL:  SOLID  LINE 

EXPERIMENTAL  MODEL:  DASHED  LINE 
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Figure  15:  Average  System  Time  by  Route  - Experimental  Model 

(0PW=PR6/70U=PR5/EAT=PR4-RN1)  versus  Basic  Model 


BASIC  MODEL:  SOLID  LINE 

EXPERIMENTAL  MODEL:  DASHED  LINE 
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Figure  16.  Average  System  Time  by  Route  - Experimental  Model 
(0PW=PR6/70U=PR5/EAT=PR4-RN5)  versus  Basic  Model 


BASIC  MODEL:  SOLID  LINE 

EXPERIMENTAL  MODEL:  DASHED  LINE 


(0PW=PR6/70U=PR5/EAT=PR4-RN1)  versus  Basic  Model 


BASIC  MODEL:  SOLID  LINE 
EXPERIMENTAL  MODEL:  DASHED  LINE 


Figure  18.  Average  Delay  by  Route  - Experimental  Model 

(0PW=PR6/70U=PR5/EAT=PR4-RN5)  versus  Basic  Model 


BASIC  MODEL:  SOLID  LINE 

EXPERIMENTAL  MODEL:  DASHED  LINE 
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0*5  >-  ui »—  us  : 


cL  = (D\,  Experimental)  - (D^,  Basic) 

= Number  of  aircraft  using  route  i 

a^  = Adjustment  for  category  mix 

E>0  is  considered  unfavorable,  and  E<0  is  considered  favor- 
able . 

The  adjustment  factor  for  category  mix  was  computed 
for  each  route  as  follows: 


a . = k 
1 


3 5 

c . + y c . 
j-i  3 3=4  ^ 


where:  j 


c ■ 


3 


k 


Category  number 

(Number  of  category  j aircraft  using  route 
i,  Experimental)  - (Number  of  category  j 
aircraft  using  route  i,  Basic) 

100 


If  the  number  of  aircraft  in  each  category  on  route  i is  the 
same  as  in  the  basic  model,  then  a^  = 100.  If  the  aircraft 
on  route  i tended  to  be  faster  (i.e.,  more  were  in  category 
1,  2,  or  3),  a^<100  and  E is  increased.  If  the  aircraft  on 
route  i tended  to  be  slower  (i.e.,  more  were  in  category  4 
or  5) , a^>100  and  E is  decreased. 

The  computation  of  E can  be  demonstrated  best  with 
an  example.  The  following  data  was  extracted  from  the 
computer  output: 
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Model : 

Basic (RNl) 

0PW=PR6/70U=PR5/EAT=PR4 (RNl) 

Route : 

OPW 

OPW 

Total  Number: 

170 

147 

Category  1: 

10 

17 

Category  2: 

60 

50 

Category  3: 

15 

8 

Category  4: 

15 

12 

Category  5: 

75 

60 

D: 

28.3 

26.8 

al  = 100- [ (10-17)  + (60-50)  + (15-8) ] + [ (15-12)  + (75-60) ]=108 
dx  = 2 6.8-28. 3=-l . 5 
fx  = 170 

E(OPW)  = = -2.36 

E (Route)  is  computed  in  a similar  manner  for  each  route,  and 
the  sum  of  these  values  is  E.  The  value  of  the  efficiency 
factor  for  the  first  experiment  is  +61.03  using  RNl , and 
-64.34  using  RN5.  The  efficiency  factor  indicates  an 
unfavorable  change  in  the  first  case,  but  a favorable  change 
in  the  second  case. 

The  second  experiment  increased  the  priority  of  air- 
craft arriving  from  OPW  from  PR3  to  PR5.  OPW  accounts  for 
45  percent  of  the  arrivals  to  Wr ight-Patterson  AEB,  and  it 
was  believed  that  giving  these  aircraft  service  priority 
would  improve  the  efficiency  of  the  system.  The  computer 
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output  from  this  experiment  is  summarized  in  Tables  17  and  18. 
A comparative  analysis  was  conducted  on  this  data  in  the  same 
manner  as  described  for  the  previous  experiment.  Changes  in 
D and  ST  were  again  inconsistent  and  showed  no  obvious  pat- 
tern. Analysis  of  ST  by  route,  D by  route,  and  ST  by 
category  provided  equally  inconsistent  results.  The  value 
of  the  efficiency  factor  was  found  to  be  +61.73  using  RNl 
and  -86.79  using  RN5. 

The  final  experiment  in  this  group  was  effected  by 
revising  the  priority  of  aircraft  arriving  from  OPW,  70U, 
and  EAT  from  PR3  to  PR5.  The  rationale  used  for  this  experi- 
ment was  that  the  hierarchy  of  priorities  used  in  the  first 
experiment  might  be  causing  complex  traffic  interactions , 
thus  creating  delays.  The  computer  output  from  this  experi- 
ment is  summarized  in  Tables  19  and  20.  The  same  analysis 
was  conducted  on  this  data  as  on  the  previous  experiments. 

As  before,  there  is  no  discernable  pattern  in  the  output. 

The  value  of  the  efficiency  factor  is  +88.44  using  RNl,  and 
-46.09  using  RN5 . This  concluded  the  experiments  involving 
priority  by  AF. 

Experiments  Involving  Priority  by  Speed  Category 

The  first  experiment  in  this  group  tested  a hierarchy 
of  service  priority  for  the  five  aircraft  speed  categories 
(V  through  Z) . The  aircraft  types  and  airspeeds  associated 
with  each  category  are  presented  in  Data  File  Construction, 
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Table  17 


1 


f 


Model 

Data:  OPW= 

:PR5  ( RN1 ) 

TOTAL 

ARRIVALS : 

505 

TOTAL 

DEPARTURES 

: 437 

MA  ( UP 

) : 11 

MA  (P ) 

: 74 

ROUTE 

# 

D 

ST 

CATEGORY 

# 

ST 

DEP 

437 

61.5 

_ 

V 

30 

732.1 

VFR 

42 

57.6 

162.6 

W 

143 

854.3 

OPW 

172 

26.7 

902.4 

X 

47 

899.1 

RODA 

20 

55.2 

1412.4 

Y 

74 

1168.6 

RODB 

19 

42.8 

1018.3 

Z 

179 

1405.1 

RID 

11 

45.7 

2289.4 

70U 

64 

36.6 

1106.9 

D 

= 105. 

9 

3ZH 

13 

28.7 

1434.6 

— _ _ 

EAT 

41 

35.3 

1684.4 

ST 

= 1037 

.4 

MA 

85 

550.8 

737.2 

HAA 

10 

354.3 

2140.4 

LAA 

28 

320.6 

1652.6 

i 

I 


Table  18. 


Model  Data:  OPW=PR5(RN5) 


TOTAL 

ARRIVALS : 

501 

TOTAL 

DEPARTURES 

: 439 

MA  (UP) 

: 16 

MA  (P)  : 

75 

ROUTE 

# 

D 

ST 

CATEGORY 

# 

ST 

DEP 

439 

60.9 

_ 

V 

35 

720.4 

VFR 

31 

23.8 

129.0 

W 

135 

802.1 

OPW 

166 

26.9 

877.3 

X 

50 

976.4 

RODA 

22 

59.1 

1502.2 

Y 

68 

1171.8 

RODB 

21 

41.3 

989.2 

Z 

182 

1465.9 

RID 

9 

40.0 

2014.6 

70U 

58 

34.5 

1101.3 

D 

= 109. 

8 

3 ZH 

15 

40.0 

1405.5 

____ 

• 

EAT 

55 

36.7 

1788.8 

ST 

= 1063 

.5 

MA 

91 

574.9 

744.8 

I 

HAA 

9 

384.4 

2202.8 

• 

LAA 

24 

346.2 

1678.1 
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Table  19 


Model  Data:  OPW,  70U,  EAT  = PR5(RN1) 

TOTAL  ARRIVALS:  541 

TOTAL  DEPARTURES:  420 

MA(UP) : 18 

MA (P) : 93 


ROUTE 

# 

D 

ST 

CATEGORY  # 

ST 

DEP 

420 

58.6 

_ 

V 

44 

729.3 

VFR 

46 

96.6 

201.2 

W 

137 

796.8 

OPW 

161 

27.7 

883.3 

X 

56 

872.7 

RODA 

15 

61.3 

1558.7 

Y 

68 

1169.1 

RODB 

30 

48.4 

1033.9 

Z 

191 

1454.3 

RID 

10 

37.2 

2015.6 



70U 

73 

37.3 

1145.8 

D = 117. 

9 

3ZH 

13 

53.6 

1608.5 

EAT 

53 

33.9 

1759.5 

If  = 1028 

.9 

MA 

111 

568.0 

738.3 

HAA 

9 

300.0 

1840.7 

LAA 

20 

308.2 

1675.9 

Table 

20. 

Model  Data 

: OPW, 

70U , EAT 

= PR5 (RN5) 

TOTAL 

ARRIVALS : 

510 

TOTAL 

DEPARTURES 

: 412 

MA  ( UP ) 

: 21 

MA  (P)  : 

72 

ROUTE 

# 

D 

ST 

CATEGORY 

# 

ST 

DEP 

412 

51.8 

— 

V 

49 

677.7 

VFR 

21 

176.8 

281.8 

W 

126 

801.7 

OPW 

181 

28.6 

883.9 

X 

71 

848.1 

RODA 

22 

62.8 

1532.8 

Y 

49 

1324.8 

RODB 

20 

48.1 

1040.4 

Z 

194 

1442.8 

RID 

6 

39.1 

1894.3 

70U 

70 

32.5 

1100.1 

D 

= 105. 

0 

3ZH 

15 

35.1 

1361.7 

_ 

EAT 

53 

35.4 

1847.6 

ST 

= 1069 

.0 

MA 

93 

544.5 

716.0 

HAA 

9 

309.1 

1898.6 

LAA 

20 

292.9 

1715.7 
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page  34,  and  Development,  page  64.  The  highest  priority,  PR7 , 
was  assigned  to  category  V aircraft  and  each  of  the  remaining 
four  categories  received  successively  lower  service  priorities 
(W=PR6,  X=PR5,  Y=PR4 , Z=PR3) . The  hierarchy  was  established 
in  response  to  the  hypothesis  that  average  system  time  would 
decrease  if  higher  speed  aircraft  were  always  given  service 
priority.  The  experiment  was  conducted  to  test  the  hypothesis 
and  to  ascertain  the  extent  of  traffic  interaction  by  iden- 
tifying the  magnitude  and  direction  of  change  patterns. 

The  computer  output  from  this  experiment  is  summa- 
rized in  Tables  21  and  22.  The  data  in  Table  21  resulted 
from  RN1  and  that  in  Table  22  from  RN5.  The  data  tables 
were  compared  against  their  counterparts  (Tables  13  and  14) 
of  the  basic  model.  ST  for  all  aircraft  decreased  somewhat 
in  the  experimental  model  but  D increased.  Category  W, 
Category  X,  and  Category  Z aircraft  experienced  a decrease 
in  ST  while  Category  Y experienced  an  increase.  The  change 
in  ST  for  Category  V was  inconsistent. 

Figures  21  and  22  are  a graphic  presentation  of  ST 
by  route,  for  the  experimental  model,  compared  to  the  same 
for  the  basic  model.  Figures  23  and  24  illustrate  D by 
route  for  the  experimental  model,  compared  to  D by  route 
for  the  basic  model.  Figures  25  and  26  graphically  compare 
ST  by  category,  for  the  experimental  model,  to  ST  by 
category  for  the  basic  model.  In  each  figure,  the  solid 
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Table  21 


Model  Data:  V=PR7/W=PR6/X=PR5/Y=PR4/Z=PR3 ( RN1 ) 


TOTAL  ARRIVALS:  504 

TOTAL  DEPARTURES:  435 

MA (UP) : 20 

MA(P) : 86 


ROUTE 

# 

D 

ST 

CATEGORY 

# 

ST 

DEP 

435 

56.2 

— 

V 

43 

723.8 

VFR 

31 

54.5 

161.3 

W 

140 

778.2 

OPW 

149 

29.4 

901.6 

X 

60 

908.0 

RODA 

17 

50.9 

1311.4 

Y 

63 

1212.4 

RODB 

29 

36.5 

892.9 

Z 

167 

1437.2 

RID 

11 

40.4 

2108.2 

__ 

70U 

61 

31.7 

1061.5 

D 

= 109. 

9 

3ZH 

12 

32.9 

1435.8 

EAT 

54 

31.4 

1725.5 

ST 

= 1023 

.7 

MA 

106 

521.9 

697.7 

HAA 

7 

342.9 

1993.0 

LAA 

28 

307.5 

1530.8 

Table  22. 

Model  Data:  V=PR7/W=PR6/X=PR5/Y=PR4/Z=PR3 ( RN5 ) 

TOTAL  ARRIVALS:  490 

TOTAL  DEPARTURES:  406 

MA  (UP)  : 31 

MA (P) : 66 


ROUTE 

# 

D 

ST 

CATEGORY 

* 

ST 

DEP 

406 

76.5 

V 

40 

729.5 

VFR 

33 

63.4 

167.9 

W 

131 

835.8 

OPW 

170 

30.6 

917.2 

X 

53 

918.3 

RODA 

17 

53.6 

1373.2 

Y 

58 

1211.7 

RODB 

17 

55.8 

1114.8 

Z 

208 

1445.7 

RID 

9 

42.3 

2111.6 

_ 

70U 

70 

35.8 

1113.5 

D 

= 128. 

5 

3ZH 

13 

31.2 

1477.9 

EAT 

60 

38.7 

1890.7 

ST 

= 1150 

.7 

MA 

97 

581.5 

754.9 

HAA 

12 

279.8 

1553.5 

LAA 

24 

399.5 

1530.6 
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BASIC  MODEL:  SOLID  LINE 

EXPERIMENTAL  MODEL:  DASHED  LINE 


ROUTE. 

Average  System  Time  by  Route  - Experimental  Model 
(V=PR7/W=PR6/X=PR5/Y=PR4/Z=PR3-RN1)  versus  Basic  Model 


BASIC  MODEL:  SOLID  LINE 

EXPERIMENTAL  MODEL:  DASHED  LINE 


Figure  22.  Average  System  Time  by  Route  - Experimental  Model 

( V=PR7/W=PR6/X=PR5/Y=PR4/Z=PR3-RN5 ) versus  Basic  Model 


BASIC  MODEL:  SOLID  LINE 

EXPERIMENTAL  MODEL:  DASHED  LINE 


Figure  23.  Average  Delay  by  Route  - Experimental  Model 

(V=PR7/W=PR6/X=PR5/Y=PR4/Z=PR3-RNl)  versus  Basic  Model 


BASIC  MODEL:  SOLID  LINE 
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ROUTE 

Figure  24.  Average  Delay  by  Route  - Experimental  Model 

(V=PR7/W=PR6/X=PR5/Y=PR4/Z=PR3-RN5)  versus  Basic  Model 
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BASIC  MODEL:  SOLID  LINE 

ISOQ  r EXPERIMENTAL  MODEL:  DASHED  LINE 


igure  26.  Average  System  Time  by  Category  - Experimental  M 
(V=PR7/W=PR6/X=PR5/Y=PR4/Z=PR3-RN5)  versus  Basic 


line  is  the  basic  model  and  the  dashed  line  is  the  experi- 
mental model. 

The  efficiency  factor  for  the  experimental  model 
using  RN1  is  -8.73.  The  experimental  model  using  RN5  has  an 
efficiency  factor  of  +50.96.  Inspection  of  the  data  failed 
to  identify  a pattern  in  the  magnitude  and  direction  of 
changes  in  D,  D,  ST,  and  ST  attributable  to  the  experimental 
mode 1 . 

The  second  experiment  in  this  group  tested  a struc- 
turing of  service  priority  by  airspeed  category  based  on  a 
grouping  of  aircraft  categories.  Categories  V,  W and  X 
were  combined  in  one  group  and  assigned  PP.5.  Category  Y 
and  Z constituted  a second  group  to  which  PR  3 was  assigned. 
The  basic  model  was  altered  to  reflect  this  new  priority 
schema  and  the  experiment  was  conducted  using  first  RN1  and 
then  RN5. 

The  hypothesis  underlying  this  experiment  was  that 
ST  would  decrease  at  the  expense  of  Categories  Y and  Z (i.e. 
ST  for  Categories  Y and  Z would  increase  while  ST  for  Cate- 
gories V,  W and  X decreased) . Tables  23  and  24  present  the 
computer  output  for  this  experiment.  Table  23  is  the  com- 
puter output  using  RN1 ; Table  24  is  the  output  using  RN5 . 

Tables  23  and  24  were  compared  to  Tables  13  and  14, 
respectively,  to  ascertain  the  nature  of  the  changed  in  D, 
S7?,  D,  and  ST  produced  by  the  experimental  model.  ST  for 
Categories  V,  X,  and  Z decreased  in  the  experimental  model 
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Table  23 


Model  Data:,  V,W,  X=PR5/Y,  Z=PR3  (RNl) 


TOTAL 

ARRIVALS : 

507 

TOTAL 

DEPARTURES 

: 420 

MA  ( UP ) 

: 20 

MA  ( P ) : 

86 

ROUTE 

# 

D 

ST 

CATEGORY 

# 

ST 

DEP 

420 

59.8 

— 

V 

35 

639.3 

VFR 

33 

50.0 

155.6 

W 

111 

855.6 

OPW 

163 

28.7 

895.4 

X 

60 

816.3 

RODA 

17 

62.1 

1542.9 

Y 

77 

1154.8 

RODB 

25 

42.8 

990.6 

Z 

191 

1408.0 

RID 

13 

39.1 

2156.4 

_ 

70U 

66 

33.8 

1090.3 

D 

= 121. 

4 

3ZH 

15 

31.1 

1383.9 

EAT 

41 

32.1 

1690.5 

ST 

= 1044 

.0 

.MA 

106 

609.4 

778.0 

HAA 

11 

351.5 

2087.7 

LAA 

17 

353.7 

1866.5 

Table  24. 

Model  Data:  V, W, X=PR5/Y, Z=PR3 (RW5) 

TOTAL  ARRIVALS:  539 

TOTAL  DEPARTURES:  407 

MA(UP) : 24 

MA  ( P ) : 94 


ROUTE 

# 

D 

ST 

CATEGORY 

# 

ST 

DEP 

407 

54.4 

— 

V 

58 

615.7 

VFR 

28 

44.7 

151.5 

W 

145 

811.8 

OPW 

177 

29.6 

895.1 

X 

60 

873.0 

RODA 

21 

60.1 

1511.9 

Y 

57 

1226.8 

RODB 

21 

40.7 

964.1 

Z 

191 

1424.1 

RID 

7 

43.4 

2156.3 

— . 

70U 

65 

30.5 

1058.8 

D 

= 114. 

4 

3ZH 

12 

25.8 

1275.0 

EAT 

64 

35.3 

1696.3 

ST 

= 1024 

.1 

MA 

118 

543.9 

708.8 

HAA 

7 

363.3 

2124.7 

LAA 

19 

312.3 

1632.0 
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(for  both  random  number  sequences)  whereas  ST  for  Categories 
W and  Y fluctuated  inconsistently.  ST  decreased  in  the 
experimental  model  (both  random  number  sequences)  and  D 
remained  the  same  or  increased.  A pattern  of  change  was  not 
identified.  The  efficiency  factors  for  the  experimental 
model  using  RN1  and  RN5,  were  computed  to  be  103.59  and 
-104.79,  respectively. 

The  third,  and  last,  experiment  involving  service 
priority  by  aircraft  speed  category  examined  the  effect  of 
granting  priority  to  only  Category  V aircraft.  Category  V 
was  assigned  PR5  and  PR3  was  assigned  to  the  remaining  cate- 
gories. The  hypothesis  being  tested  was  similar  to  those 
that  preceded:  ST  should  decrease  and  the  change  should 

be  reflected  by  a corresponding  decrease  in  ST  for  Category 
V and  an  increase  in  ST  for  the  remaining  categories. 

This  experiment  was  also  conducted  using  two  dif- 
ferent random  number  sequences  (RNl,  RN5)  and  the  resultant 
computer  outputs  appear  in  Tables  25  and  26.  The  data  tables 
were  compared  to  those  of  the  basic  model.  ST  for  Cate- 
gories V and  2 decreased  in  the  experimental  model  while 
sT  increased  for  Categories  W and  Y.  ST  for  Category  X did 
not  experience  a consistent  change.  D was  also  inconsistent: 
it  increased  for  RNl  and  decreased  for  RN5 . ST  decreased 
in  the  experimental  model  for  both  random  number  sequences. 

A pattern  to  characterize  change  in  system  parameters  was 
not  detected.  An  efficiency  factor  of  45.96  was  computed 
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Table  25 


Model  Data: 

V=PR5/W, 

X, Y, Z=PR3 (RNl) 

TOTAL 

ARRIVALS : 

508 

TOTAL 

DEPARTURES 

: 432 

MA ( UP ) 

: 19 

MA  ( P ) : 

72 

ROUTE 

# 

D 

ST 

CATEGORY 

# 

ST 

DEP 

432 

64.3 

— 

V 

46 

637.0 

VFR 

42 

50.6 

155.1 

W 

134 

869.8 

OPW 

164 

28.8 

893.2 

X 

56 

921.9 

RODA 

11 

60.5 

1501.5 

Y 

53 

1127.0 

RODB 

24 

40.0 

1012.2 

Z 

177 

1418.8 

RID 

8 

38.0 

2095.1 

70U 

65 

32.8 

1054.4 

D 

= 106.4 

3ZH 

12 

37.3 

1444.3 

EAT 

57 

36.8 

1713.8 

ST 

= 1013. 

5 

MA 

91 

527.5 

698.6 

HAA 

10 

340.5 

1913.1 

LAA 

24 

305.6 

1583.8 

Table  26. 


Model  Data:  V=PR5/W, X,Y, Z=PR3 ( RN5) 

TOTAL  ARRIVALS:  477 

TOTAL  DEPARTURES:  451 

MA ( UP ) : 18 

MA (P) : 66 


ROUTE 

# 

D 

ST 

CATEGORY 

# 

ST 

DEP 

451 

67.5 

— 

V 

41 

669.9 

VFR 

29 

71.4 

176.6 

W 

133 

862.9 

OPW 

150 

28.6 

893.6 

X 

55 

987.7 

RODA 

17 

62.6 

1564.0 

Y 

53 

1199.8 

RODB 

19 

42.1 

973.4 

Z 

166 

1380.5 

RID 

7 

37.1 

2001.0 

__ 

70U 

76 

35.8 

1071.5 

D 

= 108. 

6 

3ZH 

13 

28.7 

1337.3 



EAT 

55 

36.7 

1723.6 

ST 

= 1036 

.5 

MA 

84 

575.0 

742.3 

HAA 

9 

319.0 

1474.2 

LAA 

18 

308.9 

1496.8 
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for  the  experimental  model  using  RN1 . The  efficiency  factor 
for  RN5  is  -14.69. 

Experiment  Involving  Arrival  Route  Re-design 

One  experiment  of  this  type  was  conducted.  It  con- 
sisted of  several  alterations  to  the  basic  program  to 
structure  the  arrival  routes  and  the  final  approach  course 
differently.  A complete  listing  of  the  altered  program  is 
contained  in  Appendix  B.  In  this  revised  system,  all  arriv- 
ing aircraft  are  routed  to  an  intermediate  point  twelve 
nautical  miles  Northeast  of  the  FAF.  From  this  point  the 
aircraft  fly  a straight-in  course  to  the  FAF  and  continue 
to  the  runway.  The  revised  system  configuration  is  shown 
in  Figure  27. 

The  airspace  between  the  FAF  and  the  intermediate 
point  allows  the  controller  to  provide  altitude  separation 
for  aircraft  proceeding  to  the  FAF.  The  fastest  category  of 
aircraft  being  considered  in  the  model  can  descend  4000  feet 
between  the  intermediate  point  and  the  FAF  using  a descent 
rate  of  1000  feet  per  minute.  Therefore,  four  levels  of 
altitude  are  provided  in  the  model  for  aircraft  arriving  at 
the  intermediate  point.  Aircraft  at  the  highest  level  (6000 
feet  mean  sea  level)  must  begin  descent  immediately  after 
departing  the  intermediate  point.  Aircraft  at  the  second 
level  (5000)  can  maintain  that  level  for  up  to  three  nautical 
miles  before  beginning  descent.  Aircraft  at  each  lower  level 
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Figure  27.  Revised  Arrival  Routes 


can  maintain  their  altitude  progressively  longer.  Aircraft 
at  the  fifth,  and  lowest  level,  do  not  descent  until  reach- 
ing the  FAF.  It  should  be  recognized  that  when  an  aircraft 
departs  a particular  altitude,  the  altitude  below  must  be 
clear  of  other  aircraft.  "Facilities"  were  used  in  the  model 
to  replicate  this  situation.  An  aircraft  descending  to  the 
next  lower  level  will  "sieze"  the  "facility"  (representing 
the  airspace  it  is  entering)  thus  precluding  entry  of  other 
aircraft.  The  arrangement  of  facilities  used  in  the  model 
is  conceptualized  in  Figure  28.  The  remainder  of  the  model 
was  retained  in  its  original  form. 

The  altered  model  was  run  on  the  computer  using 
RNl  and  RN5 . Tables  27  and  28  present  a summary  of  the 
computer  output.  Comparison  of  ST  with  the  basic  model 
reveals  a decrease  using  both  RNl  and  RN5 . D decreased 
using  RN5  and  increased  slightly  using  RNl.  A graphic 
comparison  of  ST  by  route,  D by  route,  and  ST  by  category 
is  presented  in  Figures  29  through  34.  In  each  figure, 
the  solid  line  is  the  basic  model  and  the  dashed  line  is 
the  experimental  model.  As  with  previous  experiments,  no 
consistency  in  direction  or  magnitude  of  change  was  identi- 
fied, nor  was  a pattern  of  change  detected.  The  value  of 
the  efficiency  factor  for  this  model  is  +27.14  when  RNl 
is  used,  and  -103.56  when  RN5  is  used. 
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Table  27 


Model  Data:  Straight  In  (RN1 ) 

TOTAL  ARRIVALS:  502 

TOTAL  DEPARTURES:  433 

MA(UP) : 6 

MA ( P ) : 78 


ROUTE 

# 

D 

ST 

CATEGORY 

# ST 

DEP 

433 

44.7 

— 

V 

24  710.7 

VFR 

46 

91.4 

197.5 

W 

115  888.3 

OPW 

157 

28.7 

925.3 

X 

70  937.8 

RODA 

- 

- 

- 

Y 

72  1200.3 

RODB 

36 

36.3 

1048.5 

X 

175  1486.0 

RID 

10 

31.5 

1628.3 



70U 

63 

35.7 

1210.9 

D 

= 99.1 

3ZH 

14 

22.9 

1420.7 



EAT 

65 

36.8 

1905.5 

ST 

= 1076.5 

MA 

84 

587.7 

763.4 

HAA 

11 

307.6 

1977.6 

LAA 

17 

309.5 

1581.1 

Table  28. 

Model  Data:  Straight  In  (RN5) 

TOTAL  ARRIVALS:  508 

TOTAL  DEPARTURES:  412 

MA(UP) : 18 

MA ( P ) : 76 


ROUTE 

# 

D 

St 

CATEGORY 

# 

ST 

DEP 

412 

56.0 

_ 

V 

32 

699.7 

VFR 

33 

58.6 

163.6 

W 

132 

849.6 

OPW 

176 

30.2 

948.1 

X 

62 

922.7 

RODA 

- 

- 

- 

Y 

64 

1241.9 

RODB 

42 

40.4 

1041.4 

Z 

185 

1407.9 

RID 

8 

46.6 

2001.1 

__ 

70U 

58 

35.6 

1201.1 

D 

= 105. 

8 

3ZH 

15 

29.4 

1505.7 

___ 

EAT 

52 

32.6 

1816.8 

ST 

= 1057 

. 3 

MA 

94 

538.4 

716.6 

HAA 

6 

357.0 

1924.8 

LAA 

25 

321.0 

1606.1 

BASIC  MODEL:  SOLID  LINE 

EXPERIMENTAL  MODEL:  DASHED  LINE 
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Figure  29.  Average  System  Time  by  Route  - Experimental  Mode 
(Straight  In  - RNl)  versus  Basic  Model 


BASIC  MODEL:  SOLID  LINE 

EXPERIMENTAL  MODEL:  DASHED  LINE 
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Figure  30.  Average  System  Time  by  Route  - Experimental  Model 
(Straight  In  - RN5 ) versus  Basic  Model 


BASIC  MODEL:  SOLID  LINE 

EXPERIMENTAL  MODEL:  DASHED  LINE 


Figure  31.  Average  Delay  by  Route  - Experimental  Model 
(Straight  In  - RN1 ) versus  Basic  Model 


BASIC  MODEL:  SOLID  LINE 

530  | EXPERIMENTAL  MODEL:  DASHED  LINE 


Figure  32.  Average  Delay  by  Route  - Experimental  Model 
(Straight  In  - RN5 ) versus  Basic  Model 


BASIC  MODEL:  SOLID  LINE 

jUJ  T EXPERIMENTAL  MODEL:  DASHED  LINE 


! 


i 


132 


Figure  33.  Average  System  Time  by  Category  - Experimental  Model 
(Straight  In  - RNl ) versus  Basic  Model 
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Summary  of  Experiment  Results 


Analysis  of  the  computer  output  of  the  experiments 
described  above  failed  to  reveal  a pattern  in  the  changes 
produced  by  the  experimental  model.  Comparison  of  the  effi- 
ciency factors  computed  for  each  model  shows  that  models 
using  RNl  tended  to  have  a positive  efficiency  factor,  while 
models  using  RN5  tend  to  have  a negative  efficiency  factor. 
Figure  35  shows  these  tendencies  clearly.  The  solid  line 
connects  efficiency  factors  from  models  using  RNl,  and  the 
dashed  line  connects  efficiency  factors  from  models  using 
RN5.  This  dichotomy  is  apparently  the  result  of  variation 
in  the  behavior  of  the  basic  model  caused  by  changing  the 
random  number  sequence.  This  variation  in  the  basic  model 
was  compensated  for  by  "standardizing"  the  efficiency 
factors;  that  is,  the  mean  of  the  efficiency  factors  from 
models  using  RNl  was  subtracted  from  each  efficiency  factor 
in  that  group.  A similar  correction  was  applied  to  each 
efficiency  factor  from  models  using  RN5.  The  corrected 
efficiency  factors  are  shown  in  Figure  36  to  facilitate 
comparison.  Model  number  three  (OPW,  70U,  EAT  = PR5) 
has  a positive  efficiency  factor  using  both  RNl  and  RN5. 
Model  number  seven  (Straight-In)  has  a negative  efficiency 
factor  in  both  cases.  All  other  models  have  one  positive 
and  one  negative  efficiency  factor.  In  summary,  analysis 
of  data  produced  by  experimentation  failed  to  identify  any 
significant  changes  in  system  performance. 
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Figure  36.  Standardized  Efficiency  Factor  by  Model  - RNl  versus  RN5 


Chapter  5 


CONCLUSIONS  AND  RECOMMENDATIONS 

Four  research  questions  were  posed  in  Chapter  1. 
Chapter  5 addresses  these  questions  under  two  major  subareas: 
(1)  Conclusions,  and  (2)  Recommendations  for  Future  Research. 

Conclusions 

The  first  research  question  was  : Can  the  variables 

comprising  the  military  ATC  environment  be  adequately  repre- 
sented by  a computer  model?  The  development  process  of  the 
Dayton  ACA  simulation  model  and  its  subsequent  validation 
indicate  that  the  military  ATC  environment  lends  itself  well 
to  computer  simulation.  It  has  been  concluded  that  the  Dayton 
ACA  model  is  a satisfactory  representation  of  the  true  envi- 
ronment. This  conclusion  was  considered  a prerequisite  for 
subsequent  use  of  the  model  in  experimentation. 

The  second  research  question  was:  Can  changes  in  the 

terminal  environment  be  effectively  evaluated  by  the  simula- 
tion model?  Several  experiments  were  conducted  with  the 
model  to  test  various  hypotheses  related  to  the  efficiency 
of  the  total  system.  The  data  resulting  from  the  experiments 
fluctuated  with  apparent  randomness.  Possible  causes  for  the 
fluctuations  include:  (1)  low  traffic  volume  at  Wright- 

Patterson;  (2)  inadequate  correction  for  variation  in  speed 
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category  mix  within  the  simulation  model;  (3)  insufficient 
length  in  the  simulation  run,  and  (4)  treatment  of  missed 
approaches  in  the  simulation.  Each  of  these  possible  causes 
is  discussed  below. 

The  volume  of  air  traffic  using  Wright-Patterson  is 
well  below  system  capacity.  This  low  traffic  volume  could 
account  for  the  random  fluctuation  in  the  data  resulting 
from  the  experiments.  The  occurrence  of  greater  or  lesser 
amounts  of  aircraft  delay  could  be  attributable  to  chance 
conflict  of  aircraft  rather  than'  the  experimental  changes 
in  the  system. 

The  random  fluctuation  in  output  data  suggests 
the  length  of  the  simulation  runs  may  have  been  insufficient. 
This  would  mean  that  the  simulation  had  not  reached  steady 
state  conditions  with  respect  to  the  characteristics  of 
system  aircraft.  Results  from  the  experimental  models 
exhibited  wide  disparity  in  the  number  of  aircraft,  their 
categories,  and  their  routes.  The  disparity  can  be  accounted 
for,  in  part,  by  the  effect  of  the  PRIORITY  block.  As 
explained  earlier,  this  statement  causes  GPSS  to  re-initiate 
scan  of  the  current  events  chain  resulting  in  re-distribution 
of  aircraft.  The  effect  of  the  re-distribution  could  be 
dampened  by  lengthening  the  computer  simulation  run.  Alter- 
natively, the  effect  of  the  PRIORITY  block  might  be 
eliminated  by  more  sophisticated  programming  techniques. 
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A consequence  of  the  re-distribution  could  be  the 
generation  of  a higher  or  lower  percentage  of  aircraft,  in 


the  experimental  model,  of  a particular  category  than  in 
the  basic  model.  Similarly,  there  could  be  a higher  or 
lower  percentage  on  a particular  SVR.  The  impact  of  these 
changes  on  a given  experiment  is  not  known  with  certainty. 

The  correction  factor  (a^)  applied  in  the  computation  of 
the  efficiency  factor  may  not  have  adequately  compensated 
for  the  changes.  The  weight  of  the  correction  factor  can 
be  changed  by  changing  the  value  of  k (100). 

Another  possible  cause  of  the  fluctuation  in  the 
data  resulting  from  experimentation  is  the  treatment  of 
missed  approaches  (MAs)  in  the  model.  In  the  majority  of 
experiments,  the  changes  in  D for  MA  aircraft  tended  to 
dominate  the  statistics  for  the  total  system,  possibly 
obscuring  some  change  pattern  otherwise  identifiable.  In 
the  analysis  of  the  computer  output  of  the  various  models 
(described  earlier)  the  entire  time  required  by  an  aircraft 
to  execute  an  MA  was  treated  as  delay  time. 

No  conclusions  could  be  drawn  from  the  experiments 
conducted  with  the  simulation  model  because  of  the  randomness 
exhibited  by  the  data.  The  possible  causes  for  the  random 
fluctuation  in  the  data  discussed  above  indicate  that  the 
model  is  not  necessarily  ineffective  in  evaluation  of 
experiments.  Therefore,  the  research  question  remains 
unanswered . 
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The  last  two  research  questions  dealt  with  prediction 
of  conflict  and  conflict  resolution.  These  areas  were  not 
addressed  by  this  research  due  to  the  time  constraints 
imposed  by  the  academic  situation. 

Recommendations  for  Future  Research 

Further  research  is  recommended  in  the  following 
areas.  First,  addition  experiments  should  be  conducted 
with  the  model.  This  research  was  limited  to  experiments 
involving  priority  changes  and  one  attempt  at  system  redesign. 
There  is  a wide  variety  of  design  configurations  and  traffic 
control  decision  rules  that  could  be  tested  other  than  those 
chosen. 

If  further  experiments  are  conducted,  the  following 
recommendations  apply.  The  length  of  the  simulation  run 
should  be  increased  to  insure  that  steady  state  conditions 
are  reached.  In  addition,  the  simulation  model  should  be 
artificially  forced  to  contend  with  a greater  volume  of 
traffic  than  is  characteristic  of  the  real  system.  This 
could  be  accomplished  by  reducing  the  mean  time  between 
aircraft  arrivals.  The  increased  volume  of  traffic  would 
preclude  the  domination  of  the  statistical  output  by  the 
chance  occurrence  of  conflict. 

Another  recommendation  is  that  MAs  be  treated 
differently.  One  possible  approach  is  to  treat  planned  MAs 
as  a separate  traffic  classification  and  record  as  delay 
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only  that  time  in  excess  of  the  normal  time  to  execute  a 
missed  approach  procedure.  Unplanned  MAs  could  be  handled 
separately  and  treated  as  additional  delay.  The  number  of 
unplanned  MAs  could  also  be  used  as  an  independent  criterion 
of  system  efficiency. 

The  next  two  recommendations  deal  with  potential  use 
of  simulation  in  an  on-line,  real-time  context.  Toward  this 
end,  the  simulation  model  could  be  used  to  evaluate  the 
feasibility  of  identifying  conflict  in  a real  traffic  envi- 
ronment. It  is  recommended  that  the  GPSS  program  be 
modified  to  permit  evaluation  of  traffic  conflict  in  the 
model  through  comparative  analysis  of  aircraft  estimated 
arrival  times  at  points  in  the  system.  This  could  be 
accomplished  with  an  algorithm  similar  to  MACRO  eight 
which  allows  the  controller  to  resequence  aircraft  by 
using  airspeed  control.  This  recommendation  is  directed 
toward  answering  research  question  three. 

The  final  recommendation  addresses  research  question 
four  which  concerns  the  use  of  the  simulation  model  as  a 
means  of  evaluating  alternative  solutions  to  conflict 
(2:245-252).  The  feasibility  of  using  on-line,  real-time 
simulation  in  an  ATC  environment  could  be  evaluated  by  the 
model.  Once  conflict  is  identified  in  the  manner  described 
above,  the  algorithm  could  be  expanded  to  compare  various 
estimated  arrival  times  based  on  hypothetical  changes  in 
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airspeed  or  vectoring  route  of  the  aircraft  involved  in  the 
conflict.  If  the  algorithm  finds  a combination  of  changes 
that  resolves  the  conflict,  those  changes  would  be  made  to 
the  appropriate  aircraft.  The  results  obtained  from  an 
experimental  model  using  this  algorithm  could  be  compared 
with  results  from  the  basic  model  to  determine  its  effec- 
tiveness . 

Ideally,  a simulation  model  would  be  on-line  in  an 
ATC  facility  computer,  and  available  to  continuously 
evaluate  existing  traffic  using  an  heuristic  search  pro- 
cedure. The  best  solution  to  a predicted  conflict  situation, 
as  indicated  by  evaluation  of  alternatives  through  simulation, 
would  be  provided  to  the  controller.  Therefore,  it  is 
recommended  that  future  research  explore  the  possibility 
of  using  simulation  in  this  context. 
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APPENDIX  A 


This  appendix  contains  the  simulation  program  used 
to  model  the  Dayton  Approach  Control  Area.  The  program  is 
written  in  GPSS  language  and  was  executed  on  a 635  Honey- 
well Computer  located  at  Headquarters,  AFLC,  Wright- 
Patterson  APB,  Ohio.  The  program  was  stored  in  the  CREATE 
time-sharing  system  using  the  file  name  SIM03,  on  user 
identification  number  77A79  and  password  XS76. 
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26l7*3ASe  LEG  ENTRIES 
262.'* 

26373A5=«ENTER*I 
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36 60FI N AL  J F V AR I ABLE  * I / PL5*36 70*rN  I 4 
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.27  yiASSIGN*  I 4,  VSB.1L  I , PL 

27  I 0 *GATE  MU  * 3DF I , AGO 35 

272 "T .'<0*  MACRO  * ADF  1 

2 7 37  TWO  * MACRO  * 3DF2 

274*1  .•(')*  MACRO  *3DF3 

2757  *TEST  E*PL4, 227,ADD97 

2760* 

2 770+HSAV IES  FLY  FINAL 
2730* 

2790/HPEE  * MACRO *BDF4,^GF I 

2600THREE  »MACR(J*BDr5,3Dr  2 

251 0AJJ2S * PREEMPT* iQF6, PR 

2r 20  * ADVANCE  * PL  I 4 

2 " 20  JP  ET'JPN  J 3DF  3 

254  2 * ENTER  * I 0 

2657F0 JR  J MACRO «3DF7,3DF4 

26  60FOUR J MACRO  * 5 DF3 , 8DF5 

26  77FO JR  ; MACRO  *SDF9, 3DF6 

3660 AG 029 1 ENTER  * I I 

2590FOUR  * MACRO *3DF 1 2, 5DF 7 

29 00 FOUR  J MACRO  »3DF  I I ,^F5 

29  I 7F0UR  * MACRO  JriDF  1 2 , 3DF9 

2920  J PREEMPT t BDF I 3, PR 

2930  J ADVANCE  *PWI 5 
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2957  *T'EST  NE  5 PL  I 8.9,  ADD52 
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3 1 4 0 » T£  ST  LE  t P L 2 5 , ?l  26 , A JD 9 1 
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3170ADO9I  * TE  5 T LE ; PL  25  , PL2  7,  AEJ92 
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373  '?tfa:;sff:2 ; , adh 33 
37? :+ 
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APPENDIX  B 


This  appendix  contains  the  simulation  program  used 
to  model  the  "straight  in"  experimental  design  of  the 
Dayton  Approach  Control  Area.  The  program  is  written  in 
GPSS  language  and  was  executed  on  a 635  Honeywell  Computer 
located  at  Headquarters,  AFLC , Wr ight-Patterson  AFB , Ohio. 
The  program  was  stored  in  the  CREATE  time-sharing  system 
using  the  file  name  SIM02,  on  user  identification  number 
77A79  and  password  XS7o. 
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001 0#*S,R(SL)  », 8, 16?;, 8, 19, 31 

0(320$  ? I DENT  ? NP  1191,  AFIT/SLG , LORENZ  AMDGIPEAR  77A 

0030$  ? PROGRA" : PLFS,  ON  I 

004  0$  i LIMITS?  5, 39L,,  5K 

ww5 OS ? PR"FL ? ! * , R , R , AF  . L I a /GPSSH S 

0060$ * FILE? * I , A I 3,3L 

0070$ ?FILE?*2,B IR,2L 

0080$* FILE? *3, Cl R, 1 L 

0090» SI»ULATE; , ,2000 

0 1 00 ? CONTROL  ? FLO , 700, VA  P , 80  , CUE , 60 , FAC, 50 , F"S , 1 5 
0 1 1 0? UNLIST 


0 1 20* 

0 I 30*DEFI N I TION  OF  MATRICES 
0140* 


015  0 

2 1 MATRIX 

My  s 1 

0 160 

3? MATRIX 

,rn,5,  I 

0 1 70 

4 ; "AT0 IX 

l‘F,5,  1 

0180 

5 {MATRIX 

MH,5,  1 

0 190 

6 {MATRIX 

" H , 5 , 1 

0200 

7 {MATRIX 

*■'  H , 5 , 1 

0210 

8 {MATRIX 

,fH,  5 , 1 

0220 

o 5 MATRIX 

MM  1 

• * J t 1 

0230  I 0 ? MATRI  X ; "F , 12,3 

0240  I 1 r.'ATRI  X? "F, 1 5,  I 

025  0 12  ?UATRI X ? MF , 15,1 

02  63  I 3?"ATRf  X?’<F,  I 5,  I 

0273?  INITIALT-'F  ' ' ~MFI  3(1-15,1  ) ,0 

04  83? I NIT  I AL?“F 1 0 ( I - 1 2 , 1 -8 ) , 0 


3290?  I N IT  I AL?  MH2-MH9  ( I -5,  I - 11  ) ,3 
0300* 


03  I 3*FUNCTI0NS 
0320* 


0330  1 ? FUNCTION  ?RNI  ,C24  {EXPONENTIAL  C I STR I P UTION 

034  0*0,0/, 1 , . I04/.2 , .222/. 3, .3 55/. 4, .509/. 5, .69/. 6,  .9  15 

0350*. 7, I .2/. 75, 1 .38/. 8, 1 . 6/ . 84 , I . 83/ . 88,2 . 1 2/ . 9 , 2 . 3/ .92 

0360#. 94 ,2. 81/. 95, 2 .99/. 06, 3. 2/. 9 7, 3. 5/. 98, 3. 9/. 99, 4.6 

0373#. 995, 5. 3/. 998, 6. 2/. 999, 7/. 9998, 8 

0380  2 5 FUNCTI ON  ? RN 1 , C25  ? STANDARD  NORMAL  DISTRIBUTION 

0j90#0 , -5/ . 00003 , -4/ . 00  I 35 , -3/ . 0062  I , -2. 5/. 02275,  -3 

04  00*. 0 6681 ,-! .5/. 1 1507,-1 .2/. 1 5866,-1/. 2 I 1 86, -.8/. 27425 

0410* .34458, -.4/. 42074, -.2/. 5,0/. 5792*, . 2/. 6 554 2, .4 

0420*. 72575, .6/. 788 14, .8/. 84  I 34, 1/, 88493, 1 . 2/ .9 33 1 9 , 1 .5 

0430* .97725 ,2/. 993 79, 2. 5/. 99865, 3/. 99097,4/1 ,5 

0440  3; FUNCTION ? RN 1 ,D2? RULES 

0450#. 7269, ADDI /.727,ATD2 

0460  4 ; FUNCTION f RN 1 ,D5? CATEGORY 

047 0# . 062 » A DD3/ . 3 1 0 , A DD4/ ,4  34  f A DD5/ . 562 , A DD */.562  1 , AC  r'7 

0480  5 ? FUNCTION ? 9N I ,D2 lOPERATI ON  1 

0490#.  5,  ADDS/. 5 001  ,AI)D9 

0500  6 ? FUNCTION IPN I ,D3 ? OPERATION  2 

051 0# .536, ADDS/ .965 , ADD9/.965 1 , \ CD  I 0 

0520  7 ? FUNCTI ON iRNI ,D2 lOPERATION  3 


52 


6 
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053v5*  .393 , A CDS/.  3 93  I , ADD9 
0540  B ; FUNCTION  * RM I ,D3 ; OPERATION  4 
0 550# . 4 8 3 , A CD  8 / . 9 9 7 , ADD  9/ . 39  7 1 , A DD 10 
0560  9 l FUNCTION ;RNI , D3 ? OPERATION  5 
0 5 7 0#  . 3 74  , A CC  6 / . 8^8  , A DC  9 / . 5 08  1 , A DD  1 •? 

0580  1 0; FUNCTION ; RN ! , D3 ; AI RPORT 
0590#  . I 68,  ADD  I I / . 4 70,  ADD  I 2/.  4 7")  I , ADD  I 3 
0600  1 1 ; FUNCTION ; RN I , D5  ; ROUTE  1 

.'6:  0#. 056, ADD  14/.  1 67, ADD  15/. 724, ADD  1 6/. 891 , ADD  17/ .891 
0620  12  ; FUNCTION  « RN 1 ,D6;ROUTE2 

0630# . 48 1 , ADD  I 5/. 592, ATD I 6/ .629 , ADD  I 7/ .8 1 4 , ADD  I 8 
0640#.  351  , A CD  19/.  351  1 , A DD2r/ 

0650  1 3;FUN JTI  )N ; RN  I , D4 ; ROUTE3 

0660#. 222, ADD  14/. 500, ATD 1 6/.870,ADD1 7/. 8 70  I ,ADD20 
0670  1 4 {FUNCTION*  FM2,C3  ? ADVANCE 
0680#-5, .90/0, 1 .00/5, 1.10 
069  0 1 5 ; FUNCTION i FN2 , C3  ; I A 
0700#-5, .80/0, 1 . 00/5, 1 .20 
0/10* 

0 /20*MACRO  STATEMENTS 
0730* 

0740ONE iSTARTNACRO 

075  01  “ARK  5 1 PF 

0760 » ASSIGN ; 23, #A ,PL 

v'770|ASSIGN  ;24,#P,PL 

0780; TRANSFER » . 1 00, , ADD 70 

0790; ENDMACRO 

0800* 

031  OT'AO  ;STAF.TMACRO 
082O; PREEMPT  s #A , PR 
0830; ADVANCE; PL  14 

084  3;  ENDMACRO 

085  0* 

0 860T HR  EE  J STA  PT''ACR() 

0870; RPEEMPT; #A, PR 

0 880; ADVANCE; PL1  4 
0890; RETURN;#? 

0900; EMD"ACRO 
091  3* 

092OFOUR ; START’*  ACRO 
0930; PREEMPT ; #A , PR 
0940; ADVANCE; PH  15 
0950; RETURN ;#B 
0960; ENDMACRO 
0970* 

098 0F I VE;START"ACRO 
0V90;ASSIGN;#A,#B,PL 
1 000; ASSIGN ;#C,  #D,PL 
l 01 0; ASSIGN ; #E , #F , PL 

1 023;  ENDMACRO 
I 03  0* 

1 040S I X ; STAPT’'ACR() 


160 


'5-'  VAR#A;F  VARI APLE;  PL*B-dL25 
'6' ' I SA  vT£  VA LUE  * *C , P L-#  , XL 
07  ' ; A SSI  GN  ? 3 , VS  V AR?M  , PH 
'30;  ADVANCE?  V$  VAR-  A 
09  TRANSFER  ? , ACD3G 
I ‘-'0 ; END"ACR0 
] I ■* 

I20SE  VEN?START"ACR0 
1 3r,APT"A  ?NULL 

1  40RTE4P  IF  VARI  APLE?  ! 0f'/3L3*3600 
1 5 GS  VR#B  » F VARI APLE  » *C/PL4*3AO0 
160?  ADVANCE ? VS  RTEA3, FN I 4 

1 700';e?"ACR()  ;«*d,  *E 

1 8 0E I GHT  ? "AC R0  ? *3 , #G , F 
1 9 0;  END'1  ACRO 
200* 

21 TEIGfT ; START’^ACRO 

2 20ETA #A ? VARI ABLE  ? AC  1 +VS SVR^A 

23  'A ETA"  A ? VARI  APLE  ? V$ETA#A+  VSVAR6 
240CETA#A*VARIABLE ? VSETA#A+VS VAR7-MF*B (PH6  * 1 ) 
250DETA#A;  VARI  ABLE  tVSA ETA  *A+V$VAR7-’<F#3  ( PF6,  I ) 
260ADR  I <*A?TEST  G;"F*B(  PF6 , I ) , 0 , ADR4*A 
2 70 ; TEST  G ; VSCETA#A,  0,  APR2AA 

2 8 0;  TEST  G ? VS  D ETA  P A , 0 , AD  R3  - A 

29  • 'A  T R2  «A? SAVE  VALUE; 10, V$  ETA  - A , XF 

300ADR4*A ; LOOP  ;6P!',  ADR  1 * A 

3!0;TRANSFERi,ADRe#A 

320ADR3#A?SAVE  VALUE;  ! 0,  VSAETA#A,XF 

330ADR8#  A ; SAVE  VALUE  ? 11  + , ! , XF 

340 ; TEST  GE  * XH 1 1 ,16, ADR7“A 

350 ; SAVEVALUE ? 11,1 ,XH 

3<V3ADR7*Al  TEST  E ? XF  1 0 , VS  A ETA  A' A , ADR5#A 

37  0?  MS A VE VALUE  ;*B, XH 1 I , 1 , VSAETA*A, "F 

3 = 0ASVR#A? VARIABLE  ? VSSVRsLA-VS VAR6 
390 ; AD VA NICE  ? VS  AS VR #A  , FI ’ 1 4 

400 ; TRANSFER ? , ADR6-A 

4 1 0APR5#  A ? MSA  V'E  VALUE  ? -B , XH  1 1 , 1 , VS  ETA -A  , '-'B 
42?; ADVANCE; VSSVR^A, FT  1 4 

430APR6*A ; TRANSFER ? , *C 
440 ; ENDMACRO 
45  0* 

460*START 

470* 

4 8 0;  GENERATE 1 82, FN I , ,5  000,  3, 40PL.20PH,  I5PF 

490 ; A SSI GL  ? 6 , 5 , PM 

5 00;  ASSIGN;  I 7,0,  PH 

5 1G» ASSIGN; 18, 12.PL 

520;  ASSIGNS,  12,  PL 

530;ASSIGN;6, 1 2, PL 

540; ASSIGN ;7, 12, PL 

550* 

560*FLIGET  RULES 
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157.?* 

I 58Q ? TRANSFER  ;FN, 3, 0 
1 590ACD1  ; ASSIGN  ; 1 , 1 , PL 
I 600;TRANSFER;FN,4,0 
1 6 1 0 A DD2  ; A SS I GN  ; I ,2, PL 
I 620;TRANSFER ; , ADD  I 02 
i 630* 

! 64 '’'★CATEGORY 
l 65  0* 

I 660ADD3;ASSIGN;2,90,PH 
! 670; ASSIGN’; 2,  ! , PL 
1 66  0FI VE; VACRO ; 3, 550,4,250,5, I 70 
1 69n  ; TRANSFER ; FN , 5 , 0 
I 7 00A  CD4 ; A SS  I GN ; 2 , I 20 , P F 
I 71 0; ASSIGN ;2, 2, PL 
I 720FI VE; wACRO; 3,500,4, 220,5 , 1 60 
1 730; TRANSFER IFN, 6,0 
I 74 0A  DD5 ; A S S I GN ; 2 , 6 0 , P H 
I 75 0*  ASSIGN ; 2 , 3 , PL 
1 760FI VE; MACRO ;3, 300,4, 170,5,1 50 
1 770; TRANSFER ;FN, 7,0 
; 78  1ADI  6;ASSICN?2,50,PH 
1 790; ASSIGN; 2, 4, PL 
1 8 ’ VI  VE;MA  :R(  ; 3, 250,  4, 130, 5,  1 I 0 
i 81  '; TRANSFER ;FN , 8,0 
! *2  0A  PD  7 ; A SS  I GN ; 2 , 20 , PF 
1330;  \ SS  I G‘i ; 2 , 5 , PL 
1 840FI VE;MACR(  ;3, 140,4,  110,5,80 
1 350 ; TRANSFER ;FN, 9, 0 
1 860* 

1 370*OPERATION 
1 880* 

! 390ADD8; ASSIGN ;6, 1 ,DL 
19.  ' ; TRANSFER » FN, 1 1 ,0 
I 9 1 0A  PD9 ; ASSI GN ; 6 , 2 , PL 
1 92!'  ITRAMSFER  ;FN,  10,0 
1 9 3 GA  CL  1 0 ; A SS I GN ; 6 , 3 , PL 
1 940 1TRANSFER; , ADD  1 00 
I950* 

1 960*A  I RPORT 
1970* 

1 98 GATT;  I 1 ; ASSIGN;  7,  1 , PL 
1 990; TRANSFER; .ADD  1 00 

2 : '00  A DD 1 2 ; A SS  I GN ; 7 , 2 , PL 

I 01 T!  i IPL6, 1 ,A1 
2020; TRANSFER ;FN, 12,0 
2 ’30ADD 1 3;  ASSIGN; 7,  3,  PL 
i . 1 4 ;T:-.'  I : , ' 

205  '* 

2"60*ROUTE 
2 ’ 7 

■-  )A  TP  1 5 ; A SS  I GN » 3 , |o,PL 
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2 094F  I VE  { MACRO  ;9 , 2 , Hi,  2,  IS,  I 
2 I .'0{  TRANSFER  ; , ADD2  I 
2 ! I 0ADD 16 {TRANSFER; f ROCF 

2 1 23R0DAI ASSIGN  18, *2, PL 

2 t 3 OF  I VEl’<ACR0;9,  I 2 , 1 0,  1 2 , 1 6 , 2 

2 I 4 J; TRANSFER? , ADD 22 

2 15  3RODB  5 ASSIGN?  3,25.5,  PL 

2 160FIVE{'‘ACRO{9,7.5,  10,16.5,18,3 

2 170; TRANSFER ;,ADD23 

2 i SNA DD I 7; ASSIGN; 3, 57, PL 

2 I 90FI VE; MACRO; 9, 4,  I 5,7.5, I 6,4 

220-5; TRANSFER  ; , ADD24 

22 1 0 ADD  1 6 {ASSIGN ; 8, 24 .5 , PL 

2 220FI VE; MACRO ;9 ,4 , 10,4,18,5 

2*23 4; TRANSFER  ? , ACD25 

2240ADD 1 9 ; ASSIGN ; 3, 35 , PL 

22 50F I VE ; MACRO ; 9 , 2 .5, 10,6, 18,6 

2260;TRANSFER  ? , ADD26 

22 70 A CD2 0 {ASSIGN ; 8, 52 , PL 

2280FI  VE ; -MACRO;  9,4,  10,7.5,18,7 

2297 {TRANSFER ; , ADD2 7 

2300* 

231 3*0 PW 
2320* 

2 334VAR6; FVARI ABLE;  ( ( 1 0/PL4  ) - ( 1 0/PL5))*36C0 
2 340VAR7;  VARI  ABLE;  PL1  4*3 
2 35 0St VE N ; 1 ' AC RO ; 2 1 ,2,  I 5 ,36, 24, BASE,  13 
2 860* 

23/o*ROD(A) 

2 38  0* 

2 390SEVEN ; *'ACRO  ;2  2 , 3, 42 ,29 , 37 , BASE , 1 3 
2400* 

24  I 0*ROD( 3 ) 

2420* 

24 3 0SE VEN {MACRO ; 23 , 4 , 20 .5 , 29 , 37 , 5 ASE, I 3 

2 44  M* 

2453*RID 

2460* 

24  70SE VEN {MACRO {24,5,50.5,28, 44, BASE,  1 3 
2480* 

2 4 90* 70 U 
2500* 

25 1 3SEVEN ; MACRO {25,6,24.5,34, 1 5, BASE,  1 3 
2524* 

2530*3ZH 

2544* 

255 0Sc  VEN {MACRO { 26 , 7, 37 .5,28, 17, BASE,  13 
2d  60* 

2574*EAT 

2580* 

2594SEVEN I NACRO {27 ,6 , 52 ,27, 35 ,5ASE, 13 
2604* 
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26 1 3*STRAIGhT  IN  ENTRIES 
2620* 

2630BASE; ENTER; I 

2 640; ASSIGN; 16, I 6, PL 

265 H; ASSIGN ; 1 9, I , PL 

2 660; ASSIGN; 15, AC  I , PF 

2670BDL1  ;F'/ARI  ABLE;  1 /PL4*360O*FN I 4 

2680FI NAL iFVARI ABLE ; I /PL5*3600*FN I 4 

2 690 ; ASSI 3N; 1 5, V$ FINAL, PH 

2700; ASSIGN; 14, VSBDLI ipL 

2 7 1 0 ; G ATE  MU;BDF1 ,ADD85 

2 720TWO  ;UACR() ; 3DF 1 

2 730TNO ;"ACR0 ;BDF2 

2740TWO ; MACRO ; 3DF3 

2 750; TEST  E ; PL4 ,220 , ADC90 

2 760* 

2770*HEAVIES  FLY  FINAL 
2 780* 

2790THREE  ? MACRO ; 3DF4 , BDF 1 
2 8 00THR  EE ; m ACRO ; BDF  4 A , B DF2 
2 3 I 0THREE ; "ACRO ; BDF4B , FDF3 
2 82  0THR EE ; MACRO ; RDF4C , BCF4 
2 83 0THR EE ; MACRO; BDF 5, BDF4 A 
2840ADD28 ; PREEMPT ;BDF6, PR 
2850; ADVANCE; PL  14 
2 360 ; RETURN; 3PF4B 
2 370; ENTER; 10 
2 88 OF OUR ; MACRO ; BDF  7 , BDF4C 
2 390FOUR; uACRO ; BDF8 , BDF5 
2900FOUR; "ACRO ; BDF9 , BDF6 
291 0ADD29; ENTER; I 1 
2920FOUR ; ''ACRO ; BDF  1 0,  BDF7 
2v30F0UR  ;'1ACR();  BDF  1 1 , BDFS 
294  0FOUP; MACRO ; BDF I 2, BDF9 
2950; PREEMPT ;BDF 1 3, PR 

2 960; ADVANCE; PHI  5 

297 3A CD 33  ;TEST  ME ; PL1  8 , 8 , ADD30 
2980; TEST  NE; PL  1 8 ,9 , ADC62 
2990; TEST  NEIPLI8, 10.ADD80 
3 0001 TEST  ME; PL  1 8 , 11, ADD80 
301 0; SPLIT; 5, ADD95, 1 7PH 
302  4* 

3030*CO*'PUTE  DELAY  TP'ES 
3040* 

3050 ; SA  VE  VALUE ;2-6 , 0,  XL 

3060TI M I ; F VAR  I ABLE ; AC  1 -FF 1 5-( PL  1 4*PL 1 6 ) + ( ( PHI 5-PL 1 4 ) *7 ) 

JPI70I  ASSIGN  125,  VSTP'I  ,PL 

3030DEL 1 : F VAR  I ABLE ; ( ( PL9/PL4 ) *36  00) *FN  1 4 

3090; A SSI GN ; 2 6, VS DELI ,PL 

3 1 00DEL2 ; F VARIABLE;  ( ( PL  1 0/PL4 ) *36  00 ) *FN 1 4 

3 I 1 i;  A SSI  JN  ;27 , VSDEL2  , PL 

3 1 2 UDEL3;  F VAR  I ABLE ;((((!  0/PL5 ) - ( I 0/PL4  ) ) *36  00 ) *FN  1 4 )+\/$  DEL2 
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3 1 305  ASSIGN  ? 28,  VSDEL3,PL 
314,5DEL4;FVARIABLE:  12'7i*FHI4 
3 15  05 ASSIGN 5 29 , V$DEL4,PL 
3 1 6 3; TEST  GE; PL25 , I , A CD 30 
3 1 7057EST  LE5  PL25  , PL26,  ADD9  I 
31S0SIX;”ACRO; 1 ,26,2 
3190* 

32O0ADD91  ;TE57  LE ; PL25 , PL27 , A CD92 
32  I OS  I X ;"ACnO  52  ,27,3 
3 22  3* 

323';'ACD92  ;TEST  LE ; PL25 , PL28 , ACD93 

3244SIX  S>'ACRO?  3,23,4 

J25U* 

3263ACD93  5 TEST  LE;PL25, PL29.ACD94 
32  7051 X J MACRO  ? 4 , 29 , 5 
32S0ACD945  SA VE VALUE 56 , PL25 , XL 
3290* 

3300*RUNUAY 

3310* 

3320ACD30; LEA  VE  5 PL  1 9 
33305  LEAVE; 10 
3340;7RA:!SFER? . 1 00,  ,ADD60 
33505 GATE  MU  5 RVIY,  ADC60 
3360  5 PREEMPT 5 PWY , PR 
3370 ;TES7  E 5 PL4 , 220 , ADD 3 1 
3330  5 RETURN 5 BDF 1 0 
3390ADD3 1 5 RETURN  5 BDF I 1 
3400 ; RETURN  5 3DF  12 
341  0 ; RETURN  5 BDF 1 3 
3 4 2 0 5 A D VAMC  E 5 P H2 , FN  1 4 
34305  RETURN  5 RV.Y 
3440; LEAVE  5 1 1 
34505TRANSFER5.070, ,ADD64 
346  0VAR5  5 VARIABLE 5 PF2-PF1 
34705VAPK52PF 

34305 "SAVE  VALUE  5 10+.PL1  , 1 , 1 ,MF 
34905 "SAVEVALUE  5 1 0+ , PL2 , 2 , I ,UF 
3500  5 MSA  VE VALUE  5 1 0+ , PL6 ,3 , 1 ,».<F 
35  I 05  MSA VEVALUE 5 I0+,PL7,4, | ,"F 
3520  5 "SAVE VALUE  5 1 0+ , PL  1 8 , 5 , I ,MF 
3 530  5 " SA  vE  VALUE  5 1 0+ , PL I 8 , 6 , V$  VAR5 , MF 
35405  MSA  VE  VALUE  5 1 0+ , PL2 , 7,  V$VAR5,MF 
355  0 5 MS AVEVALUE  5 I 0+ , PL2 , 8 , I ,"F 
3 560 5 SAVE  VALUE 5 I 4 + , VS  VAR5,  XF 
3 57 0A CD  1 00 5 TERM  I NATE  5 1 
3530* 

J590*OTHERS  FLY  FINAL 
3600* 

361 0ACD90 5 RETURN  5 BCF1 
3620TKREE  5 "ACRO  5RDF4, 8DF2 
3630TUREE  5 "ACPO ;?DF4A,RCF3 
3640THREE  5 "ACRO  5BDF43 , BDF4 
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365  OTHR  EE ; M ACRO ; BDF4C , BDF4A 
3 660;  PREEMPT ; 3DF5  , PR 
3670; AD  VANCE? PL  I 4 
36S3ADD32  ? RETURN; 3DF4B 
3690THREE;  ”ACRO  ;3DF6, 8PF4C 
3 70  3;  ENTER;  10 
371 0FOUR  ; MACRO ; BDF7 , BCF5 
3720FOUR; MACRO ; BDF8 , BDF6 
3 730; PREEMPT?BDF9, PR 
3 740; ADVANCE; PH  1 5 
3 750ADD34 ; RETURN ; BDF7 
3763; ENTER;  1 1 
3 770FOUP?  MACRO; BDF 1 3, BCFB 
3783F0UR;MACR0; 3DF ! 1 , BDF9 
3793ADD36 ; PREEMPT; BCFI 2 , PR 
3333; ADVANCE; PH  15 
381 0; RETURN ;3DF1 0 
3323; PREEMPT ; BCF 1 3 , PR 
j830; ADVANCE; PH  I 5 
3343;TRAMSFER; ,ADD33 
3850* 

3863*1  NTERCEPT  FI N'AL 
3870* 

3 9801 AP ;ENTER  5 1 3 
3S9  3;  ADVANCE » PH  1 5 
3930 ; P P EE” PT ; E DF8 , P R 
39 1 3;  AD  VANCE ; PHI  5 
3923;PREE”PT;BDF9,PR 
3933;AD VANCE?  PHI  5 

3°40; TEST  E ; PL-1 , 223 , ADD34 
3 y5  0; TRANSFER  ? , At  D29 
3963* 

3973*QUEUES 

3980* 

J993ADD95; ASSIGN; 15, PHI  7, PL 
4 000ADD96; ASSIGN; 1 7, PHI  7, PL 
401 0VARI0;  VARIABLE?  PL  I 7 +PL1 8*5 
4 32 3? QUEUE? VS VAR  10 
403O;TEST  NE;XL*PL1 5,0, APD97 
4340; MSA  VE VALUE  ? PL  1 5+ , PL2 , PL1  3, 1 ,WH 

4 35  3A  DD97;  ASSIGN;  22  , VS  VAR  I 3,  PL 
4060 ; ADVANCE  5 XL*PL 1 5 

4 070;  DEPART;  PL22 
4080  TERMINATE 
4390* 

4 I 33*DE PARTURES 
4 113* 

4 1 20ADD50; QUEUE ;2 

4 I 30;PP IOR ITY ;2 

4 140; TEST  G;02,5, ADD5I 

4 I 53;  PR  I OR  ITY?  3 

4 1 63ADD5I ;TEST  G; 02 , 1 0, ADD52 
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4 ! 7 0;PR IORITY »4 
4 180ADD52 ?GATE  SE;1I 


4 1 OQ  ; PREEMPT  ; RV.'Y  , PR 

4200; DEPART? 2 

421  3;  AD  VANCE  ? PH2,  Ff!  1 4 

4 223;  RETURN' ;R'..Y 

4 2 33; TRANSFER  *, ADD  I 30 

4240* 

4 25  U*"  I SSED  APPROACH 
426  '* 

427 3A DD60 ; LEAVE ♦ I I 

4 2S0ADD63 ;TEST  E ; 3L4, 220, ADDA  I 

4290 ; RETURN; BDF 10 

4300ADD6! ; RETURN ; BDF 1 I 

43 t 0; RETURN; BDF 12 

4320?  RETURN ;BDF  1 3 

4 333ADD64;UARfC?2PF 


4 J40 

MSA VE VALUE? 

4 35  0 

MSAVEVALUE; 

4 360 

"SAVE  VALUE? 

4 37  0 

"SAVEVALUE? 

4 380 

"SAVEVALUE; 

4 390 

"SAVEVALUE; 

4400 

"SAVEVALUE? 

44  l 0 

MSAVEVALUE? 

442  0 

SA  VE  VALUE?  1 

443  0 

a;  ARK ; 1 PF 

444  0 

TEST  NEIPL2 

44  5 -5 

ENTER? 6 

446  OF I VE? macro  ; 1 

447  0! 

ASSIGN; 15,  A 

I0+.PLI  ,1  , 

I0+,PL2,2, 

1 0+ , PL6 ,3,1 
I 0+ , DL7 , 4 , 1 
1 0+ , PL  1 8 , 5 , 

I 0+.3LI  5,6,  V$VAR5,"F 
1D+,PL2,B, 1 , "F 
10+,PL2,7, V$VAR5,NF 
4 + , VS  VAP5.XF 


,"F 
,"F 
,HF 
fF 
,"F 


444  U‘<AP;FVARI  ABLE;  1 8/?L4*3600 


440 0 ; -\D  VA NCE  ? VS  MAP , FN  1 4 
4.  J;F  REEMPT?BDF7,PR 
45 1 C ;TRAN SFER ? , IAP 
4 52 MA nD62 ; ASSIGN  t I 7 , 0 , P H 
4533; ASSIGN; 1 5,9, PL 
4540;  SPLIT*  1 ,ADD96,  1 7Pf: 

4 550? SA  DEVALUE? 9,  VSTI’.'l  ,XL 
4560; TRANSFER  ? , ADD30 
4573* 


4580*1  NSTPUNENT  APPROACHES 
4 590* 


4600ADD70 ; PR  I OR ITY ; 2 

46I0HEST  NE ;PL4 , 80 , ADD7 1 

4 62  3?  TRAN’ SFER  ?.  7R0 ,,  A CD  7 1 

4 6 3 0 VA  R 1 1 IFVARI ABLE ;PL23/PL3*3600 

4^4  jv  SI  2 IFVARI  ABLE*  ( (24/PL3)  + ( 1 5/PL4  ))*3< 

465  JVA  ? 1 3 * F VAR  I ABLE  * 6/P  L4*3  60! 

4 66  0;  ADVANCE  IVSVARI  1 , FN  1 5 
4*70; ENTER ; 4 

463 WF I V5? MACRO? 16,7,19,4,18,(0 
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469  0; ASSIGN;  1 5 , AC  1 , PF 
4700BDL4; F VAR  I ABLE ? I /PL4*36t 0*FN 1 4 
471  01  A SSI  GN  ? I 4 , VS5CL4  , PL 
4 / 2 0 ; A SS I GT  ; I 5 , V $ F I ' I A L , P I ' 

47341  PREEMPT! HAA, PR 
4 740  | ADVANCE?  v'S  VAR  I 2,  FN  15 
4 /50;  RETURN  5HAA 
4760?  ADVANCE?  VS  VAR  1 3,  FN  15 
4 770  ? PREEMPT  ? BDF7 , PR 
4 750? TRANSFER?, IAP 


a r'  ,-i  . 

U / y yr 

48O0APD71 ?ASSIGN?!8, II, PL 
4 8 I 7 V A R I 4 ? F VA  RI ABLE ? PL24/PL4*3oO0 
4 82 OVA R 15  ?F VARIABLE ? I 7/PL4*3600 


4 837  VAR  I 6 ? FVARI ABLE? 2 /PL5 *3600 

4840; ADVANCE ? Vs VAR  I 4, FN 1 5 

4850 ? ENTER ?5 

4860  ? A 5SIGN  ? 16, 7, PL 

4 870?  ASSIGN? 19,5, PL 

4 880? ASSIGN? I 5, AC!  , PF 

4 89  0BDL5?FVARIABLE;  1 /PL4*360O->FN1  4 

4900;ASSIGN;l4tV$SDL5,PL 

49 1 0 ? A SSI GN  ? 1 5 , VS  F I NA  L , PH 

4920 ? PREEMPT? LAA,  PR 

493  0?  AD  VANCE?  VS  VAR  I 5,  FM  15 

4940? RETURN? LAA 

495  0?  A D VANCE  ? VS  VA  R I 6,  FN  1 5 

496  0?  PREE"PT?  PDF7 , P R 
4970  ? TRANSFER  ? , IAP 
4980* 


4 990APD80 ? TEST  NE?PLI  = , ll,ADD8l 
5 000?  SA  VE  VALUE  ? 7, 0,  XL 
d 1 1 0 ? A SS  I Gf  ■ ? I 5 , 7 , PL 
5020?  SPLIT? 1 , ADD96 , 17PF 

DO3  0TI M2? FVARI A3LE ? VSTI V I -VSVAR I 2- VS VAR  I 3 

5 04 0 ? AS S I GN  ?30, VSTI  M2 , PL 

505 3? TEST  GE  ? PL30  , I , ADD 30 

5 060? TEST  LE? PL3C, 300 , ADD82 

5070?  SA  VEVALUE?  7, 300, XL 

508 OVA R I 7? FV API  ABLE ? 300-PL 30 


5090?  AD  VANCE?  VS  VAR  1 7 
5 1 OO  ? TRANSFER  ? , ACD30 
5 I I OA DD82 ? SAVE  VALUE  ? 7 , PL30 , XL 
5 1 20  ?TRANSFER  ? , ADD30 
5 I 3 0* 

5 I 40ADD8 1 ? SA  VE  VALUE  I 8 , 0 , XL 
5 I 50; ASSIGN? 1 5 , 8, 3L 
5 I 60; SPLIT? i ,APD96, I 7PH 

5 I / OT  I"  3 ? F V AP  I A RLE  ? VSTI"  I -V  S VA  R I 5 - VS  VA  R I * 
5 1 80?  ASSIGN ;3 I , VSTI "3 , PL 
5 I 9 0? TE  ST  GE  ? PL3 1 , I t ADD30 
92 00? TEST  LE?  PL3 I , 300,ADD83 
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i:i  'jsaveval'J  ,3 
5 33  V \ ? I > 1 FV A.~.  I L2  * 3«  -JL3  1 
2311  ADVATfC  :fV$VA  • l 3 
24  Jt  RA  ISFES1,/  503  5 
: 35  !A  ~003 ; S W37AL  J3;~  ,?L31  . HO 
~ 5 ' ; TH  A.  .'Sr  EH  ; , A r ’3  ’ 

527  J* 

p 9 '.yy-.j 

f 2 5 -v 

' • APOt  52tT.-5ANSFE3t.943,,  A HD  I X 
531  lADOt  3 : I Rl 
537  JFIVEt  a:  t5,3G,l6,2,  18,5 

533  '-173;  "AC..':  ; 17,7.2  , 7,21  ,3 

534  It  ASSIGN  ; 1 5,  VS  HI..'  AL,  PH 

535  1?  ASSIGN*  I 2, 2'-  , PH 

535  );TPAHSF33 : .5  , , ADD 5 3 

537  ; 1 PF 

533  ' ; E’JTE. . ; 7 

539-1  ;C‘J.2UH;  i 

54 '0;  GATE  53;  1 

54!  *;  2. ’73  : ; 1 '■) 

542  ';  EHT32;  1 1 

543-i;  7333  :PT;FCF1  3 i 

544  ■);  P3P  i.  T;  1 

545  -3’PTi-  'Fll  ,23 

545  i ; THANSFEr. ; ,A7~35 

547"* 

54?  l** : 1 3HEP  A LTITUOE 
540  ;-v 

55  ' 5A CO 35 ; 3 A7E  NUtHIGHI  , ADD  105 

551  1790;  'AC 30; HIGH  I 

552  '7:  );  'AC HO  : HICH2 


553 

;t 

;o 

• ■ > 

1 

in 

IG 

" 

3 

554 

“3T 

cj*  rr 

; P 

-r 

0 

• • - 

2 

' 

« 

37 

555 

r~* 

10 

77  • ’/ 

A.C 

cO 

T 

i 

3H4 

* 

in 

HI 

JT! 

2 

77  • 1 1 

:.) 

r 

0* 

» /• 

• 

' J T 

- i 

'>*  r 

■j:  ' 

^57 

.'7' 

12 

E3 ; " 

A 

•'  * 

T 

g ; 

•'4 

• • 

HI 

(J’’ 

=;r:  a 

j .) 

3 3 1 

• "> 

,0 

; ' I 

I 

u 

•4 

c. 

HI 

0’ : 

^ J>  * > L ...  .-Z  4 ../»  I v j 4 j j f . , I vj ; *r  /\ 

56  ■ JTHSSEl  'A5  ; ;■■■:  3H6,HIOH4 3 
55!  J A no  S3  ? :HTE  . ; 1 
562  TFOUiil  ~F7,"I  ?H4C 

5*'  2 )F0  ; <AC  ’0?FDFS,HIGH5 
>64GF0U  t )lFDF9,HIG*-i6 

5 6 5 ^ * 7 .3 A; i 5 F E ;,A 7*729 
5 45  J* 

567  5AP05  7 ; 3 375  . ; HIGH  1 
36  JTHSEc*  AC  1 0lHIGH4,HI  5H2 

56  ■ JT  IS  St  'AH  *HI  JH4A.HIGH3 

57  ' '7':  ? - I AC  - IHI  JH4B.HIGH4 

571  '7"  HEtMACf  0IHIGH4C.MI0K4A 

572  'T"H  53*  'AC,/'  »HIJV'5,HIGH43 


i 


5 7 3 T’  ■ i = E ; M ACE 0 ; H 1 .3 --'6 , F I GrA C 

574  3AD099»ENTH  ?t10 

575  3E0JE ; 'ACRO  ; F0F7 ,H IGH5 

576  ’FOUR ; ,(AC20  ; ;1DF3  ,H I'JH6 


5 77  );RREE 


:9 , PE 


578  1;A0VA:;C  = ?PU5 
5 79JJT«A:;SFE!? ; , ADD34 
581'1  H 

531  >N£X7  HIGHER  ALT  ITU  HE 
582;* 


58  3 'LA 

E'H 

1 85  ; 

gate 

HU; 

:!EXT 

1 , ADD 

55  4. T 

WO 

• ' 

VC 

\ J fit 

EXTI 

55  5 3T 

WO 

• f 

AC 

20  ;fl 

EX  72 

5 3607 

ViO 

• ■/ 

AC 

20  »N 

EXT3 

587'; 

T'” 
1 — 

57 

ET 

;?L4 

,220 

, ATP 

1 07 

553  57 

W 

ACRO 

;:lex 

T4/; 

T Y T*  ] 

5 89  '7 

» f 

ACT) 

;::ex 

75,  N 

EXT2 

59  ■ 17 

TEE 

'( 

ACRO 

;.'IEX 

76,:. 

t a r 3 

591  3 A 

DD 

l . • ; 

MULL 

592  !T 

•JO 

TT 

H 

ACRO 

;HIG 

H4C, 

' I EXT 4 

59  387 

■j  i 

M 

ACRO 

*HIG 

H5,U 

EXT5 

594  JT 

HE 

EE 

y 

AC  RO 

; H I o 

H6,F 

EXT6 

E2 : , ADD33 

:'ETUHN;MEXT! 

acro  EX  T4,:; 

EX  72 

AC ?:.)  ;NEXT5,.. 

EXT  3 

HULL 

' CRO ; NEXT6 , N 

EX  T 4 

ACRO ; H IGH4C , 

ME  XT 

6 7 ' A CO  1 


6C30T-R  EE ; MACRO  ? H IOH5  .TEXTS 
6 84 OTHER EH;’1  AC EO  ? V I GH 6 , H I GH4C 
60501 T3  A1ISFZE 1 » A D039 
6 ; 6 H' 

6-'770*LAST  ALTITUDE 
603-7* 

689  lAnDl  1 5? MULL 

6 1 0 3T <VO  ; MAC  EO  ; LAST  I 

61  1 TWO  ; MACED  ;LAS72 

6 1 2 UTWO ; MACRO  ; LAST3 

61 37; TEST  E;?L4,220,ADD1 1 7 

6 1 4 THREE ; MACRO  I . . EXT* , L ASTI 

6 1 507HE EE ; V AC  DO f K EXT5 , L A5T2 

6 1 6 3TH  2 EE ; ’ ( AC  RO  ; : I EXT6 , L AST 3 

61 701 T3 A ISFE2 1.AC Dl 06 

6 1 3 '* 

6 1 9 J A CD  I I 7?  IETU  H. ; L AST  1 
2 3 THREE*  LACRO INZXT4  f LAST2 

621  XT'  TR EE ; MACRO  j:;EXT5 , LAST3 

622  ?1  TEA. FIFE:  ; , ADD!  09 

623  ':  LIST 

624  M START;  I 'TO 
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